Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,442,017 G→C 7.6% G479A (GGG→GCG)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,442,0170G→C 7.6% 49.9 ‑2.9 ‑0.5249G479A (GGG→GCG) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (14/5):  ref base (113/116):  total (127/122)
Fisher's exact test strand distribution p-value = 5.51e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.45e-01

Only 100 of 263 total aligned reads displayed.

AATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAG  >  NC_007779/3441917‑3442117
                                                                                                    |                                                                                                     
aaTATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCggg                                                                                                       <  1:7839392/101‑1 (MQ=255)
aaTATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCggg                                                                                                       >  1:13601237/1‑101 (MQ=255)
  tatGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGa                                                                                                     >  1:1771545/1‑101 (MQ=255)
      aCCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTgg                                                                                                 <  1:577996/101‑1 (MQ=255)
       ccGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGt                                                                                                <  1:10526054/101‑1 (MQ=255)
          aggagTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcg                                                                                             <  1:5740170/101‑1 (MQ=255)
          aggagTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcg                                                                                             >  1:6431114/1‑101 (MQ=255)
           ggagTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcgg                                                                                            <  1:7740261/101‑1 (MQ=255)
             agTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcggcg                                                                                          <  1:2733951/101‑1 (MQ=255)
             agTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcggcg                                                                                          <  1:4639894/101‑1 (MQ=255)
             agTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcggcg                                                                                          >  1:5669508/1‑101 (MQ=255)
              gTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCga                                                                                         >  1:765730/1‑101 (MQ=255)
              gTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCga                                                                                         >  1:12049183/1‑101 (MQ=255)
               ttAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgat                                                                                        >  1:11443268/1‑101 (MQ=255)
                tAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatg                                                                                       <  1:13844153/101‑1 (MQ=255)
                tAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatg                                                                                       >  1:8811194/1‑101 (MQ=255)
                tAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatg                                                                                       >  1:10805477/1‑101 (MQ=255)
                 agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                                                                      >  1:13737532/1‑101 (MQ=255)
                    gtggtgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTg                                                                                   >  1:3397699/11‑101 (MQ=255)
                    gCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTg                                                                                   <  1:5588192/101‑1 (MQ=255)
                      ggtgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcg                                                                                 <  1:14022737/93‑3 (MQ=255)
                      aCTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGc                                                                                >  1:11303355/1‑101 (MQ=255)
                       cTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCt                                                                               >  1:3016062/1‑101 (MQ=255)
                       cTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCt                                                                               >  1:3550012/1‑101 (MQ=255)
                        tGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATAATGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTg                                                                              >  1:9854712/1‑101 (MQ=255)
                            aCGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGttt                                                                          >  1:10302757/1‑101 (MQ=255)
                            aCGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGgtt                                                                          >  1:7877899/1‑101 (MQ=255)
                             cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTg                                                                         >  1:3737109/1‑101 (MQ=255)
                             cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATGGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                                                         <  1:4235302/101‑1 (MQ=255)
                              gAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctgg                                                                         <  1:8335209/101‑11 (MQ=255)
                              gAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTgg                                                                        <  1:10173848/101‑1 (MQ=255)
                               aaGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGc                                                                       >  1:725740/1‑101 (MQ=255)
                               aaGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGc                                                                       >  1:898013/1‑101 (MQ=255)
                               aaGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGc                                                                       >  1:1342785/1‑101 (MQ=255)
                                  cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTTGCTgc                                                                    >  1:12687009/1‑101 (MQ=255)
                                  cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                    >  1:4590960/1‑101 (MQ=255)
                                  cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                    <  1:2531938/101‑1 (MQ=255)
                                  cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                    <  1:8397261/101‑1 (MQ=255)
                                   gcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcg                                                                   <  1:7326077/101‑1 (MQ=255)
                                    ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                                                  <  1:9810357/101‑1 (MQ=255)
                                    ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                                                  <  1:1074944/101‑1 (MQ=255)
                                     cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                 <  1:6130310/101‑1 (MQ=255)
                                     cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                 >  1:6275881/1‑101 (MQ=255)
                                      gTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGa                                                                >  1:9936514/1‑101 (MQ=255)
                                      gTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGa                                                                >  1:14291294/1‑101 (MQ=255)
                                      gTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGCTTGGCTGCCCGa                                                                >  1:11743051/1‑101 (MQ=255)
                                       ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                               >  1:8163898/1‑101 (MQ=255)
                                          tACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGct                                                            <  1:7868142/101‑1 (MQ=255)
                                          tACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGct                                                            <  1:4711840/101‑1 (MQ=255)
                                            cacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctct                                                          <  1:9544552/101‑1 (MQ=255)
                                             acTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTg                                                         >  1:12269720/1‑101 (MQ=255)
                                             acTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTg                                                         <  1:7159086/101‑1 (MQ=255)
                                                  ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                    <  1:9479649/101‑1 (MQ=255)
                                                   gggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGt                                                   >  1:2141669/1‑101 (MQ=255)
                                                     gCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTaa                                                 <  1:6623327/101‑1 (MQ=255)
                                                     gCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTaa                                                 <  1:10001353/101‑1 (MQ=255)
                                                     gCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTaa                                                 <  1:1934415/101‑1 (MQ=255)
                                                      ccGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAAc                                                >  1:494737/1‑101 (MQ=255)
                                                       cgctggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACg                                               >  1:10686588/9‑101 (MQ=255)
                                                       cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                               <  1:8859499/101‑1 (MQ=255)
                                                        gctggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGc                                              >  1:12858266/8‑101 (MQ=255)
                                                        gctggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGc                                              >  1:3999553/8‑101 (MQ=255)
                                                        gCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGc                                              <  1:6222893/101‑1 (MQ=255)
                                                          tAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCa                                            >  1:11650273/1‑101 (MQ=255)
                                                           aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                           >  1:7592619/1‑101 (MQ=255)
                                                           aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                           >  1:13191689/1‑101 (MQ=255)
                                                            aCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaa                                          <  1:138767/101‑1 (MQ=255)
                                                              tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                        <  1:290236/101‑1 (MQ=255)
                                                              aaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCgaagg                                        >  1:6639615/2‑99 (MQ=255)
                                                               aCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAaga                                       <  1:4808717/101‑1 (MQ=255)
                                                                ccGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                      <  1:1575776/101‑1 (MQ=255)
                                                                 cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTTGCTGCGCGATGCTCTGTTACGTAACGCCAAAAAAGc                                     >  1:4320469/1‑101 (MQ=255)
                                                                 cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc                                     >  1:1355064/1‑101 (MQ=255)
                                                                 cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc                                     >  1:5314540/1‑101 (MQ=255)
                                                                 cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc                                     >  1:744740/1‑101 (MQ=255)
                                                                   atatTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCAt                                   >  1:9158479/1‑101 (MQ=255)
                                                                     atTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCt                                 >  1:10349323/1‑101 (MQ=255)
                                                                      ttGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTg                                <  1:4047166/101‑1 (MQ=255)
                                                                      ttGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTg                                <  1:2376843/101‑1 (MQ=255)
                                                                       tGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg                               >  1:701161/1‑101 (MQ=255)
                                                                            cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                          >  1:10150298/1‑101 (MQ=255)
                                                                            cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                          <  1:3418703/101‑1 (MQ=255)
                                                                            cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                          >  1:2771767/1‑101 (MQ=255)
                                                                            cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                          <  1:6820897/101‑1 (MQ=255)
                                                                             gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                         >  1:5170825/1‑101 (MQ=255)
                                                                                cTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCa                      <  1:1514057/101‑1 (MQ=255)
                                                                                cTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCa                      <  1:1248712/101‑1 (MQ=255)
                                                                                   tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                   <  1:8006877/101‑1 (MQ=255)
                                                                                   tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                   <  1:13091576/101‑1 (MQ=255)
                                                                                    gACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa                  >  1:5016839/1‑101 (MQ=255)
                                                                                       cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga               <  1:9689648/101‑1 (MQ=255)
                                                                                          ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGGTGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga            >  1:9105133/1‑101 (MQ=255)
                                                                                          ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga            <  1:446320/101‑1 (MQ=255)
                                                                                           tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt           <  1:14225772/101‑1 (MQ=255)
                                                                                              cGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt        >  1:5522945/1‑101 (MQ=255)
                                                                                               gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa       <  1:11139370/101‑1 (MQ=255)
                                                                                                tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTaa      >  1:10921996/1‑101 (MQ=255)
                                                                                                 cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg     <  1:9226665/101‑1 (MQ=255)
                                                                                                   gggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa   >  1:1156249/1‑101 (MQ=255)
                                                                                                    ggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAg  <  1:3903893/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
AATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAG  >  NC_007779/3441917‑3442117

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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.