Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,487,141 C→T 8.9% intergenic (+77/‑70) 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,1410CT 8.9% 309.4 / 17.2 116intergenic (+77/‑70)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (46/57);  new base T (2/8);  total (49/67)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

Only 100 of 131 total aligned reads displayed.

TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAA  >  NC_002947/5486992‑5487290
                                                                                                                                                      |                                                                                                                                                     
tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGGATACGGGGGTGGGTTGAAcgccgc                                                                                                                                                       >  2:2430910/1‑150 (MQ=255)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                                 <  2:1297014/150‑1 (MQ=37)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCATAAGATTGTGGTTTCACcgccgccgtcccg                                                                                                                                                >  1:1843923/1‑143 (MQ=255)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                                <  2:3254426/150‑1 (MQ=255)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                                <  1:718999/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                               <  2:3978050/150‑1 (MQ=37)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                               <  2:302141/150‑1 (MQ=37)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                               <  2:1223951/150‑1 (MQ=37)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:1815277/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:1883869/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:2149980/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                               <  1:2503146/150‑1 (MQ=25)
         cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCggg                                                                                                                                              <  2:903717/150‑1 (MQ=37)
         cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                                              <  2:844110/150‑1 (MQ=255)
         cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                                              <  2:760849/150‑1 (MQ=255)
          gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCggg                                                                                                                                              <  2:1273143/150‑1 (MQ=14)
            tgctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCaa                                                                                                                                           <  1:288824/150‑1 (MQ=255)
             gctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCCCGCAAAAGAGGGTGGTTTAAACCCCCCCGTCGCGGGCacg                                                                                                                                          >  1:3002191/1‑148 (MQ=255)
              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                                         <  1:1246499/150‑1 (MQ=255)
              ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                                         <  2:3133775/150‑1 (MQ=255)
                 ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATTCGATGGTGGAATCACCCCCGCAATCGCGGGCAAGCCCg                                                                                                                                      >  2:3747229/1‑150 (MQ=255)
                      gaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCATGCCCGCGCATACGCTGGTGGCTTCACCGCCGCATTCGCGGGCACGCCCGCTccc                                                                                                                                 >  2:4477444/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)
                                              aCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCCCGCATACCAGGGGGGGGTCACCGCCCCAATCGCCGGCCACCCCGCGCCCCCCCGGTTAGGCGCCAGCCCCgcg                                                                                                         >  2:4325111/1‑148 (MQ=255)
                                                     ggcgcCCCTCCCGCCACTACTTCCGGGGATCCGGCGGCACGTCGCGCTCCGGTGGGCGCGGGCTTGCCGGCTAATCCGATGGTGTATTCCCGGCCGCCTGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                                                  <  1:4119238/149‑1 (MQ=255)
                                                            ctgccGCCACTAATCCCGGGGATCCGGCGGCACGGCGATATCCTGTGGGTGGGGGCTTGCCCGCCAATCCGTGGGTGGATACCCCGCCGCTTGGGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCGCTTGAATCGa                                                                                           <  2:2264030/147‑1 (MQ=255)
                                                            cTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCCCATAAGATTGTGGGTTAACCGCCCCCCTACCCGGCCAGCCCCCGCCCCCCGGGTTGGGCGCAAGCCACGGGGGCACATGaaaaga                                                                                           >  2:2749636/1‑146 (MQ=255)
                                                            cTTCCGCCACTAATGACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAACCGATGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                           <  2:4106306/150‑1 (MQ=255)
                                                             ttCCGCCACTACTTACGCGGATCCCGCGGAAAGTCGAGATCCTGTGGGCGCGGGCTTGCCCGCGAATACGAGGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                          <  2:1070681/150‑1 (MQ=38)
                                                                cGCCACTATTTAGGCGTATCCCGTGGAATGTCGATATCGTGTGGGTGGGGGCTTGCCCGCTAATACGTTGGTGTTTTCCCCGCCGCTTTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtt                                                                                       <  2:637132/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)
                                                                                 aTCCGGAGGAATGTCTAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                                      <  2:2758908/150‑1 (MQ=38)
                                                                                  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)
                                                                                                       cTGGGGGCGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                                 <  2:3886259/150‑1 (MQ=255)
                                                                                                               gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCCCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGc                                         >  2:2891317/1‑150 (MQ=255)
                                                                                                                   gCTTGCCCGCGAATACTATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                                     <  2:3744092/150‑1 (MQ=255)
                                                                                                                     ttGCCCGCGAATACGATGGTGGATTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                                   <  1:4080672/150‑1 (MQ=255)
                                                                                                                     ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                                   >  1:315952/1‑150 (MQ=255)
                                                                                                                     ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                                   >  1:3860159/1‑150 (MQ=255)
                                                                                                                     tggCCCGCGAATACGATGGTGGATTCACCCCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                                   <  1:3372811/148‑1 (MQ=255)
                                                                                                                      tGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcc                                  <  2:3346734/150‑1 (MQ=255)
                                                                                                                       gCCCGCGACTACGATGGTGGATTCACCGCCGCCGGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                 <  2:2066627/150‑1 (MQ=255)
                                                                                                                       gCCCGCGAATCCGCTGGTGGATTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                 <  2:2066436/150‑1 (MQ=255)
                                                                                                                       gCCCGCGAATACGCTGGTGGATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                 <  2:309534/150‑1 (MQ=255)
                                                                                                                       gCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                 <  2:3691173/150‑1 (MQ=255)
                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:2906591/1‑150 (MQ=37)
                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:3148962/1‑150 (MQ=37)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:2470864/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:406721/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:3469460/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:2720185/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:2238921/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  1:511302/1‑150 (MQ=255)
                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                               >  2:3886257/1‑150 (MQ=255)
                                                                                                                          cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                              >  1:1452804/1‑150 (MQ=37)
                                                                                                                          cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                              >  2:3044680/1‑150 (MQ=255)
                                                                                                                            ggAATACTCTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                            <  2:2617374/149‑1 (MQ=255)
                                                                                                                            cgAATACGCTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGCTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                            <  2:2617459/150‑1 (MQ=255)
                                                                                                                             gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg                           <  2:2595389/150‑1 (MQ=255)
                                                                                                                              aaTACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGt                          <  2:1843071/150‑1 (MQ=255)
                                                                                                                                tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCg                        >  1:3478316/1‑150 (MQ=255)
                                                                                                                                 aCGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc                       <  2:1843923/150‑1 (MQ=255)
                                                                                                                                  cgttggtggATTCCCCGCCGCTTTGGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc                      <  2:272727/147‑1 (MQ=255)
                                                                                                                                   gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                     >  2:4025098/1‑150 (MQ=255)
                                                                                                                                   gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                     >  2:4025065/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:226375/1‑150 (MQ=255)
                                                                                                                                       gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa                 >  1:1292298/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)
                                                                                                                                            ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc            >  1:4490181/1‑150 (MQ=255)
                                                                                                                                            ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc            >  1:3318336/1‑150 (MQ=255)
                                                                                                                                            tcccccgcctCTTTCGGGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc            <  1:2745592/146‑1 (MQ=255)
                                                                                                                                             tCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcg           >  1:2995129/1‑150 (MQ=255)
                                                                                                                                               aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACGTGACTCGAGGTTTTGACATGAAAGGTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg         >  2:2911739/1‑150 (MQ=255)
                                                                                                                                                  gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg      >  2:3842588/1‑150 (MQ=255)
                                                                                                                                                  gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg      >  2:3842342/1‑150 (MQ=255)
                                                                                                                                                  gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg      >  1:2432906/1‑150 (MQ=255)
                                                                                                                                                   ccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCga     >  1:1958748/1‑150 (MQ=255)
                                                                                                                                                    cgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCgag    >  2:3772980/1‑150 (MQ=255)
                                                                                                                                                     gcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCgaga   >  1:4147172/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAA  >  NC_002947/5486992‑5487290

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: