Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,442,021 T→C 11.7% I480I (ATT→ATC thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,442,0210T→C11.7% 73.0 ‑3.0 ‑0.2247I480I (ATT→ATCthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (19/10):  ref base (99/119):  total (118/129)
Fisher's exact test strand distribution p-value = 4.87e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.61e-01

Only 100 of 264 total aligned reads displayed.

TGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGC  >  NC_007779/3441921‑3442121
                                                                                                    |                                                                                                     
tGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAtt                                                                                                       <  1:13289212/101‑1 (MQ=255)
tGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAtt                                                                                                       <  1:520305/101‑1 (MQ=255)
 gACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTg                                                                                                      <  1:6429762/101‑1 (MQ=255)
 gACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTg                                                                                                      <  1:1791722/101‑1 (MQ=255)
  aCCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTgg                                                                                                     >  1:11260956/1‑101 (MQ=255)
   ccGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGt                                                                                                    >  1:2376257/1‑101 (MQ=255)
    cGAGGAGTTAGAGCACTGCCACGAAGCGCCTTTTTACACTCTGGGGCCGCTAACTGCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                                                                                   <  1:13462738/101‑1 (MQ=255)
    cGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                                                                                   >  1:1669206/1‑101 (MQ=255)
            tAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatg                                                                                           <  1:4076703/101‑1 (MQ=255)
             agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                                                                          >  1:10972058/1‑101 (MQ=255)
              gagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatgat                                                                                         <  1:4577631/101‑1 (MQ=255)
                gCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTg                                                                                       >  1:1407081/1‑101 (MQ=255)
                 tggtgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgc                                                                                      >  1:8341683/10‑100 (MQ=255)
                 cACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTgg                                                                                     >  1:10226117/1‑101 (MQ=255)
                  aCTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGc                                                                                    >  1:11314013/1‑101 (MQ=255)
                   cTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCt                                                                                   <  1:634202/101‑1 (MQ=255)
                     gCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTgg                                                                                 <  1:8626153/101‑1 (MQ=255)
                      tgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcagg                                                                                 >  1:553835/5‑95 (MQ=255)
                       cACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGtt                                                                               >  1:3876346/1‑101 (MQ=255)
                       cACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGtt                                                                               >  1:890446/1‑101 (MQ=255)
                         ggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctg                                                                              >  1:7270025/2‑92 (MQ=255)
                          gAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctgg                                                                             >  1:13708744/1‑91 (MQ=255)
                          gAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTgg                                                                            <  1:7011556/101‑1 (MQ=255)
                          gAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTgg                                                                            <  1:13821639/101‑1 (MQ=255)
                           aaGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGc                                                                           >  1:417474/1‑101 (MQ=255)
                            aGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCt                                                                          >  1:2984189/1‑101 (MQ=255)
                              cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                        >  1:6749588/1‑101 (MQ=255)
                              cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                        <  1:7117127/101‑1 (MQ=255)
                                ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCCCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                                                      <  1:11992107/101‑1 (MQ=255)
                                ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                                                      <  1:2604202/101‑1 (MQ=255)
                                 cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                     >  1:9025261/1‑101 (MQ=255)
                                 cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                     <  1:10712738/101‑1 (MQ=255)
                                 cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                     >  1:5061392/1‑101 (MQ=255)
                                 cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                     >  1:8278608/1‑101 (MQ=255)
                                   ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                                   <  1:9631677/101‑1 (MQ=255)
                                   ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                                   <  1:8412876/101‑1 (MQ=255)
                                    tttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATg                                                                  >  1:1499119/1‑101 (MQ=255)
                                      tACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGct                                                                <  1:79210/101‑1 (MQ=255)
                                       acacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTAGCTGTGCGactgtc                                                               >  1:4954278/1‑96 (MQ=255)
                                        cacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctct                                                              >  1:1325540/1‑101 (MQ=255)
                                         acTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTg                                                             <  1:10364610/101‑1 (MQ=255)
                                           tctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGtt                                                           <  1:3270501/101‑1 (MQ=255)
                                              ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                        <  1:13936531/101‑1 (MQ=255)
                                              agggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                        >  1:10058201/2‑101 (MQ=255)
                                                ggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTa                                                      >  1:5694349/1‑101 (MQ=255)
                                                   cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                                   <  1:9736532/101‑1 (MQ=255)
                                                   cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAACg                                                   >  1:7682134/1‑101 (MQ=255)
                                                    gCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGc                                                  <  1:573080/101‑1 (MQ=255)
                                                     cTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGcc                                                 <  1:760251/101‑1 (MQ=255)
                                                      tggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGccg                                                >  1:8968330/6‑100 (MQ=255)
                                                      tAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTTGCTGCGCGATGCTCTGTTACGTAACGCCa                                                >  1:5147623/1‑101 (MQ=255)
                                                       aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                               >  1:10644174/1‑101 (MQ=255)
                                                         cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGAGTGGCTGGTTTGGGTGCGCGATGCTCTGTTACGTAACGCCaaaa                                             <  1:13033198/101‑1 (MQ=255)
                                                          tACCGATATTGCGCCGGGCTATGTCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                            <  1:8721583/101‑1 (MQ=255)
                                                          tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                            <  1:7199490/101‑1 (MQ=255)
                                                          tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                            <  1:7395214/101‑1 (MQ=255)
                                                          tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                            <  1:2806844/101‑1 (MQ=255)
                                                          aaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCgaagg                                            >  1:12275970/2‑99 (MQ=255)
                                                           aCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAaga                                           <  1:765434/101‑1 (MQ=255)
                                                            ccGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAggag                                          >  1:474629/1‑101 (MQ=255)
                                                                tatTTCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTTCGGCGATGATTGGCTGGTTTTGCTGCGCAATGGCCCGCTACGCAATGCCAAAAGAGCaac                                      >  1:13649561/1‑99 (MQ=255)
                                                                 atcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCt                                     <  1:5969454/98‑1 (MQ=255)
                                                                 atTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCt                                     <  1:4610107/101‑1 (MQ=255)
                                                                  ttGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTg                                    <  1:13141453/101‑1 (MQ=255)
                                                                  ttGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTg                                    <  1:12898592/101‑1 (MQ=255)
                                                                   tGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCTGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg                                   <  1:3691954/101‑1 (MQ=255)
                                                                   tGCGCCGGGCTATGACCACTTCACGTCGGCGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg                                   >  1:8142239/1‑101 (MQ=255)
                                                                     cgcCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGt                                 <  1:11730317/101‑1 (MQ=255)
                                                                        cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGGTGCTCTGTTACGCAACGACAAAAGAGCATCTgggagtg                              >  1:9298916/1‑97 (MQ=255)
                                                                        cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                              >  1:11791363/1‑101 (MQ=255)
                                                                        cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                              >  1:8530109/1‑101 (MQ=255)
                                                                           gCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                           <  1:5421140/101‑1 (MQ=255)
                                                                           gCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                           >  1:12103995/1‑101 (MQ=255)
                                                                              aTGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTTGGTCTGCCCAAt                        >  1:4246650/1‑101 (MQ=255)
                                                                              aTGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAAt                        >  1:12506214/1‑101 (MQ=255)
                                                                               tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                       <  1:1049386/101‑1 (MQ=255)
                                                                                gACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa                      >  1:2204692/1‑101 (MQ=255)
                                                                                  ccACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCAAAAAGAGCATCTGGGTCTGCCCAATAaag                    <  1:12191374/101‑1 (MQ=255)
                                                                                   cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga                   <  1:13332643/101‑1 (MQ=255)
                                                                                   cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga                   <  1:53333/101‑1 (MQ=255)
                                                                                   cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga                   >  1:7730998/1‑101 (MQ=255)
                                                                                      ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGTCCAATAAagaaga                <  1:1947835/101‑1 (MQ=255)
                                                                                      ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                >  1:528132/1‑101 (MQ=255)
                                                                                       tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt               <  1:13559411/101‑1 (MQ=255)
                                                                                        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg              >  1:3347185/1‑101 (MQ=255)
                                                                                         aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt             <  1:5627652/101‑1 (MQ=255)
                                                                                           gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa           <  1:6704658/101‑1 (MQ=255)
                                                                                           gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa           <  1:352040/101‑1 (MQ=255)
                                                                                            tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTaa          <  1:5795671/101‑1 (MQ=255)
                                                                                            tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTaa          <  1:4894091/101‑1 (MQ=255)
                                                                                             cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg         >  1:802055/1‑101 (MQ=255)
                                                                                             cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg         >  1:9149052/1‑101 (MQ=255)
                                                                                              ggggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc        <  1:9592434/101‑1 (MQ=255)
                                                                                               gggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa       >  1:12786883/1‑101 (MQ=255)
                                                                                               gggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa       <  1:2871870/101‑1 (MQ=255)
                                                                                                 gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgg     <  1:5165789/101‑1 (MQ=255)
                                                                                                  aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAggg    >  1:12138044/1‑101 (MQ=255)
                                                                                                   ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg   >  1:3135513/1‑101 (MQ=255)
                                                                                                   ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg   <  1:10026019/101‑1 (MQ=255)
                                                                                                    tGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGc  >  1:1585312/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
TGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGC  >  NC_007779/3441921‑3442121

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.