Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,586 G→C 6.2% intergenic (‑51/+18) guaA ← / ← guaB GMP synthetase (glutamine aminotransferase)/IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,5860GC 6.2% 425.5 / 16.0 142intergenic (‑51/+18)guaA/guaBGMP synthetase (glutamine aminotransferase)/IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (62/71);  new base C (4/5);  total (66/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

Only 100 of 144 total aligned reads displayed.

GTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCC  >  NC_000913/2632437‑2632735
                                                                                                                                                     |                                                                                                                                                     
gTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCCTTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATg                                                                                                                                                       <  2:644280/150‑1 (MQ=255)
     cacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTc                                                                                                                                                  >  2:1718113/1‑150 (MQ=255)
     cacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTc                                                                                                                                                  >  1:1787743/1‑150 (MQ=255)
      acCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCg                                                                                                                                                 >  1:1679115/1‑150 (MQ=255)
        ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCggg                                                                                                                                               <  1:954058/150‑1 (MQ=255)
         cctcTCCCGCCCGCGGCGCGCAACCAGGTGCGTGTCCTGAGACCCGAAGTCCAGAATGAGGCTGCGAGGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGAGAAACAGATAATATAAATCGCCCGACATGAAGGCGGGc                                                                                                                                              <  2:120839/147‑1 (MQ=255)
         cAGCTCACGCCGTGGGTGCGGAACCAGTTGAGTGCACTGAGAACCTAAGTCAAGTATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGTAACAGATAATATAAATCGCCCGACATGAAGTCGGGc                                                                                                                                              <  2:1511280/150‑1 (MQ=255)
         cAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGc                                                                                                                                              <  1:1639877/150‑1 (MQ=255)
         cAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGAGGCGATGCTTATGAATTTTTTCCGTCATAGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                                                                              <  2:1825488/150‑1 (MQ=255)
         cAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACGGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCATATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                                                                              <  1:361574/150‑1 (MQ=255)
          cgCTCACGCACGCGGCGCGCAACCCGTTGAGTGTACTGCGCACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATCTAAATCGCCCGACTTCATGTCGGGCg                                                                                                                                             <  1:1540646/149‑1 (MQ=255)
           gCTCAAGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGa                                                                                                                                            <  1:211105/150‑1 (MQ=255)
            cTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa                                                                                                                                           >  1:1293368/1‑150 (MQ=255)
                 gcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGCCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCAGGCGAAGAGaa                                                                                                                                      >  2:1426042/1‑150 (MQ=255)
                   acgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATc                                                                                                                                    <  1:404346/150‑1 (MQ=255)
                    cgcgGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCa                                                                                                                                   >  1:1636722/1‑150 (MQ=255)
                    cgcgGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCa                                                                                                                                   >  1:513799/1‑150 (MQ=255)
                      cggcgCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAgg                                                                                                                                 <  1:543178/150‑1 (MQ=255)
                        gcgcgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAg                                                                                                                               <  1:142020/150‑1 (MQ=255)
                            gcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCa                                                                                                                           >  2:461709/1‑150 (MQ=255)
                            gcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCa                                                                                                                           >  1:466683/1‑150 (MQ=255)
                                 cAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACgg                                                                                                                      >  1:697946/1‑150 (MQ=255)
                                 cAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACgg                                                                                                                      >  1:1517781/1‑150 (MQ=255)
                                   gTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTa                                                                                                                    >  1:1251406/1‑150 (MQ=255)
                                      gAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGGCATTGACGCTTATTCCGAGTCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGtt                                                                                                                 <  2:1800610/150‑1 (MQ=255)
                                        gtgtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAACATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCg                                                                                                               >  2:461723/1‑150 (MQ=255)
                                          ctACTTAGAACCGAACTCCAGAATGAGTCTGGTATGCATATGAATGTTTTCCGTCATTGACGCTTATTCAGAGGCAATTGAAACAGATAATATAAATCGCCCTACATGAAGTCTGGCGAAGAGAATAAGGAGCCCAGACGGTAGTTCggg                                                                                                             <  2:1174411/149‑1 (MQ=255)
                                             cTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                                                                          >  1:1355301/1‑150 (MQ=255)
                                             cTGAGAACCGAAGTCCAGAATGAGGATGAGATGCTTATGAACGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                                                                          <  1:1755515/150‑1 (MQ=255)
                                               gagaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGCCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGActc                                                                                                        >  1:1383339/1‑150 (MQ=255)
                                                 gaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATCAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCtt                                                                                                      <  1:17457/150‑1 (MQ=255)
                                                 gaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCtt                                                                                                      >  2:1036404/1‑150 (MQ=255)
                                                  aaCCGAAGTCCAGAATGCGGATGCGATGCTTATGAATGTTTTCCGTAATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCttt                                                                                                     <  2:717853/150‑1 (MQ=255)
                                                    ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                                                                                   <  2:1518758/150‑1 (MQ=255)
                                                    ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGAATAATCCGAGGCAAATCAAACAGATTATATAAATCGCCCGACTTGAAGTCGGGCGAAGAGAATCAGGAGCCAAGACGGTAGGTCGGGGACTCTTTAg                                                                                                   >  1:578302/1‑150 (MQ=255)
                                                      gAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTa                                                                                                 >  1:1797255/1‑150 (MQ=255)
                                                      gAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTa                                                                                                 <  1:1561006/150‑1 (MQ=255)
                                                      gAAGCCCAGAATGAGGATGCGATGCTATAGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATCTAAAGCGCCCGACAGGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTa                                                                                                 <  2:314472/150‑1 (MQ=255)
                                                         gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                              <  2:1293368/150‑1 (MQ=255)
                                                          tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATgg                                                                                             <  2:554868/150‑1 (MQ=255)
                                                          tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATgg                                                                                             <  1:1006055/150‑1 (MQ=255)
                                                           ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGt                                                                                            >  2:1779631/1‑150 (MQ=255)
                                                           ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGt                                                                                            >  1:1658258/1‑150 (MQ=255)
                                                            cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                                                                           >  1:798951/1‑150 (MQ=255)
                                                            cAGAATGAGGATGCGATGCTTATGAATGTGTTCCGGCATTGACGCTTCTTCCGAGCCACGTGAAACAGATAATCTAAATCGCCCGACATGAAGTCGGGCGAAGAGTATCAGGAGCCCCGACGGTAGTTCGGGGACTCTTTAGTACTggcc                                                                                           <  2:771553/150‑3 (MQ=255)
                                                                    ggATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGTTCACGTCGTg                                                                                   >  2:1253828/1‑150 (MQ=255)
                                                                    ggATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                                                                                   <  2:649232/150‑1 (MQ=255)
                                                                    ggATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                                                                                   <  1:461709/150‑1 (MQ=255)
                                                                    ggATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGGCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                                                                                   <  2:1791339/150‑1 (MQ=255)
                                                                       tgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAAc                                                                                >  1:1025726/1‑150 (MQ=255)
                                                                        gcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACg                                                                               >  2:199426/1‑150 (MQ=255)
                                                                          tatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTg                                                                             <  2:825534/149‑1 (MQ=255)
                                                                           atgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTgg                                                                            <  1:1426632/150‑1 (MQ=255)
                                                                            tgcTTATGATTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGc                                                                           <  2:1765456/150‑1 (MQ=255)
                                                                                  tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGTGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                                                                     >  1:48897/1‑150 (MQ=255)
                                                                                  tGAATGTTTTCCGTCATAGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                                                                     <  2:311949/150‑1 (MQ=255)
                                                                                   gAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTg                                                                    <  1:1035771/150‑1 (MQ=255)
                                                                                     aTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa                                                                  <  1:1134209/150‑1 (MQ=255)
                                                                                         tttCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc                                                              <  2:1439804/150‑1 (MQ=255)
                                                                                           tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg                                                            >  2:1427424/1‑150 (MQ=255)
                                                                                           tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg                                                            <  1:1292530/150‑1 (MQ=255)
                                                                                            ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGc                                                           <  2:1256816/150‑1 (MQ=255)
                                                                                            ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGc                                                           >  1:1121731/1‑150 (MQ=255)
                                                                                             cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCa                                                          >  1:1449383/1‑150 (MQ=255)
                                                                                                  ttGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGACGGCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt                                                     <  1:593676/150‑1 (MQ=255)
                                                                                                  ttGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt                                                     <  2:1355301/150‑1 (MQ=255)
                                                                                                   tGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTg                                                    <  2:1111910/150‑1 (MQ=255)
                                                                                                    gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGa                                                   >  2:1364032/1‑150 (MQ=255)
                                                                                                       gCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATAc                                                <  2:498454/150‑1 (MQ=255)
                                                                                                         ttATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGTGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGt                                              <  1:76127/150‑1 (MQ=255)
                                                                                                          tATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTa                                             <  1:203288/150‑1 (MQ=255)
                                                                                                            ttCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCATGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACa                                           <  1:1527866/150‑1 (MQ=255)
                                                                                                              ccGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACaaa                                         >  1:1259466/1‑150 (MQ=255)
                                                                                                                  ggCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTAAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTcc                                     >  1:1632966/1‑150 (MQ=255)
                                                                                                                  ggCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCACACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTcc                                     <  2:256517/150‑1 (MQ=255)
                                                                                                                    cAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGc                                   <  2:1514682/150‑1 (MQ=255)
                                                                                                                    aaaGTGAAACATATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGc                                   <  2:215842/149‑1 (MQ=255)
                                                                                                                     aacttaaaCAGTTAATATTAATCGCCAGACATGAATTCGGGCGAAGAGAATCTGGATCCCAGACGGTAGATCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCGTGAATGCCCGCACCGCTGATACTTACAAACTCCGCt                                  <  2:1679115/145‑1 (MQ=255)
                                                                                                                       gTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGCTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCttt                                >  2:1333147/1‑150 (MQ=255)
                                                                                                                        tGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTa                               <  2:615257/150‑1 (MQ=255)
                                                                                                                         gAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAg                              >  1:587070/1‑150 (MQ=255)
                                                                                                                          aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACTGTAGTTCGGGGACTCTTTAGTAATTGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGt                             <  1:1208376/150‑1 (MQ=255)
                                                                                                                          aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGt                             >  2:1666848/1‑150 (MQ=255)
                                                                                                                                 aaataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGt                      <  1:1451092/149‑1 (MQ=255)
                                                                                                                                  aataTAAATCGCACGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGtt                     <  1:1718113/150‑1 (MQ=255)
                                                                                                                                   atatAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCGTTAGGAATGGTCACGTCGTGCACGTGGCTTTACTGAATGCCCGCACCGCTGATAGGTACAAAATCCGCGGTAGTAAGccggtg                    >  1:1517922/1‑145 (MQ=255)
                                                                                                                                    tataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACTTACAAACTCCGCTTTAGTACGCAGTtcg                   >  2:689794/1‑150 (MQ=255)
                                                                                                                                      taAATCGCCCGAAATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtcgtc                 <  1:1696938/150‑1 (MQ=255)
                                                                                                                                       tcttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCAGCTTTAGTACGCAGTttgtcg                >  1:26703/4‑150 (MQ=255)
                                                                                                                                       aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtcgtcg                <  1:446365/150‑1 (MQ=255)
                                                                                                                                         ttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGAt              >  2:1609248/2‑150 (MQ=255)
                                                                                                                                           cGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAg            <  2:159437/150‑1 (MQ=255)
                                                                                                                                            gCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGt           >  1:477671/1‑150 (MQ=255)
                                                                                                                                             cccGACATGAAGTCGGGCGAAGAGAATCAGCAGCCCAGACGTTAGTTCGGGGACTCGTTAGTAATGTTCACCTCGTGAGCGAGGCTTTCCTGAATGCCAGCACCGCTGAAACGTACAAACTCCGCTTTAGTACGCAGTTCGACGATAgga          >  1:767580/1‑148 (MQ=255)
                                                                                                                                                  cATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACcaca     <  1:461723/150‑1 (MQ=255)
                                                                                                                                                  cATGAAGTCGGGCGAAGAGAATCAGGAGCACAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTTAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTACTACGCAGTGCGTCGATAGTACcacc     >  2:1559545/1‑149 (MQ=255)
                                                                                                                                                   aTGACGTCGGGCGAAGAGACTCAGGAGCCCAGCCGGTCGGTCGGGGACTCTATAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGACTGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAg    <  2:246396/150‑1 (MQ=255)
                                                                                                                                                    tGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCCCCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGc   >  1:52901/1‑150 (MQ=255)
                                                                                                                                                    tGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGc   >  1:434537/1‑150 (MQ=255)
                                                                                                                                                     gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGcc  <  1:705916/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCC  >  NC_000913/2632437‑2632735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: