Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,487,137 G→A 10.6% intergenic (+73/‑74) purH → / → purD bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,487,1370GA10.6% 302.3 / 20.2 115intergenic (+73/‑74)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (43/56);  new base A (3/8);  total (49/66)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 134 total aligned reads displayed.

TGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCG  >  NC_002947/5486989‑5487286
                                                                                                                                                    |                                                                                                                                                     
 gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                                     <  2:2324072/150‑1 (MQ=255)
 gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                                     <  1:3408122/150‑1 (MQ=255)
 gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                                     <  2:2324244/150‑1 (MQ=255)
 gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                                     <  2:2324281/150‑1 (MQ=255)
 gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                                     <  2:3030842/150‑1 (MQ=255)
   tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATAAGGTGGTGGATTAACcgccgc                                                                                                                                                   >  2:1454071/1‑150 (MQ=255)
   tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGGATACGGGGGTGGGTTGAAcgccgc                                                                                                                                                   >  2:2430910/1‑150 (MQ=255)
         cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                             <  2:1297014/150‑1 (MQ=37)
          aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCATAAGATTGTGGTTTCACcgccgccgtcccg                                                                                                                                            >  1:1843923/1‑143 (MQ=255)
          aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                            <  1:718999/150‑1 (MQ=255)
          aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                            <  2:3254426/150‑1 (MQ=255)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                           <  2:302141/150‑1 (MQ=37)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                           <  2:1223951/150‑1 (MQ=37)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                           <  2:3978050/150‑1 (MQ=37)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                           <  1:2149980/150‑1 (MQ=255)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                           <  1:1883869/150‑1 (MQ=255)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                           <  1:1815334/150‑1 (MQ=255)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                           <  1:3396476/150‑1 (MQ=255)
           tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                           <  1:2503146/150‑1 (MQ=25)
            cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCggg                                                                                                                                          <  2:903717/150‑1 (MQ=37)
            cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                                          <  2:760849/150‑1 (MQ=255)
                gctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCCCGCAAAAGAGGGTGGTTTAAACCCCCCCGTCGCGGGCacg                                                                                                                                      >  1:3002191/1‑148 (MQ=255)
                 ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc                                                                                                                                     <  2:3049898/150‑1 (MQ=25)
                 ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                                     <  1:1246499/150‑1 (MQ=255)
                 ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                                     <  2:3133775/150‑1 (MQ=255)
                    ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATTCGATGGTGGAATCACCCCCGCAATCGCGGGCAAGCCCg                                                                                                                                  >  2:3747229/1‑150 (MQ=255)
                              cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTTGCCGCGTATAAGATTGTTGTTTTAACGCCGGATTGTCGGGGAAGACCGCTACCACAgt                                                                                                                        >  2:2421173/1‑150 (MQ=255)
                                     gcgTTGGTCTCCACTGGCTTGCGCCACTCCCGCCACTAATTACGCGGATCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgc                                                                                                                 <  2:1070934/150‑1 (MQ=37)
                                       gagggTCTTCACTGGCATGCGCCACTCCCGCCACTACTTACGCGGATCCCGAGGCAAGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgcgc                                                                                                               <  2:3348367/147‑1 (MQ=255)
                                             cTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCACCCCCCGTGGTCGGCGCAAGccc                                                                                                         >  2:4140900/1‑149 (MQ=255)
                                                       aTGCGCCACTTCCGCCCCTAATTCCGCGGATCCCGCGGCACGTCGAGACCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGCATCCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACt                                                                                               <  2:4460459/150‑1 (MQ=37)
                                                        ggcgcCCCTCCCGCCACTACTTCCGGGGATCCGGCGGCACGTCGCGCTCCGGTGGGCGCGGGCTTGCCGGCTAATCCGATGGTGTATTCCCGGCCGCCTGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                                              <  1:4119238/149‑1 (MQ=255)
                                                               ctgccGCCACTAATCCCGGGGATCCGGCGGCACGGCGATATCCTGTGGGTGGGGGCTTGCCCGCCAATCCGTGGGTGGATACCCCGCCGCTTGGGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCGCTTGAATCGa                                                                                       <  2:2264030/147‑1 (MQ=255)
                                                               cTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCCCATAAGATTGTGGGTTAACCGCCCCCCTACCCGGCCAGCCCCCGCCCCCCGGGTTGGGCGCAAGCCACGGGGGCACATGaaaaga                                                                                       >  2:2749636/1‑146 (MQ=255)
                                                                ttCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGTATTCACCGCCGCATTGGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                      <  2:3282484/150‑1 (MQ=255)
                                                                          aaTTCCGCGGATCCCGCGGCACGTCGAGATCGGGTGGGTGCGGGCTGGCCGGCTAATCCGTGGGTGGATCCCCGGCCGCTTTGGCGGGCATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAt                                                                            <  2:3894091/150‑1 (MQ=255)
                                                                             tccgggtaTCCCGTGGCATGTCGCGATCGTGTGGGCGGGGGCTTGCCCGCGAATCCGTGGGTGTATCCCCCGCCGCCTTCGGGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaa                                                                         <  2:1920221/143‑1 (MQ=255)
                                                                              ccgcgGATCCCGTGGCATGTCGTGCTCGTGGGGGTGCGGGCTTGCCCGCTAATCCGTTGGTGGATTCCCCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                        <  2:4148284/149‑1 (MQ=255)
                                                                               cgcgTATCCCGTGGCATGTCGAGATCGTGTGGGAGGGGGTTGGCCCGCTAATACGTTGTTGTATTCCCCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAg                                                                       <  2:458182/150‑1 (MQ=255)
                                                                               cgcgTATCCCGAGCCATGTCTTTATCCGGTGGGCGGGGGCTGGCCCGCGAATACGTGGGTGTTTACACCGCCGCTTTCGGGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAg                                                                       <  2:458235/150‑1 (MQ=255)
                                                                                  ggATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAAAACGATGGTGGGTTAAACCCCGCATTCGCGGGCCAGCCCCCGCCCACACTGTTCGGCGCAAACCCCCGGATCACATGAATCGAGGTTTTGACCAGAAAGttt                                                                    >  2:4060196/1‑150 (MQ=255)
                                                                                   gaccccGCGGCATGTCGATATCCGGGGGGCGCGGGGTTGCCCGCGAATACGAGGGTGTATTCACCGCCGATTTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                   <  2:4184556/147‑1 (MQ=255)
                                                                                     ccccGCGGCACGCCGCCATCGTGGGGGCGCGGGCTGCCCCGCTAATCGGTGGGTGTATCCCCCGCCGCTTGCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                                 <  1:1744151/149‑1 (MQ=255)
                                                                                        cGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCCAATACCATGGTGGATTAACCGCCCCATTCGCGGGCAACCCCGCTCCCCCCCTGTTCGGCGCAAGCCACGCGCTCACGTGAAAAGAGGTTTTTAAAAGAAATGTTTTatca                                                              >  2:3761338/1‑150 (MQ=255)
                                                                                         gAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACTATGGTGGATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat                                                             <  2:3002191/150‑1 (MQ=255)
                                                                                           ggCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg                                                           <  2:4424045/150‑1 (MQ=255)
                                                                                                     gttcgtgtgGTTGCGGGCTGGCCCGCTAATCGGTTGGTGTATTCACCGCCGAATCGGGGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg                                                 <  1:702553/145‑1 (MQ=255)
                                                                                                       tCCTGTGGGCGGGGGTTTGCCCGCCAATCCTCTGTTGGATCCACCGCCGCCTGCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                               <  1:1813422/150‑1 (MQ=255)
                                                                                                         cggggggAGCGGGCTTGCCCGCGAACCCGTTGGTGGATACCCCGCCGCTTGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                             <  2:887992/146‑1 (MQ=255)
                                                                                                            tgGGAGCGGGCTTGCCCGCGAATACGAGGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCCCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCgccgc                                          <  2:2925506/150‑5 (MQ=255)
                                                                                                                 gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCCCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGc                                     >  2:2891317/1‑150 (MQ=255)
                                                                                                                   gggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGccc                                   >  1:4481331/1‑150 (MQ=255)
                                                                                                                     gCTTGCCCGCGAATACTATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                                 <  2:3744092/150‑1 (MQ=255)
                                                                                                                       ttGCCCGCGAATACGATGGTGGATTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                               <  1:4080672/150‑1 (MQ=255)
                                                                                                                       ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                               >  1:1584723/1‑150 (MQ=255)
                                                                                                                       ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                               >  1:3860159/1‑150 (MQ=255)
                                                                                                                       ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                               >  1:315952/1‑150 (MQ=255)
                                                                                                                         gCCCGCGACTACGATGGTGGATTCACCGCCGCCGGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                             <  2:2066627/150‑1 (MQ=255)
                                                                                                                         gCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                             <  2:3691173/150‑1 (MQ=255)
                                                                                                                           ccGCTAATACGTTGGTGGTTACACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           <  1:2996772/150‑1 (MQ=255)
                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:2906591/1‑150 (MQ=37)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:2720185/1‑150 (MQ=255)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  1:511302/1‑150 (MQ=255)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:406721/1‑150 (MQ=255)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:464909/1‑150 (MQ=255)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  1:2470864/1‑150 (MQ=255)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:3469460/1‑150 (MQ=255)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  1:2238921/1‑150 (MQ=255)
                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:3886257/1‑150 (MQ=255)
                                                                                                                            cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                          >  2:3044680/1‑150 (MQ=255)
                                                                                                                              ggAATACTCTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                        <  2:2617374/149‑1 (MQ=255)
                                                                                                                              cgAATACGCTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGCTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                        <  2:2617459/150‑1 (MQ=255)
                                                                                                                               taaTCCGCTGGTGGATTCCCCGCCGCATGCGCGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGatca                       <  2:2018678/149‑4 (MQ=255)
                                                                                                                               gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg                       <  2:2595389/150‑1 (MQ=255)
                                                                                                                                aaTACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGt                      <  2:1843071/150‑1 (MQ=255)
                                                                                                                                    cgttggtggATTCCCCGCCGCTTTGGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc                  <  2:272727/147‑1 (MQ=255)
                                                                                                                                     gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                 >  2:4025098/1‑150 (MQ=255)
                                                                                                                                     gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                 >  2:817004/1‑150 (MQ=255)
                                                                                                                                     gATGGTGGACTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                 >  2:2507300/1‑150 (MQ=255)
                                                                                                                                      aTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGACAGTCGcccc                >  2:1912162/1‑149 (MQ=255)
                                                                                                                                       tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg               >  2:723563/1‑150 (MQ=255)
                                                                                                                                         gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa             <  2:4386876/150‑1 (MQ=255)
                                                                                                                                         gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa             >  1:1292298/1‑150 (MQ=255)
                                                                                                                                         gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa             >  1:226375/1‑150 (MQ=255)
                                                                                                                                         gtggATACCCCCCCCCTTGCGCGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa             <  1:3798898/150‑1 (MQ=255)
                                                                                                                                          tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc            >  1:338713/1‑150 (MQ=255)
                                                                                                                                          gggATTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc            <  1:3721166/149‑1 (MQ=255)
                                                                                                                                           ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc           >  1:1347043/1‑150 (MQ=255)
                                                                                                                                            gATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                   <  1:288727/93‑1 (MQ=255)
                                                                                                                                            gATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                   >  2:288727/1‑93 (MQ=255)
                                                                                                                                            aatccACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAccc          >  2:3958181/5‑150 (MQ=255)
                                                                                                                                              ttCACCGCCGCTTTCGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc        <  1:3230158/150‑1 (MQ=255)
                                                                                                                                              ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc        >  1:3318336/1‑150 (MQ=255)
                                                                                                                                              tcccccgcctCTTTCGGGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc        <  1:2745592/146‑1 (MQ=255)
                                                                                                                                               tCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcg       >  1:2995129/1‑150 (MQ=255)
                                                                                                                                                 aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTGTGACCTGAAAGTTGTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg     >  2:3870621/1‑150 (MQ=255)
                                                                                                                                                   cgccgcTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc   <  1:2885867/150‑1 (MQ=255)
                                                                                                                                                   cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc   >  2:109656/1‑150 (MQ=255)
                                                                                                                                                    gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg  >  1:2432906/1‑150 (MQ=255)
                                                                                                                                                    gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg  >  2:3842588/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCG  >  NC_002947/5486989‑5487286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: