Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,810,335 T→G 8.0% intergenic (+96/‑73) eno → / → ftsB enolase/cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,810,3350TG 8.0% 219.4 / 8.4 81intergenic (+96/‑73)eno/ftsBenolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base T (32/41);  new base G (3/4);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

AACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGT  >  NC_002947/1810190‑1810482
                                                                                                                                                  |                                                                                                                                                   
aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCgg                                                                                                                                                  >  4:99431/1‑149 (MQ=255)
aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCgg                                                                                                                                                  <  3:187267/149‑1 (MQ=255)
aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGAAGCAGCCGGAGTTGTAGGACCGTTCGTACTGACCTTTCCTACGTAATAACGGGTTTCTGTCGACGAAGCCTGGCCTCgg                                                                                                                                                  >  3:228339/1‑149 (MQ=255)
 aCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGCTTTCTGTCGACGAAGCCTGGCCTCGGc                                                                                                                                                 <  1:16052/149‑1 (MQ=255)
      gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGGCAgg                                                                                                                                             >  1:318516/1‑148 (MQ=255)
      gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAg                                                                                                                                              <  5:15021/147‑1 (MQ=255)
      gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGc                                                                                                                                            <  3:315925/149‑1 (MQ=255)
       cgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCt                                                                                                                                           <  3:371734/149‑1 (MQ=255)
        gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCtt                                                                                                                                          >  6:167116/1‑149 (MQ=255)
         ccAATGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTc                                                                                                                                         <  7:151943/149‑1 (MQ=255)
          cAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGGACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCg                                                                                                                                        >  4:233271/1‑149 (MQ=255)
          cAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGGCAGGCTTCg                                                                                                                                        >  3:5634/1‑149 (MQ=255)
            aaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTg                                                                                                                                      >  4:16458/1‑149 (MQ=255)
                 gTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTAtt                                                                                                                                 <  5:204056/149‑1 (MQ=255)
                         cgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGccc                                                                                                                         <  3:162398/149‑1 (MQ=255)
                            tCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCg                                                                                                                       >  4:225777/1‑148 (MQ=255)
                                gCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGGCTCGGCCAGGCTTCGTGCTATTAGAGGGCCCGaact                                                                                                                   >  3:182621/1‑146 (MQ=255)
                                gCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                                                                                   >  5:68651/1‑148 (MQ=255)
                                  cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTg                                                                                                                                      >  5:340344/1‑127 (MQ=255)
                                  cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTg                                                                                                                                      <  6:340344/127‑1 (MQ=255)
                                   gAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTGGTGCTATTCGAGGCCCCGAAGTAAc                                                                                                                >  7:339194/1‑148 (MQ=255)
                                    aGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACaa                                                                                                              <  5:154005/149‑1 (MQ=255)
                                    aGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACaa                                                                                                              <  4:297109/149‑1 (MQ=255)
                                        tCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGCACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTa                                                                                                           <  3:233271/148‑1 (MQ=255)
                                          gcgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAcgc                                                                                                        >  1:66117/1‑149 (MQ=255)
                                           cgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAcgcg                                                                                                       <  5:304689/149‑1 (MQ=255)
                                                 aGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGc                                                                                                  <  5:250520/148‑1 (MQ=255)
                                                    aaGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCttt                                                                                              <  7:188197/149‑1 (MQ=255)
                                                    aaGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCtt                                                                                               <  8:319338/148‑1 (MQ=255)
                                                      gagaTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCttttt                                                                                            <  1:364982/149‑1 (MQ=255)
                                                        gaTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTaa                                                                                          <  6:32049/149‑1 (MQ=255)
                                                         aTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAc                                                                                         <  2:43579/149‑1 (MQ=255)
                                                          tGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAcc                                                                                        <  3:16458/149‑1 (MQ=255)
                                                             tAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg                                                                                      >  5:373601/1‑148 (MQ=255)
                                                             tAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg                                                                                      >  3:162029/1‑148 (MQ=255)
                                                                aaaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCACGGCCACGCTTCGGGCGATTAGAGGGCCCCCAGTAAAAAGTACGCGGCAGCCTTTTTAAACTGGaa                                                                                   >  8:40039/1‑147 (MQ=255)
                                                                    gcagcaGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGt                                                                              <  5:5878/149‑1 (MQ=255)
                                                                      agcagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                                                                            <  4:241220/149‑1 (MQ=255)
                                                                      agcagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                                                                            <  2:5257/149‑1 (MQ=255)
                                                                      agcagcCGGAGTTGAAGGACCGTTCGTACTTAACTTTCCTACGTACTAACGGGTTTCTGTCGACGATGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGa                                                                            <  2:80702/149‑1 (MQ=255)
                                                                         agcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgc                                                                         >  5:383025/1‑149 (MQ=255)
                                                                          gcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGTTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACGTGGATAACGTGATgcg                                                                        <  7:64804/149‑1 (MQ=255)
                                                                           ccGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgcg                                                                        >  7:324366/1‑148 (MQ=255)
                                                                              gAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGt                                                                    <  3:83021/149‑1 (MQ=255)
                                                                                tacgTAGGAAAGTTCAGTACGAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTcc                                                                  >  3:81573/4‑149 (MQ=255)
                                                                                    tAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                                                              <  6:383025/149‑1 (MQ=255)
                                                                                     aGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                                                              <  7:19965/148‑1 (MQ=255)
                                                                                     aGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCATCGTGCTATTAGAGGCCCCGAAGTAACAAGTACGCGGCAGGCGTTTTAACCTGGTTAACGTGATGCGCAGTCACTAtt                                                             >  8:217116/1‑149 (MQ=255)
                                                                                      ggAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg                                                            <  7:128780/149‑1 (MQ=255)
                                                                                      ggAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg                                                            <  8:324366/149‑1 (MQ=255)
                                                                                            atCGAACAGAACATACCTACGAACTAACGGGTAACTGTCGACGAACCCTGGCCTCGGCCAGGCTACGAGGTATTCCAGGCCCCGAAGTAACAACTACGCGGCCGGCTTTTTAACATGGATAACGTGATGCGCAGTCGCTTTTGggtgg                                                       >  6:342074/2‑147 (MQ=255)
                                                                                                   cTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCg                                                 <  4:5634/148‑1 (MQ=255)
                                                                                                   cTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGt                                                >  5:182557/1‑149 (MQ=255)
                                                                                                    tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGt                                                >  6:83654/1‑148 (MQ=255)
                                                                                                     gAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCGCGGCCCGGCATGGTGCGATTCGAGGCCCCGAAGTAACAAATAACCGGCAGCCATTTTTACCAGGTTAAACTGATGCGCAGTACCGATTTGTTTTTTCACGtta                                              >  4:7133/1‑148 (MQ=255)
                                                                                                      aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGtt                                               <  2:67274/147‑1 (MQ=255)
                                                                                                       aCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTg                                            <  3:370928/149‑1 (MQ=255)
                                                                                                         tttCCTACGTACTAACGGGTTTCTGTCGACGAATCCTGGCCTCGGCCAGGATTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTTACCTGGATAACCTGAAGCGCAGTACCTATTGGTTTTTCATCTTTATGc                                           >  6:78604/1‑148 (MQ=255)
                                                                                                           tCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGctc                                         <  2:262405/148‑1 (MQ=255)
                                                                                                             cTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTg                                      <  2:318516/149‑1 (MQ=255)
                                                                                                              tACGTACTAACGGGTTTCTGTCGACGAAGCCTTGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTTACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCGATTGGTTGTTCCTCGTTCTGCTCTTGc                                     >  7:140192/1‑149 (MQ=255)
                                                                                                              tACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGc                                     >  1:26456/1‑149 (MQ=255)
                                                                                                                     tAACGGGTTTCTGTCGACGAAGCCAGGCCGCGGCCAGGCGTAGTGGGATTAGAGGGCCCCAAGTTACAAATAAGCGGCAGCCGTTTTAAAAGGGTTTAAGTGTTGCGCAGTGCCGATTTGTTTTTACTAGTTATGCTGTTGGTGTTtgg                              >  4:182565/1‑149 (MQ=255)
                                                                                                                       aCGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTg                            >  3:245046/1‑149 (MQ=255)
                                                                                                                           gTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTcc                                                    <  5:364340/121‑1 (MQ=255)
                                                                                                                           gTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTcc                                                    >  6:364340/1‑121 (MQ=255)
                                                                                                                           gTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCt                        <  5:83654/149‑1 (MQ=255)
                                                                                                                             ttCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGGGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCGTTGTAACCTGGATAACGTGATGCGCAGTCCCGATTGGGTGTTCCTCGTTCGGAGCGTGCTGTTGGGGGGGCTGc                      >  6:260299/1‑149 (MQ=255)
                                                                                                                              tCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGc                      <  3:7133/148‑1 (MQ=255)
                                                                                                                                  tcgtctaaGCCTGGCCTCGGCCATGCTTCGTTCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGGGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGTGTGGCCTGCAg                    <  5:78604/140‑1 (MQ=255)
                                                                                                                                     acgaAGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGgtt                                                         >  4:255028/1‑106 (MQ=255)
                                                                                                                                     acgaAGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGgtt                                                         <  3:255028/106‑1 (MQ=255)
                                                                                                                                        aaGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGATGTAACAAGTACGCTGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtg           >  3:48069/1‑149 (MQ=255)
                                                                                                                                             tGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTgg                               <  3:318122/124‑1 (MQ=255)
                                                                                                                                             tGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTgg                               >  4:318122/1‑124 (MQ=255)
                                                                                                                                             tGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGAGATGCGCAGTCCCTATTGGTTGTTACTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggt      >  6:161285/1‑149 (MQ=255)
                                                                                                                                               gCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAAGAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTACTGTTGGGTGGCCTGCAGTACCGCCTgtgggtgg    <  5:342074/149‑1 (MQ=255)
                                                                                                                                               gCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggtgg    <  4:162029/149‑1 (MQ=255)
                                                                                                                                               gCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggtgg    <  4:182621/149‑1 (MQ=255)
                                                                                                                                               gCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggtgg    >  3:126726/1‑149 (MQ=255)
                                                                                                                                               gCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggt      <  5:167116/147‑1 (MQ=255)
                                                                                                                                                  tCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGGTGTTGGGTGGCCTGCAGTACCGCCTGTgggtgggt  >  2:275076/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGT  >  NC_002947/1810190‑1810482

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