Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,061,264 G→A 7.4% A35V (GCC→GTC)  hupA ← DNA‑binding protein HU‑alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,061,2640GA 7.4% 237.2 / 7.1 81A35V (GCC→GTC) hupADNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base A (4/2);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

TTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCT  >  NC_002947/6061120‑6061409
                                                                                                                                                 |                                                                                                                                                 
ttCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTAGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGCCCTGGACAAGGACAcc                                                                                                                                               <  2:335402/149‑1 (MQ=255)
 tCTTGCCAGGCTTGACCGCCACGGGAGTGCTGGCCTGGATCTTGACCGGTGCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCg                                                                                                                                              <  4:136662/149‑1 (MQ=255)
        aGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGc                                                                                                                                       <  8:144004/149‑1 (MQ=255)
            ttGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCATTGTCAAGGGCACCGGTGATGCTGTc                                                                                                                                   >  8:150031/1‑149 (MQ=255)
                   ccACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATc                                                                                                                            <  3:155953/149‑1 (MQ=255)
                    cACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATcg                                                                                                                           >  6:107526/1‑149 (MQ=255)
                    cACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATcg                                                                                                                           >  2:96448/1‑149 (MQ=255)
                     aCGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATcgc                                                                                                                          >  1:353767/1‑149 (MQ=255)
                             gCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCaaaa                                                                                                                  >  6:161738/1‑149 (MQ=255)
                                 gCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGGCCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAAc                                                                                                                 <  1:265269/147‑1 (MQ=255)
                                 gCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCt                                                                                                               >  1:332934/1‑148 (MQ=255)
                                   cTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGAt                                                                                                            >  6:369630/1‑149 (MQ=255)
                                    ttGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGAtt                                                                                                           >  8:231116/1‑149 (MQ=255)
                                      gATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGAtt                                                                                                           <  2:353767/147‑1 (MQ=255)
                                         cTTGACAGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCt                                                                                                      >  5:267282/1‑149 (MQ=255)
                                           tGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCttt                                                                                                    >  2:12855/1‑149 (MQ=255)
                                            gACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCTCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACAGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGcc                                                                                                       <  5:233982/145‑1 (MQ=255)
                                            gACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCtttt                                                                                                   >  3:265346/1‑149 (MQ=255)
                                            gACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTAAGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATAGCGTTCAAAACCTGATTGGCCtttt                                                                                                   >  5:220403/1‑149 (MQ=255)
                                              ccGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTcc                                                                                                 <  6:267282/149‑1 (MQ=255)
                                               cGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCt                                                                                                <  1:12855/149‑1 (MQ=255)
                                               cGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCt                                                                                                <  2:142786/149‑1 (MQ=255)
                                                 gTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTg                                                                                              <  5:143764/149‑1 (MQ=255)
                                                 gTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTg                                                                                              <  5:107526/149‑1 (MQ=255)
                                                    ccccAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGt                                                                                            <  6:254563/148‑1 (MQ=255)
                                                     cccAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGt                                                                                            >  1:86394/1‑147 (MQ=255)
                                                     cccAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGt                                                                                            >  5:173479/1‑147 (MQ=255)
                                                      ccAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGt                                                                                            >  1:27478/1‑146 (MQ=255)
                                                             gCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGcc                                                                                  >  3:102685/1‑149 (MQ=255)
                                                             gCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGc                                                                                   <  3:142973/148‑1 (MQ=255)
                                                               gggTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCt                                                                                 >  1:91458/1‑148 (MQ=255)
                                                                    cTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCg                                                                            <  1:96448/148‑1 (MQ=255)
                                                                     ttGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGc                                                                          >  6:338283/1‑149 (MQ=255)
                                                                       gCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCgg                                                                           <  2:346149/146‑1 (MQ=255)
                                                                                 gcACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACggcggc                                                              >  6:206640/1‑149 (MQ=255)
                                                                                     cgcgATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTc                                                                                           >  1:51914/1‑116 (MQ=255)
                                                                                     cgcgATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTc                                                                                           <  2:51914/116‑1 (MQ=255)
                                                                                          tGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTgg                                                     <  5:161738/149‑1 (MQ=255)
                                                                                             cGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGttt                                                  >  2:205883/1‑149 (MQ=255)
                                                                                              gTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTg                                                 >  1:144147/1‑149 (MQ=255)
                                                                                                 tttCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTgcg                                               >  8:242533/1‑148 (MQ=255)
                                                                                                   tCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCa                                             >  7:251625/1‑148 (MQ=255)
                                                                                                   tCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCa                                             >  2:253966/1‑148 (MQ=255)
                                                                                                     gAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAg                                          >  4:194136/1‑149 (MQ=255)
                                                                                                          tGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAg                                     >  6:86217/1‑149 (MQ=255)
                                                                                                             cGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACggcggc                                                              <  2:192/121‑1 (MQ=255)
                                                                                                             cGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACggcggc                                                              >  1:192/1‑121 (MQ=255)
                                                                                                             cGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCct                                  >  1:208110/1‑149 (MQ=255)
                                                                                                              gAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCTGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCctc                                 >  7:165973/1‑149 (MQ=255)
                                                                                                                aaCCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCtt                               >  8:255789/1‑149 (MQ=255)
                                                                                                                   cGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGa                            <  7:231116/149‑1 (MQ=255)
                                                                                                                    gACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGAc                           >  4:283690/1‑149 (MQ=255)
                                                                                                                    gACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGAc                           <  5:253372/149‑1 (MQ=255)
                                                                                                                    gACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGAc                           <  2:27478/149‑1 (MQ=255)
                                                                                                                     aCCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACg                          >  4:294535/1‑149 (MQ=255)
                                                                                                                     aCCAGGGTCACAGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACg                          >  1:319093/1‑149 (MQ=255)
                                                                                                                     aCCAGGGTCACAGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACg                          >  1:324169/1‑149 (MQ=255)
                                                                                                                         gggTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAAt                      >  1:351593/1‑149 (MQ=255)
                                                                                                                          ggTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAAtt                     >  6:317200/1‑149 (MQ=255)
                                                                                                                          ggTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAAtt                     <  3:294535/149‑1 (MQ=255)
                                                                                                                            caccGGTGCCCTGGACAAGGACACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCt                   >  5:196162/4‑149 (MQ=255)
                                                                                                                             cACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGt                                                    <  6:156724/115‑1 (MQ=255)
                                                                                                                             accGGTGCCCTGGACAAGGACACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCtt                  >  8:152735/3‑149 (MQ=255)
                                                                                                                              ccGGTGCCCTGGACAAGGACACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCa                                                                                          >  5:353719/2‑76 (MQ=37)
                                                                                                                              ccGGTGCCCTGGACAAGGACACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCa                                                                                          <  6:353719/75‑1 (MQ=37)
                                                                                                                              aCGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCtt                  <  8:165973/148‑1 (MQ=255)
                                                                                                                                ggTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCttgtt               <  5:206640/149‑1 (MQ=255)
                                                                                                                                ggTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCttgtt               <  4:265346/149‑1 (MQ=255)
                                                                                                                                ggTGCCCTGGACAAGGACACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCTTAGTGAAGCGTCTTTGAGGGAATTCttgtt               >  2:267855/1‑149 (MQ=255)
                                                                                                                                     tCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCctct                                >  3:147733/1‑128 (MQ=255)
                                                                                                                                     tCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCctct                                <  4:147733/128‑1 (MQ=255)
                                                                                                                                       cTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTAt          <  2:324169/148‑1 (MQ=255)
                                                                                                                                       cTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTAt          <  2:319093/148‑1 (MQ=255)
                                                                                                                                       cTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTAt          <  2:91458/148‑1 (MQ=255)
                                                                                                                                       cTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATg         >  4:195391/1‑149 (MQ=255)
                                                                                                                                          gTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGcc       <  7:242533/148‑1 (MQ=255)
                                                                                                                                           tCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGcctt     >  4:22210/1‑147 (MQ=255)
                                                                                                                                              aGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCttgtt               >  1:218481/1‑136 (MQ=255)
                                                                                                                                              aGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCttgtt               <  2:218481/136‑1 (MQ=255)
                                                                                                                                              aGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGtgct  <  4:102685/149‑1 (MQ=255)
                                                                                                                                              aGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGtgct  <  6:173479/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGTGCT  >  NC_002947/6061120‑6061409

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