Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,164 G→A 16.2% 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,490,1640GA16.2% 384.3 / 56.7 174intergenic (+73/‑74)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (53/86);  new base A (13/12);  total (75/99)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-02

Only 100 of 185 total aligned reads displayed.

ATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCG  >  NC_002947/5490015‑5490313
                                                                                                                                                       |                                                                                                                                                      
aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                                                        <  1:6533838/150‑1 (MQ=255)
    tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                                    <  1:2110179/150‑1 (MQ=255)
        gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTc                                                                                                                                                <  1:584186/150‑1 (MQ=255)
         tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGCAAACGAAGGTTGTTTAAACGCCCCACTcg                                                                                                                                               >  2:2291341/1‑150 (MQ=255)
          caccGCTGCTGCCGACGAGGCCGGCAACGCGATGCTCTTCACTGGCATGCGCCACTTCCGCCACTAATGACGCGGATCCCGAGTCAAGTCGAGATCCTGTCGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                              <  1:6027056/147‑1 (MQ=255)
          cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                              <  1:497805/150‑1 (MQ=37)
          cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                              <  2:2209090/150‑1 (MQ=255)
          cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                              <  1:3011337/150‑1 (MQ=255)
           aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                             <  2:5639204/150‑1 (MQ=255)
           aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCAAACCATTGTGGGTTCACCCCCCCAGTcgcg                                                                                                                                             >  2:1540545/1‑150 (MQ=255)
            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                            <  2:2573542/150‑1 (MQ=37)
            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                            <  1:2817720/150‑1 (MQ=255)
            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                            <  2:1171947/150‑1 (MQ=255)
            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                            <  1:3825426/150‑1 (MQ=255)
            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                            <  1:1138718/150‑1 (MQ=255)
            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                            <  1:6755544/150‑1 (MQ=255)
            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                            <  2:3165218/150‑1 (MQ=255)
             cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                            >  2:612090/1‑149 (MQ=37)
             cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                            <  1:612090/149‑1 (MQ=37)
                tgctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCTAGTCGAGATCCTGTGGGAGCGGGCGTGGCCGCGCATGCCTGGGGGGCGTGACCGCCGGAGTCGCGGGcgc                                                                                                                                        >  2:4263945/1‑148 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                                     <  1:6454458/150‑1 (MQ=255)
                      cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGc                                                                                                                                 <  2:2433062/150‑1 (MQ=21)
                             gcgggCAGCGGGTTGGGCTTCAGTGGCATGCGCCACTTCCGCCAATATTTACGCGTTGCCCGTGGCTTGTGTAGTGCCGTTGGGAGGGGGCTTGCCCGCTATTTCGATGGTGTTTTCCCCGCCGCATTCGGGGCCATGCCCGCTCCcact                                                                                                                          <  2:5034711/147‑2 (MQ=255)
                             gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCCAAACATGGTGGATTCACCGCCCCCTTAGCGGGCAAGCCCCCGCCcaca                                                                                                                          >  2:6627116/1‑150 (MQ=255)
                             gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGAAAAAGATGGTGGTTTGAACGCCGCAGTGGCGGGGCAGCCCGCGCCcaaa                                                                                                                          >  2:6557934/1‑148 (MQ=255)
                               cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTGCCCCGCAAAAGCTGGTGGGTTAAACCCCGCCATAGCGGGGACGCCCGCTCCCACAgg                                                                                                                        >  2:5798557/1‑149 (MQ=255)
                                  cATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCCAACGATGGGGGGTTAACCCCCCGCATCGCGGGCACGCCCGCGCCCACCGGGtt                                                                                                                     >  2:4871634/1‑150 (MQ=255)
                                  cATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCGGTGGGAGCGGGGGTGCCCGCGAATACGTTGGGGGTTTCACCGCCGCCGTCGCGGGCAAGCCCGCGCCCACAGTGtt                                                                                                                     >  2:3362010/1‑150 (MQ=255)
                                    tCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGc                                                                                                                                          >  1:4625890/1‑128 (MQ=255)
                                    tCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGc                                                                                                                                          <  2:4625890/128‑1 (MQ=255)
                                     cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                              >  1:6057519/1‑123 (MQ=21)
                                     cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCgg                                                                                                                  <  2:6269125/150‑1 (MQ=255)
                                                tCAGTGGCTTGCGCCACTCCCGCAACTAATTACGGGGAACCCGAGGCATGTCGAGATCCTGTGGGAGGGGGCTTGCCCGCTAATAGGTTGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                                       <  2:1995069/150‑1 (MQ=255)
                                                 cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGCATACGCTGGTGGGTTGACCCCCCCCGTAGCGGGCCACCCCCCGCCCACCCGGTTCGGCGCAAGcccccg                                                                                                      >  1:5623480/1‑146 (MQ=255)
                                                 cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGGGGGAGCGGGCTGGCCCGCGAAACCGATGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                                      <  2:6713245/150‑1 (MQ=255)
                                                  ccTGGCAGGCGCCACTCCCGCCACTATTTACGCGGATCCCGCGGAAAGTCGAGATCCTGTGGGTGGGGGCTGGCCCGCTAATCCTTTGGTGGATTCACCGCCGCAGGCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACggg                                                                                                     <  2:2962784/149‑1 (MQ=255)
                                                      gCATGCGCCACTACCGCCACTAATCCCGCGGATCCCGAGGCACGTCGAGATCCTGTGGGAGCGGGGTTGTCCGCGAATCCGATGGTGGATTCCCCGCGGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCa                                                                                                 <  2:1062410/150‑1 (MQ=255)
                                                         tGCGCCACTCCCGCCACTAAGTACGCGGATCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                                              <  2:2698030/150‑1 (MQ=255)
                                                            gcCACTTCCGCCACTAATTACGGGGATCCCGAGGCATGCCGAGATCCTGTGGGAGGGGGCGTGCCCGCTAATCCTTTGGTGTATACCCCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGaa                                                                                           <  2:5810548/150‑1 (MQ=255)
                                                              cACTCCCGCCACTAATTCGGCGGATCCCGAGGCAAGTCGAGAGCGTGTGGGCGCGGGCTGGCCCGCTAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATc                                                                                         <  2:4734542/150‑1 (MQ=38)
                                                                cTTCCCCCAATTATTCCGCGGAGCCCGAGGCACGTCTCGATCCTGTGGGAGGGGGCTGGCCCGCGAAAAGGTTGGTGTATTCCCCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                       <  2:5949536/150‑1 (MQ=255)
                                                                 ttCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                      <  2:4461148/150‑1 (MQ=255)
                                                                      ccACTACTTCCGCGGAGCCCGCGGCACGTCGCGTTCCTGTGGGGGGGGGGTGGCCCGGGAATCCGTTGGTGGATTCCCCGCCGCTTTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                 <  2:715474/150‑1 (MQ=255)
                                                                      ccAATAATTACGCGGATCCCGAGGCAAGTCGATATCCCGTGGGAGCGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                 <  2:420074/150‑1 (MQ=38)
                                                                        aCTAATTAGGCGGATCCCGAGGCATGTCGTGTTCCTGTGGGAGGGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                                               <  2:6077532/150‑1 (MQ=255)
                                                                                 gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCCCATAACATTGTGGATTAACCGCCGCAATCGCCGGCCAGCCCCCGCCCACACTGGTAGGCGCCAACCCCCGGATCCAGTTAATCGAGGGTTTTACAAGaaact                                                                      >  2:5525835/1‑148 (MQ=255)
                                                                                 gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCCAAAACCCTGGGGGGATCACCCCCCCAGTCGCCGGCAAGCCCGCACCCACAGGGTTCGGCGCCACCCACCGGGTGACGTGGATCGAGGGTGTGGAATGAAAGt                                                                      >  2:6637500/1‑150 (MQ=255)
                                                                                  cgGATCCCTCGCCAAGTCGATTTCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGTTTTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtt                                                                     <  2:2398478/150‑1 (MQ=37)
                                                                                    gATCCCGAGGCACGTCGATATCGTGTGGGAGCGGGCTTGCCCGCTAATCCGAGGGTGGATTCACCGCCGCATTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGAAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                   <  2:4087096/150‑1 (MQ=38)
                                                                                      tCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCTAATACGAGGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                                 <  2:2672038/150‑1 (MQ=255)
                                                                                          gAGGCATGTCGATATCCTGTGGGTGGGGGCTTGCCCGCTAATCCTTTGGTGTATTCACCGCCGCTTTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat                                                             <  2:1813115/150‑1 (MQ=255)
                                                                                           cggaacgTCGAGATCTTGTGGGTGGGGGCTGGCCCGCTAATACTTTGGTGTTTTCCCCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                                            <  1:1619379/144‑1 (MQ=255)
                                                                                            ggCACGTCGAGATCCGGGGGGCGGGGGCTTGCCCGCTAATCCGTTGGTGGATACACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg                                                           <  2:6330682/150‑1 (MQ=255)
                                                                                                  tCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg                                                     <  2:3986755/150‑1 (MQ=255)
                                                                                                        tCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                               <  2:3942398/150‑1 (MQ=255)
                                                                                                          gtgggggTGCGGGCTGGCCCGCAAATAGTTGGTGGTATCCCCCGCCCAAGTCGGGGCAACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                             <  1:2643032/146‑1 (MQ=255)
                                                                                                            gtgGGCGCGGGCTGGCCCGCTAATACGTTGGTGGATTCCCCGCCGCATGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                           <  2:5355940/150‑1 (MQ=255)
                                                                                                            gtgGGAGGGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                           <  2:4325862/150‑1 (MQ=255)
                                                                                                            gtgGGAGCGGGCTTGCCCGCTAATACGATGGTGGATCCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                           <  2:3270326/150‑1 (MQ=255)
                                                                                                            gtgGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                           <  2:500013/150‑1 (MQ=255)
                                                                                                            gtgGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                           <  2:6190197/150‑1 (MQ=255)
                                                                                                              attcgCGGGCAAGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGc                                         >  1:1753942/5‑150 (MQ=255)
                                                                                                               ggAGCGGGCTTGCCCGCTAATACGTTGGTGTATCCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCa                                        <  2:944350/150‑1 (MQ=39)
                                                                                                                       cTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                                >  2:508340/1‑150 (MQ=255)
                                                                                                                        ttGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                               >  1:3622929/1‑150 (MQ=25)
                                                                                                                        ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                               >  2:4394767/1‑150 (MQ=255)
                                                                                                                          gCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                             >  1:1805348/1‑150 (MQ=255)
                                                                                                                          gCCCGCCAATACGTTGTTGTATACACCGCCGACTGCGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                             <  1:543441/150‑1 (MQ=255)
                                                                                                                            ccacgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCATGAGCACGCcctggcctgg                           >  1:2383298/4‑150 (MQ=25)
                                                                                                                            ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:383699/1‑150 (MQ=255)
                                                                                                                            ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:1862786/1‑150 (MQ=255)
                                                                                                                            ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  1:5439640/1‑150 (MQ=255)
                                                                                                                            ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                           >  2:6230590/1‑150 (MQ=255)
                                                                                                                            ccGCGAATACGATGGAGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGTcctggcctgg                           >  1:4071692/1‑150 (MQ=255)
                                                                                                                              gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                         >  2:1229975/1‑150 (MQ=37)
                                                                                                                              gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                         >  2:6784992/1‑150 (MQ=37)
                                                                                                                              gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                         >  2:825610/1‑150 (MQ=255)
                                                                                                                              gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                         >  2:4431699/1‑150 (MQ=255)
                                                                                                                              gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                         >  1:4130094/1‑150 (MQ=255)
                                                                                                                              gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                         >  2:3483852/1‑150 (MQ=255)
                                                                                                                               cgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                        <  1:2137497/150‑1 (MQ=37)
                                                                                                                                 aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGt                      <  2:6207158/150‑1 (MQ=255)
                                                                                                                                  atgcggcggtgaATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc                     <  1:3806842/143‑1 (MQ=37)
                                                                                                                                     cGAGGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc                  <  2:6042725/150‑1 (MQ=255)
                                                                                                                                      gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCTCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                 >  2:67397/1‑150 (MQ=255)
                                                                                                                                      gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCTCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                 >  2:1364766/1‑150 (MQ=255)
                                                                                                                                      gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg                                                     >  2:5690887/1‑114 (MQ=255)
                                                                                                                                      gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg                                                     <  1:5690887/114‑1 (MQ=255)
                                                                                                                                      gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                 >  1:954111/1‑150 (MQ=255)
                                                                                                                                       ttggtggATTCACCGCCGCCTGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa                <  1:6572881/149‑1 (MQ=255)
                                                                                                                                        tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGACAGTCGCCCAg               >  2:5334083/1‑150 (MQ=255)
                                                                                                                                            ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc           >  1:5738080/1‑150 (MQ=255)
                                                                                                                                            ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCAGGCCTGGAAAGTCGCCCAGGAcc           >  1:1225150/1‑150 (MQ=255)
                                                                                                                                                ttCCCCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc        <  2:3760384/150‑1 (MQ=255)
                                                                                                                                                ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc        >  1:4877936/1‑150 (MQ=255)
                                                                                                                                                tacACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc        <  1:1115036/148‑1 (MQ=255)
                                                                                                                                                   aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg     >  2:3288572/1‑150 (MQ=255)
                                                                                                                                                     ccgccgCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGt    >  2:1618556/1‑150 (MQ=255)
                                                                                                                                                      cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc   <  2:1702338/150‑1 (MQ=255)
                                                                                                                                                       gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTGTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCg  >  2:4709947/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
ATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCG  >  NC_002947/5490015‑5490313

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: