Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,168,588 G→T 5.4% intergenic (+99/+154) 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,5880GT 5.4% 263.5 / 3.9 94intergenic (+99/+154)eda/leuAKHG/KDPG aldolase/2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base G (48/40);  new base T (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 = 2.64e-01

ATCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCA  >  NC_002947/1168442‑1168737
                                                                                                                                                   |                                                                                                                                                      
aTCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAAAcc                                                                                                                                                     >  2:919161/1‑150 (MQ=255)
 tCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGGAGGccc                                                                                                                                                    >  1:1013501/1‑150 (MQ=255)
      gCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTcc                                                                                                                                               >  1:1190763/1‑150 (MQ=255)
        gTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTa                                                                                                                                             <  1:2230635/150‑1 (MQ=255)
        gTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTa                                                                                                                                             >  1:209726/1‑150 (MQ=255)
             gTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCAc                                                                                                                                        <  2:442709/150‑1 (MQ=255)
             gTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCAc                                                                                                                                        <  2:227722/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                                                                                                                                              >  1:1052587/1‑113 (MQ=255)
                                        ccAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTCTACCCGcgaaa                                                                                                                                              <  2:1052587/117‑5 (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                                                                                                 >  2:1333188/1‑145 (MQ=255)
                                                         tgtgCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGc                                                                                                 <  1:1333188/145‑1 (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                                                                 <  2:773891/111‑1 (MQ=255)
                                                                                                                            ttGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                 >  1:773891/1‑111 (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                                                                                  <  2:1057351/86‑1 (MQ=255)
                                                                                                                                    gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCggg                                                                                  >  1:1057035/1‑86 (MQ=255)
                                                                                                                                    gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCggg                                                                                  >  1:1057351/1‑86 (MQ=255)
                                                                                                                                    gCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCggg                                                                                  <  2:1057035/86‑1 (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)
                                                                                                                                                   gACCCGCTCCTCCGCCCCGATCTGCGGTGGGCACACACGGGGTGTTGGACACAGGCCGTAGGAGGAGCGGGTTTTCCCGTGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCTCATCAAACGATGCGGGCTTTCTCATCCCTCATttca  <  1:910517/150‑1 (MQ=255)
                                                                                                                                                   gACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATttca  >  2:804637/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ATCGAAGCGTGCAGTGCCGAGGCGATGGCACTGCTGGACGCCAACTGAAAACAGTTGTGTGCTTTACGGCTTGTGGGGCGCTTGGTCGGCGCCCCTTTTTTTTGCCTGGGGTATTTGCAGTGGTTGTGCTGGCCCTTTCGCGGGTAGACCCGCTCCTACGCACCGATCTGCGGTGGGCACACACCGGGTGTTGGACACAGGCCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGTCCAGGCCAACAAAAAAAGCCCGCATCAAACGATGCGGGCTTTCTCATCCCTCATTTCA  >  NC_002947/1168442‑1168737

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: