Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,168,587 A→C 5.3% intergenic (+98/+155) eda → / ← leuA KHG/KDPG aldolase/2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,168,5870AC 5.3% 261.7 / 4.4 94intergenic (+98/+155)eda/leuAKHG/KDPG aldolase/2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base A (49/40);  new base C (3/2);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GCGGATCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCA  >  NC_002947/1168438‑1168732
                                                                                                                                                      |                                                                                                                                                  
gcgGATCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTa                                                                                                                                                    >  2:2230635/1‑150 (MQ=255)
    aTCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAAAcc                                                                                                                                                >  2:919161/1‑150 (MQ=255)
     tCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGGAGGccc                                                                                                                                               >  1:1013501/1‑150 (MQ=255)
          gCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTcc                                                                                                                                          >  1:1190763/1‑150 (MQ=255)
            gTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTa                                                                                                                                        >  1:209726/1‑150 (MQ=255)
            gTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTa                                                                                                                                        <  1:2230635/150‑1 (MQ=255)
                 gTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCAc                                                                                                                                   <  2:227722/150‑1 (MQ=255)
                 gTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCAc                                                                                                                                   <  2:442709/150‑1 (MQ=255)
                        ggCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTg                                                                                                                            >  2:2085562/1‑150 (MQ=255)
                           gATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTGTTTTTGCCTGGGGTGTTTGCAGTGGTTGTGCTGGCCCGTTCCGGGAGAGACCCGCCCCTACACCCCGATCAGCgg                                                                                                                         >  1:946928/1‑150 (MQ=255)
                             tGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTg                                                                                                                       <  1:2085562/150‑1 (MQ=255)
                               gCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTggg                                                                                                                     <  1:832547/150‑1 (MQ=255)
                                cACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGc                                                                                                                    >  1:747011/1‑150 (MQ=255)
                                cACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGc                                                                                                                    <  1:644975/150‑1 (MQ=255)
                                 aCTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTTCGGTGGGcc                                                                                                                   >  2:1572831/1‑149 (MQ=255)
                                       ggACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACAcc                                                                                                             <  2:1190763/150‑1 (MQ=255)
                                         aCGCCAACTGAAAGCAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCgg                                                                                                           <  2:901674/150‑1 (MQ=255)
                                            ccAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTCTACCCGcgaaa                                                                                                                                         <  2:1052587/117‑5 (MQ=255)
                                            ccAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTCTACCCGcgaaa                                                                                                                                         >  1:1052587/1‑113 (MQ=255)
                                             cAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGGCTCGCTCCTACGCACCGTTCTGCGGTGGGCACACACCGGgtgt                                                                                                       >  2:1189668/1‑150 (MQ=255)
                                              aaCTTAAAACAGTTGTGTGCTTTACGGCTTGTTGGGCGCTTGGTCGGTGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGTTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGtt                                                                                                      <  2:946928/150‑1 (MQ=255)
                                              aaCTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATc                                                                                                                              <  2:1331151/126‑1 (MQ=255)
                                              aaCTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCGCCTACGCACCGATc                                                                                                                              >  1:1331151/1‑126 (MQ=255)
                                                 tGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGa                                                                                                   <  2:209726/150‑1 (MQ=255)
                                                 tGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGa                                                                                                   <  2:1920615/150‑1 (MQ=255)
                                                    aaaCAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGacac                                                                                                >  2:1075594/1‑150 (MQ=255)
                                                     aaCAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGAc                                                                                                                                                 <  1:1953022/100‑1 (MQ=255)
                                                     aaCAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGAc                                                                                                                                                 >  2:1953022/1‑100 (MQ=255)
                                                     aaCAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGAcaca                                                                                               >  1:350768/1‑150 (MQ=255)
                                                      aCAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAg                                                                                              <  1:1075594/150‑1 (MQ=255)
                                                        aGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGtt                                                                                                      >  1:2107075/1‑140 (MQ=255)
                                                        aGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGtt                                                                                                      <  2:2107075/140‑1 (MQ=255)
                                                          ttGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTa                                                                                                                                        <  1:1638803/104‑1 (MQ=255)
                                                          ttGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTa                                                                                                                                        >  2:1638803/1‑104 (MQ=255)
                                                             tgtgCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGc                                                                                            <  1:1333188/145‑1 (MQ=255)
                                                             tgtgCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGc                                                                                            >  2:1333188/1‑145 (MQ=255)
                                                               tgCTTTACGGCTTGTGGTGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGGCCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCgtagta                                                                                     <  2:2189800/150‑1 (MQ=255)
                                                                 cTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTCTACCCGcgaaa                                                                                                                                         >  2:513455/1‑92 (MQ=255)
                                                                 cTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTCTACCCGcgaaa                                                                                                                                         <  1:513455/96‑5 (MQ=255)
                                                                  tttACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGa                                                                                  >  1:1212641/1‑150 (MQ=255)
                                                                      cGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCTGCTGGTTGTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGACACACACCGGGTGTTGGAAACAGGCCGTAGTAGGAgggg                                                                              >  2:1667381/1‑150 (MQ=255)
                                                                      cGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCgg                                                                              <  2:747011/150‑1 (MQ=255)
                                                                            gtgGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGTGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTAc                                                                        >  2:1446899/1‑150 (MQ=255)
                                                                               gggCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCcg                                                                     >  2:145426/1‑150 (MQ=255)
                                                                                     ttGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTTGCGGGTAGACGCGCTCCTACGCGCCGATCTGCGGTGGGCACACACCGGGTGTTGGACACATGCCTTAGTATGAGCGGGTTTACCCGCGAAAg                                                               >  2:2018390/1‑150 (MQ=38)
                                                                                     ttGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGAGACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAg                                                               >  2:1292374/1‑150 (MQ=255)
                                                                                     ttGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCcgc                                                                    <  1:1332827/145‑1 (MQ=255)
                                                                                     ttGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCcgc                                                                    >  2:1332827/1‑145 (MQ=255)
                                                                                     ttGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGACAg                                                               >  1:583651/1‑150 (MQ=255)
                                                                                     tctgTCGGCGCCCCTTTTTGTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGTGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAg                                                               <  1:2051805/147‑1 (MQ=255)
                                                                                      tGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAgg                                                              <  1:1292374/150‑1 (MQ=255)
                                                                                      tGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAgg                                                              <  1:1581919/150‑1 (MQ=255)
                                                                                       ggTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCCCCACTCTGCGGAGGGAAAGCAGGTACTGTTTGACACAGGCGGTAGTAGGAGCGGGTTTACCCGCGAAAggg                                                             <  1:657604/150‑1 (MQ=255)
                                                                                         tCGGCGCCCCTTTTTTTTTTTTTTTGTTTTTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTTGTAGGATCGGGTTTTCCCGCGAAAGGGc                                                            >  2:1490470/1‑150 (MQ=255)
                                                                                         tCGGCGCCCCTTTTTTTTTTTTTTTTTTTTTTCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCTTACGCACCGATCTGCGGTGTGCACACACCGGGTGTTGGCCACATGCCGTTGTAGTTGCGGTTTTACCAGCGaaatgaca                                                           >  2:2051805/1‑145 (MQ=255)
                                                                                         tCGGCGCCCCTTTTTTTTGCCTGGTGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCg                                                                                                            <  2:1002987/101‑1 (MQ=255)
                                                                                         tCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCg                                                                                                            >  1:1002987/1‑101 (MQ=255)
                                                                                              gcCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTAGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACCCCGAGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTcc                                                      >  1:2302833/1‑150 (MQ=255)
                                                                                                cccTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGTGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAg                                                    <  1:1446899/150‑1 (MQ=255)
                                                                                                cccTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTg                                                                                                     <  1:1001846/101‑1 (MQ=255)
                                                                                                cccTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTg                                                                                                     >  2:1001846/1‑101 (MQ=255)
                                                                                                 ccTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGAGGGTTTACCCGCGAAAGGGCCCGTACAgg                                                  >  1:501539/1‑150 (MQ=255)
                                                                                                  cTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGTCACGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGc                                                 >  1:1501816/1‑150 (MQ=255)
                                                                                                         tttGCCTGGGTTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCcaacaa                                           <  1:1189668/150‑1 (MQ=255)
                                                                                                          ttGCCTGGGGTATTTGCAGTGGTTGGGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACAGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGCCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACaaa                                          <  2:1189888/150‑1 (MQ=255)
                                                                                                              cTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACaaaaaaa                                      >  2:2093587/1‑150 (MQ=255)
                                                                                                              cTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACaaaaaaa                                      >  1:505803/1‑150 (MQ=255)
                                                                                                               tGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAg                                     >  2:169382/1‑150 (MQ=255)
                                                                                                                          cAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCaaa                          <  1:1303742/150‑1 (MQ=255)
                                                                                                                               gTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCAGCAGAAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATg                     >  1:41956/1‑150 (MQ=255)
                                                                                                                                ttGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTCCCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGc                    <  1:1302554/150‑1 (MQ=255)
                                                                                                                                ttGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                            >  1:773891/1‑111 (MQ=255)
                                                                                                                                ttGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                            <  2:773891/111‑1 (MQ=255)
                                                                                                                                ttGTGCGGGCCCTTTGGCGAGTAGAACCGCTCCTACGCACGGATCTGCGGTGGGCACAAACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCAAGGCAAAAAAAAAAAGCCCGCATCAAACGATGc                    <  2:505803/150‑1 (MQ=37)
                                                                                                                                  gtgCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCACTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCgg                  <  2:1391339/150‑1 (MQ=255)
                                                                                                                                   tgCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCggg                 >  2:1549325/1‑150 (MQ=255)
                                                                                                                                   tgCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCAAGGCAAACAACAAAAGCCCGCATCAAACGATGCggg                 >  2:910517/1‑150 (MQ=255)
                                                                                                                                       ggCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCttt             >  1:2054762/1‑150 (MQ=255)
                                                                                                                                        gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCggg                                                                             >  1:1057351/1‑86 (MQ=255)
                                                                                                                                        gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCggg                                                                             <  2:1057035/86‑1 (MQ=255)
                                                                                                                                        gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCggg                                                                             <  2:1057351/86‑1 (MQ=255)
                                                                                                                                        gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCggg                                                                             >  1:1057035/1‑86 (MQ=255)
                                                                                                                                          ccTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTtctc          >  2:1210342/1‑150 (MQ=255)
                                                                                                                                          ccTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTtctc          >  1:1934897/1‑150 (MQ=255)
                                                                                                                                           cTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGCGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCa         >  2:2010871/1‑150 (MQ=255)
                                                                                                                                           cTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCa         >  2:1156750/1‑150 (MQ=255)
                                                                                                                                            tttCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCAt        <  2:1212641/150‑1 (MQ=255)
                                                                                                                                              gagcgGGTCTACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGTCAGTCCAGAACCAACCTGACGGGTCGTATCAAACGATGCCGGCTTTCTCATcc      >  2:804226/3‑150 (MQ=18)
                                                                                                                                                  gggTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGaaa                                                                >  1:1206758/1‑89 (MQ=255)
                                                                                                                                                  gggTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGaaa                                                                <  2:1206758/89‑1 (MQ=255)
                                                                                                                                                  gggTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                            <  1:980047/93‑1 (MQ=255)
                                                                                                                                                  gggTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                            >  2:980047/1‑93 (MQ=255)
                                                                                                                                                  gggTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGCGTTGACCCGCGAAAGGGCCAGTACAGGCCAACCCGACAGACCCGCATCAAACGATGCGGGCTTTCTCATCCCTCa  >  1:540331/1‑150 (MQ=255)
                                                                                                                                                    gTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCGTTTGATGCGGGCtttttttgt       >  2:617300/1‑137 (MQ=37)
                                                                                                                                                    gTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCGTTTGATGCGGGCtttttttgt       <  1:617300/143‑7 (MQ=37)
                                                                                                                                                      |                                                                                                                                                  
GCGGATCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCA  >  NC_002947/1168438‑1168732

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