Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,984,087 A→C 5.3% intergenic (+37/+153) ilvE → / ← PP_3512 branched‑chain‑amino acid aminotransferase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,984,0870AC 5.3% 302.0 / 5.2 104intergenic (+37/+153)ilvE/PP_3512branched‑chain‑amino acid aminotransferase/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (51/47);  new base C (3/3);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGG  >  NC_002947/3983940‑3984231
                                                                                                                                                   |                                                                                                                                                  
cACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGACACGCATCAGCTGTATGAAAACGGGGTGGGCCAGCAATg                                                                                                                                                  <  2:3050512/150‑1 (MQ=255)
 aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATgg                                                                                                                                                 <  2:1821200/150‑1 (MQ=255)
    aCCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCAt                                                                                                                                              >  1:110492/1‑150 (MQ=255)
    aCCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCAt                                                                                                                                              <  2:1978802/150‑1 (MQ=255)
     ccTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGAAGCGGCAGGCTGGATCGTCAAGGTCGCATGAAAAGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATg                                                                                                                                             >  1:1452680/1‑150 (MQ=255)
        ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGTGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGccc                                                                                                                                          >  2:1776220/1‑150 (MQ=255)
        ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGccc                                                                                                                                          >  2:2822811/1‑150 (MQ=255)
        ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGccc                                                                                                                                          >  2:1669948/1‑150 (MQ=255)
              ggTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGtttt                                                                                                                                    <  2:1696182/150‑1 (MQ=255)
               gTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                                  >  2:1504004/1‑135 (MQ=255)
               gTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                                  <  1:1504004/135‑1 (MQ=255)
                 cGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTc                                                                                                                                 <  1:2822811/150‑1 (MQ=255)
                  ggCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCa                                                                                                                                >  1:485328/1‑150 (MQ=255)
                   gCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTcct                                                                                                                               >  1:120816/1‑148 (MQ=255)
                   gCCCGGTTACCCAGAAGATATACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAt                                                                                                                               >  2:2574776/1‑150 (MQ=255)
                     ccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAtgt                                                                                                                             <  1:2911982/150‑1 (MQ=255)
                       ggTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATgtgt                                                                                                                           <  1:1776220/150‑1 (MQ=255)
                        gTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGtt                                                                                                                          <  1:2574776/150‑1 (MQ=255)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATGAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCTCCGTTTTTTCATGTGTTTataa                                                                                                                     >  2:2552926/1‑147 (MQ=255)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGtt                                                                                                                                      <  2:2159913/133‑1 (MQ=39)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGtt                                                                                                                                      >  1:2159913/1‑133 (MQ=39)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTacac                                                                                                                     >  1:1020642/1‑150 (MQ=255)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATCCCCCGTTTTTTCATGTGtttcccc                                                                                                                     >  2:137716/1‑146 (MQ=255)
                              cAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACg                                                                                                                    >  2:654253/1‑150 (MQ=255)
                                gAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGc                                                                                                                  <  1:1669948/150‑1 (MQ=255)
                                gAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGc                                                                                                                  <  2:1466136/150‑1 (MQ=255)
                                  aGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGccc                                                                                                                >  1:1365039/1‑150 (MQ=255)
                                      ctACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTc                                                                                                            <  2:120816/150‑1 (MQ=255)
                                                cTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCg‑‑g                                                                                                <  2:1020642/150‑1 (MQ=255)
                                                  gatcGGTATCCAGAGAGAAGATGTAGAGGCACCGGCAGGCTGTATCGTCAAGGTCGGCTGAAACGCATCAGCTGTATGAAAAAGGGGAATGCCAGCAATGGCATGCTGCGTGTTTTGATGTGTTTATAGGGCGCCGTCGTAGGAGCG‑‑GGt                                                                                              >  2:1573587/4‑150 (MQ=255)
                                                       gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTAcc                                                                                         <  2:2793000/150‑1 (MQ=255)
                                                        tATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTAccc                                                                                        <  1:654253/150‑1 (MQ=255)
                                                          tCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCcgc                                                                                      <  1:173936/150‑1 (MQ=255)
                                                           ccAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCg‑‑g                                                                                                <  2:1249069/139‑1 (MQ=255)
                                                           ccAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCg‑‑g                                                                                                >  1:1249069/1‑139 (MQ=255)
                                                            cAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGc                                                                                                                                                >  1:1392580/1‑92 (MQ=255)
                                                            cAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGc                                                                                                                                                <  2:1392580/92‑1 (MQ=255)
                                                             agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGaa                                                                                   <  2:2490199/150‑1 (MQ=255)
                                                             agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGaa                                                                                   >  1:1458936/1‑150 (MQ=255)
                                                             agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGaa                                                                                   <  1:70787/150‑1 (MQ=255)
                                                                  ggcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGc                                                                              <  2:110492/150‑1 (MQ=255)
                                                                    cgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCca                                                                            <  2:1452680/150‑1 (MQ=255)
                                                                       tttAGAGGCACCGGCAGGCTGGACCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCacac                                                                         <  1:819622/148‑1 (MQ=255)
                                                                           gagGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCgc                                                                     >  2:3058867/1‑150 (MQ=255)
                                                                           gagGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTCTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCgc                                                                     >  2:343557/1‑150 (MQ=255)
                                                                              gCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAt                                                                                                                               <  1:853988/91‑1 (MQ=255)
                                                                                aCCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCca                                                                            >  2:3047644/1‑138 (MQ=255)
                                                                                aCCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCca                                                                            <  1:3047644/138‑1 (MQ=255)
                                                                                aCCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCgccagcc                                                                >  2:397103/1‑150 (MQ=255)
                                                                                  cGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGc                                                                              >  1:482502/1‑134 (MQ=255)
                                                                                  cGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGc                                                                              <  2:482502/134‑1 (MQ=255)
                                                                                  cGGCAGGCTGGATCGTCAAGGTCGCCGGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCACGGGCATGCCCCGTTTTTTCATGTGTTCACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCAAACCCGCCAGCCaa                                                              <  2:537969/150‑1 (MQ=255)
                                                                                    gCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAtg                                                                                                                              <  2:2232961/86‑1 (MQ=255)
                                                                                    gCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAtg                                                                                                                              >  1:2232961/1‑86 (MQ=255)
                                                                                    gCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGt                                                                                              >  2:2131274/1‑116 (MQ=255)
                                                                                    gCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGt                                                                                              <  1:2131274/116‑1 (MQ=255)
                                                                                        gCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGc                                                        >  1:1538307/1‑150 (MQ=255)
                                                                                                 tCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTgg                                                    <  1:1698429/145‑1 (MQ=255)
                                                                                                 tCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTgg                                                    >  2:1698429/1‑145 (MQ=255)
                                                                                                      gTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAtgt                                                                                                                             <  1:1515753/69‑1 (MQ=255)
                                                                                                      gTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAtgt                                                                                                                             >  2:1515753/1‑69 (MQ=255)
                                                                                                      gTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCg                                          >  2:853954/1‑150 (MQ=255)
                                                                                                       tCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCg‑‑g                                                                                                >  2:1927699/1‑95 (MQ=255)
                                                                                                       tCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCg‑‑g                                                                                                <  1:1927699/95‑1 (MQ=255)
                                                                                                          ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGc                                                                              >  1:2938338/1‑110 (MQ=255)
                                                                                                          ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGc                                                                              <  2:2938338/110‑1 (MQ=255)
                                                                                                          ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGt                                      <  2:485328/150‑1 (MQ=255)
                                                                                                          ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGt                                      >  2:2390775/1‑150 (MQ=255)
                                                                                                              aaaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATgtgt                                                                                                                           >  1:1959882/1‑63 (MQ=255)
                                                                                                              aaaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATgtgt                                                                                                                           <  2:1959882/63‑1 (MQ=255)
                                                                                                              aaaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGt                                  >  1:800713/1‑150 (MQ=255)
                                                                                                                 cGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCa                               <  2:800713/150‑1 (MQ=255)
                                                                                                                   cATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAgg                                                                               >  1:537358/1‑100 (MQ=255)
                                                                                                                   cATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAgg                                                                               <  2:537358/100‑1 (MQ=255)
                                                                                                                    aTCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGcc                            <  1:397103/150‑1 (MQ=255)
                                                                                                                    aTCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCAAGCGCCGGGTGTTGATGTGTTTACACGGCCCCCTCGTCCGACCGCAGCCATGGCCGTCAAAAG‑‑CACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTTTCAGTGCAGcc                            >  1:701601/1‑150 (MQ=255)
                                                                                                                     tCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCg                           >  2:2211368/1‑150 (MQ=255)
                                                                                                                           atgaaaaaaCGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                  >  1:627785/6‑121 (MQ=37)
                                                                                                                           atgaaaaaaCGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                  <  2:627785/116‑1 (MQ=37)
                                                                                                                            tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGc                    <  2:2467614/150‑1 (MQ=255)
                                                                                                                              aaaaaaCGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                  >  2:2408724/3‑150 (MQ=39)
                                                                                                                               gAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGc                                                        >  1:1434848/1‑111 (MQ=255)
                                                                                                                               gAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGc                                                        <  2:1434848/111‑1 (MQ=255)
                                                                                                                               gAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCggg                 >  2:1474579/1‑150 (MQ=255)
                                                                                                                               gAAAACGGGGCATGCCAGCAATCGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCCAAAAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCggg                 <  1:2552926/150‑1 (MQ=255)
                                                                                                                                aaaaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGc                >  2:1464936/1‑150 (MQ=255)
                                                                                                                                     ggggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                  <  2:2380362/143‑1 (MQ=255)
                                                                                                                                     ggggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                  >  1:2380362/1‑143 (MQ=255)
                                                                                                                                      gggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGt                                      >  2:790075/1‑122 (MQ=255)
                                                                                                                                      gggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGt                                      <  1:790075/122‑1 (MQ=255)
                                                                                                                                      gggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc          <  1:583884/150‑1 (MQ=255)
                                                                                                                                       ggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcg         >  2:2824727/1‑150 (MQ=255)
                                                                                                                                        gCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcgc        <  2:1538307/150‑1 (MQ=255)
                                                                                                                                          aTGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCgcgc      >  1:2691313/1‑150 (MQ=255)
                                                                                                                                          aTGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCgcgc      >  1:1239697/1‑150 (MQ=255)
                                                                                                                                            gCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                  >  1:2309384/1‑136 (MQ=255)
                                                                                                                                             ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAg                                   >  2:1112513/1‑118 (MQ=255)
                                                                                                                                             ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAg                                   <  1:1112513/118‑1 (MQ=255)
                                                                                                                                             ccAGCAATGGCATGCCCCGTCTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg   <  1:343557/150‑1 (MQ=255)
                                                                                                                                              cAGCAATGGCATGCCCCGTTTTTTTTTTTTTTTTTACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGTGTATAGGCCATACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCTGGCAAGCCCGCGCCggg  >  2:1828216/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCG‑‑GGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGG  >  NC_002947/3983940‑3984231

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