Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,061,260 C→A 8.6% L36L (CTG→CTT 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,2600CA 8.6% 227.2 / 8.5 78L36L (CTG→CTThupADNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base A (5/2);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGT  >  NC_002947/6061115‑6061406
                                                                                                                                                  |                                                                                                                                                  
gCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAgg                                                                                                                                                 <  2:197806/149‑1 (MQ=255)
     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                                                                                                                        >  2:96448/1‑149 (MQ=255)
                         cACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATcg                                                                                                                        >  6:107526/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                                                                                             <  2:142786/149‑1 (MQ=255)
                                                    cGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCt                                                                                             <  1:12855/149‑1 (MQ=255)
                                                      gTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTg                                                                                           <  5:107526/149‑1 (MQ=255)
                                                      gTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTg                                                                                           <  5:143764/149‑1 (MQ=255)
                                                         ccccAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGt                                                                                         <  6:254563/148‑1 (MQ=255)
                                                          cccAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGt                                                                                         >  5:173479/1‑147 (MQ=255)
                                                          cccAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGt                                                                                         >  1:86394/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                                          >  2:253966/1‑148 (MQ=255)
                                                                                                        tCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCa                                          >  7:251625/1‑148 (MQ=255)
                                                                                                          gAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAg                                       >  4:194136/1‑149 (MQ=255)
                                                                                                               tGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAg                                  >  6:86217/1‑149 (MQ=255)
                                                                                                                  cGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACggcggc                                                           >  1:192/1‑121 (MQ=255)
                                                                                                                  cGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACggcggc                                                           <  2:192/121‑1 (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                        <  2:27478/149‑1 (MQ=255)
                                                                                                                         gACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGAc                        <  5:253372/149‑1 (MQ=255)
                                                                                                                         gACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGAc                        >  4:283690/1‑149 (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                  <  3:294535/149‑1 (MQ=255)
                                                                                                                               ggTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAAtt                  >  6:317200/1‑149 (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            <  4:265346/149‑1 (MQ=255)
                                                                                                                                     ggTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCttgtt            <  5:206640/149‑1 (MQ=255)
                                                                                                                                     ggTGCCCTGGACAAGGACACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCTTAGTGAAGCGTCTTTGAGGGAATTCttgtt            >  2:267855/1‑149 (MQ=255)
                                                                                                                                          tCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCctct                             <  4:147733/128‑1 (MQ=255)
                                                                                                                                          tCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCctct                             >  3:147733/1‑128 (MQ=255)
                                                                                                                                            cTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTAt       <  2:324169/148‑1 (MQ=255)
                                                                                                                                            cTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTAt       <  2:91458/148‑1 (MQ=255)
                                                                                                                                            cTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTAt       <  2:319093/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)
                                                                                                                                                  |                                                                                                                                                  
GCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCCCAGTTTGCGGGTTCTTGCCGGTACGCGCACCGCGATGACGTTTTTCGAAGGTGCCGAAACCGACCAGGGTCACGGTGTCCTTGTCCAGGGCACCGGTGATGCTGTCGAGAATCGCGTTCAAAACCTGATTGGCCTTTTCCTTGGTCAGATCAGCCTTATCGGCGATGACGGCGGCGAGTTCTGGTTTGCGCATAGTGAAGCCTCTTTGACGGAATTCTTGTTGTTATGCCGT  >  NC_002947/6061115‑6061406

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: