Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,519 G→A 7.5% H132H (CAC→CAT) ‡ amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,5190GA 7.5% 1100.7 / 52.2 361H132H (CAC→CAT) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (152/181);  new base A (11/16);  total (163/198)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 372 total aligned reads displayed.

GTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGA  >  NC_002947/6005373‑6005667
                                                                                                                                                  |                                                                                                                                                    
gTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATgtgg                                                                                                                                                   >  2:5490943/1‑150 (MQ=255)
gTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATgtgg                                                                                                                                                   >  1:8836216/1‑150 (MQ=255)
        aGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCgg                                                                                                                                           <  2:8908160/150‑1 (MQ=255)
        aGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCgg                                                                                                                                           <  1:7546974/150‑1 (MQ=255)
          tcttctCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGgc                                                                                                                                         >  2:4056642/1‑150 (MQ=255)
            ttctCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGgcgc                                                                                                                                       >  2:6779274/1‑150 (MQ=255)
             tctcGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCt                                                                                                                                      >  1:4790434/1‑150 (MQ=255)
              ctcGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTc                                                                                                                                     <  2:6118475/150‑1 (MQ=255)
                 ggACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGa                                                                                                                                  >  2:8074514/1‑150 (MQ=255)
                 ggACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGa                                                                                                                                  >  1:5591712/1‑150 (MQ=255)
                 ggACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGa                                                                                                                                  <  1:4118885/150‑1 (MQ=255)
                 ggACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGa                                                                                                                                  >  2:6587039/1‑150 (MQ=255)
                 ggACGGGTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGGGCCAGGTTTCACGCAGGTGGTGGCCCGGGCGCTCGGa                                                                                                                                  <  2:5449326/150‑1 (MQ=255)
                     ggTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCg                                                                                                                              <  1:5947220/150‑1 (MQ=255)
                      gTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGa                                                                                                                             >  2:199741/1‑150 (MQ=255)
                        ttCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATg                                                                                                                           <  2:6904084/150‑1 (MQ=255)
                        ttCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATg                                                                                                                           <  1:3707720/150‑1 (MQ=255)
                         tCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGGTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATgg                                                                                                                          >  2:2966051/1‑150 (MQ=255)
                          ccATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGt                                                                                                                         <  1:931182/150‑1 (MQ=255)
                          ccATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGt                                                                                                                         <  2:6219140/150‑1 (MQ=255)
                           cATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTc                                                                                                                        >  1:112864/1‑150 (MQ=255)
                            aTACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCa                                                                                                                       >  1:447484/1‑150 (MQ=255)
                            aTACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCa                                                                                                                       >  2:3103831/1‑150 (MQ=255)
                             tACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTcga                                                                                                                      <  2:869902/150‑3 (MQ=255)
                             tACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAg                                                                                                                      >  2:628923/1‑150 (MQ=255)
                             tACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAg                                                                                                                      <  2:7473937/150‑1 (MQ=255)
                                caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGcc                                                                                                                   >  2:4973838/1‑150 (MQ=255)
                                 accGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCg                                                                                                                  >  1:9066150/1‑150 (MQ=255)
                                    gAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTAc                                                                                                               <  1:628923/150‑1 (MQ=255)
                                      gTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAg                                                                                                             <  1:5373998/150‑1 (MQ=255)
                                             gCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGgcgc                                                                                                      <  2:1459383/150‑1 (MQ=255)
                                                accaGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGTGATGGTCAGGCCGTACAGCTGGCGCGAc                                                                                                   <  2:5742993/150‑1 (MQ=255)
                                                accaGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGAc                                                                                                   <  1:6563874/150‑1 (MQ=255)
                                                 ccagggccagggCGGTGTTCCAGGCACACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACa                                                                                                  <  1:742236/150‑1 (MQ=255)
                                                  cagggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAg                                                                                                 >  1:6332053/1‑150 (MQ=255)
                                                  cagggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAg                                                                                                 >  2:1755086/1‑150 (MQ=255)
                                                  cagggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAg                                                                                                 <  2:1151541/150‑1 (MQ=255)
                                                    gggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGc                                                                                               <  1:7320375/150‑1 (MQ=255)
                                                    gggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGc                                                                                               <  1:8471537/150‑1 (MQ=255)
                                                                  ttCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCg                                                                                 >  1:5784461/1‑150 (MQ=255)
                                                                         cccACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGt                                                                          <  1:804603/150‑1 (MQ=255)
                                                                          ccACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTc                                                                         <  1:1767457/150‑1 (MQ=255)
                                                                           cacaCGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCg                                                                        >  1:8065485/1‑150 (MQ=255)
                                                                            acacGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGa                                                                       >  1:6615311/1‑150 (MQ=255)
                                                                              acGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATc                                                                     >  1:3762884/1‑150 (MQ=255)
                                                                              acGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATc                                                                     >  1:2812513/1‑150 (MQ=255)
                                                                                   aaCCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTc                                                                <  2:1839720/150‑1 (MQ=255)
                                                                                   aaCCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTc                                                                <  1:4027466/150‑1 (MQ=255)
                                                                                   aaCAGGGAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTc                                                                <  1:2911208/150‑1 (MQ=255)
                                                                                     cagggAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCt                                                              <  2:5245896/148‑1 (MQ=255)
                                                                                      cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                                             <  1:1945841/150‑1 (MQ=255)
                                                                                        ggAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTgta                                                           >  1:7796533/1‑148 (MQ=255)
                                                                                        ggAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTgta                                                           >  2:3725814/1‑148 (MQ=255)
                                                                                              tGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCAc                                                     >  1:2470956/1‑150 (MQ=255)
                                                                                                  ggccgagATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAg                                                 <  1:6082535/144‑1 (MQ=255)
                                                                                                  ggccgagATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAg                                                 <  2:2653672/144‑1 (MQ=255)
                                                                                                  ggCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAg                                                 <  2:2812513/150‑1 (MQ=255)
                                                                                                       tGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCtt                                            >  2:6920478/1‑150 (MQ=255)
                                                                                                        gATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTc                                           <  1:5354851/150‑1 (MQ=255)
                                                                                                          tCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCt                                         >  1:2617689/1‑150 (MQ=255)
                                                                                                           aaCGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTccga                                        >  1:6660977/2‑148 (MQ=255)
                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGACGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACGTCACCCAGGCCTTCCTgg                                       <  1:1313821/150‑1 (MQ=255)
                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                       <  2:4480921/150‑1 (MQ=255)
                                                                                                            aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                       <  1:4930894/150‑1 (MQ=255)
                                                                                                             cGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTccgaag                                      <  2:348707/150‑5 (MQ=255)
                                                                                                             cGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTccgaag                                      >  1:3534613/1‑146 (MQ=255)
                                                                                                              gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGt                                     >  1:3414575/1‑150 (MQ=255)
                                                                                                                 gACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTGaa                                  <  2:3823442/150‑2 (MQ=255)
                                                                                                                 gACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTGaa                                  <  2:3459052/150‑2 (MQ=255)
                                                                                                                 gACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTGaa                                  <  2:5880238/150‑2 (MQ=255)
                                                                                                                  aCCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATc                                 >  2:8541205/1‑150 (MQ=255)
                                                                                                                   ccACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCt                                >  2:5333236/1‑150 (MQ=255)
                                                                                                                   ccACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCt                                <  2:219932/150‑1 (MQ=255)
                                                                                                                    cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                               >  2:4196451/1‑150 (MQ=255)
                                                                                                                    cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                               >  1:5625327/1‑150 (MQ=255)
                                                                                                                       acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                            >  2:4596497/1‑150 (MQ=255)
                                                                                                                       acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                            >  2:6368799/1‑150 (MQ=255)
                                                                                                                       acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                            >  1:6215591/1‑150 (MQ=255)
                                                                                                                       acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                            >  1:1593209/1‑150 (MQ=255)
                                                                                                                         ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCgg                          >  2:5972691/1‑150 (MQ=255)
                                                                                                                            gCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCtt                       >  1:5111084/1‑150 (MQ=255)
                                                                                                                             cGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtc                      <  1:2688651/150‑1 (MQ=255)
                                                                                                                              gggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcg                     >  1:3092489/1‑150 (MQ=255)
                                                                                                                               ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                    >  2:4935568/1‑150 (MQ=255)
                                                                                                                                 tGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc                  <  1:8910914/150‑1 (MQ=255)
                                                                                                                                 tGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc                  <  1:6389823/150‑1 (MQ=255)
                                                                                                                                  gCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcg                 <  2:9100440/150‑1 (MQ=255)
                                                                                                                                   ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                <  1:1003614/150‑1 (MQ=255)
                                                                                                                                   ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                <  1:5490943/150‑1 (MQ=255)
                                                                                                                                   ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                <  2:8240025/150‑1 (MQ=255)
                                                                                                                                   ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                <  1:5284411/150‑1 (MQ=255)
                                                                                                                                   ccAGGATTCACGCATGGGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                <  2:7367853/150‑1 (MQ=255)
                                                                                                                                    cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag               >  1:7052696/1‑150 (MQ=255)
                                                                                                                                     aGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgaga              >  2:2351589/1‑150 (MQ=255)
                                                                                                                                      ggTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGAc             <  1:199761/150‑1 (MQ=255)
                                                                                                                                          tCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGc         <  1:6906315/150‑1 (MQ=255)
                                                                                                                                          tCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGc         <  1:2287577/150‑1 (MQ=255)
                                                                                                                                           cACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCt        <  2:5373447/150‑1 (MQ=255)
                                                                                                                                            aCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGc                             <  2:7853110/128‑1 (MQ=255)
                                                                                                                                                 tgtgGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGa  <  2:4729787/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGA  >  NC_002947/6005373‑6005667

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