Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,7190GA21.6% 187.9 / 17.9 88intergenic (+70/+32)zntA/tusAZn(2(+))/Cd(2(+))/Pb(2(+)) exporting P‑type ATPase/sulfur transfer protein TusA
Reads supporting (aligned to +/- strand):  ref base G (47/22);  new base A (0/19);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 2.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTA  >  NC_000913/3608582‑3608853
                                                                                                                                          |                                                                                                                                       
tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                      >  1:317941/1‑140 (MQ=255)
      cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                                >  1:985405/1‑140 (MQ=255)
      cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                                >  1:14654/1‑140 (MQ=255)
       aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcgg                                                                                                                               >  1:1465942/1‑140 (MQ=255)
               gggCGAGGGTCTTGGTGAAAGGGAAGGTGTTAAGATTGTCGCGAAGGAAATAAGAAAAACCGATCGAAACATTGAGCGCGATCGGTCCCCCCGCCCCTTTGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                       <  2:2030905/140‑1 (MQ=25)
               gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAGACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                       >  2:692429/1‑140 (MQ=255)
                    ggggTGCTGGTGACAGCGATTGCGTTAAGTTTGTTGCGCAGGAGAAAAGGCAACCCGATCGCAACTTTGAGCGCGATCGTCCCCCTCGCCCCTTTGGGGAGGGGGTTGGGGTGAGGGAAAAAGGCGGCATCGAAGCCAAt                                                                                                                  <  2:492524/138‑1 (MQ=255)
                    gcGGTGCGGGTGACAGCGAATGCTTAAAGTTTGTGGCGCGGGAAAAAAGGCAACCCGATCGCAACATTGAGCGCGATCGGCCCCCCCCCCCCTCGGGGGAGAGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAAt                                                                                                                  <  2:14654/139‑1 (MQ=255)
                    aCGGTGCTGGTAACAGCGAATGCGTAAAGATTGTTGCGCAGGAGAAAAGAAAACCCGATCGAAACATGGAGCGCGATCGCCCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAAt                                                                                                                  <  2:2128167/140‑1 (MQ=37)
                     cGGTGCTGGAGCCGGCAAATGCGTAAAGATTGTTCCGCAGGAAATAAGGCAACCCGATCGCAACTTTGAGCGGGATCGGCCCCCCCCCCCCTCGGGGGAGAGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                                 <  2:1465942/140‑1 (MQ=255)
                       gTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGAGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAg                                                                                                               <  2:915206/140‑1 (MQ=255)
                        tGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                              <  2:1272900/140‑1 (MQ=255)
                           tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                           >  2:1022868/1‑140 (MQ=255)
                                aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                                      >  2:1089299/1‑140 (MQ=255)
                                  aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                    >  1:343810/1‑140 (MQ=255)
                                  aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGCTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                    >  2:605658/1‑140 (MQ=255)
                                    cGAATGCTTAAAGTTTCTGGCGAAGGAGATAAGGAAAACGGTGCGCAACTTGGAGGGCGCGCGGCCCCCCCCCCCTTCGGGGGAGGGGGTGGGGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                                  <  2:343810/140‑1 (MQ=255)
                                       aTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                               >  2:1970980/1‑140 (MQ=255)
                                        tGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                              >  1:1164702/1‑140 (MQ=255)
                                              aaGATTGTGGCGCAGGAGATAAGGCAAACCGATCGCAACATTGGGGGCGATCGGCCCCCCCGCCCCTCGGGGGGGGGGGTTAGGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                        <  2:1997389/140‑1 (MQ=255)
                                                   tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaaa                                                                                   >  1:1305950/1‑140 (MQ=255)
                                                   tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaaa                                                                                   >  1:871930/1‑140 (MQ=255)
                                                    tttgGGCAGGAAAAAAGGAAACCCGATCGCAAATTGGAGCGGGACGGGCCCCCCCCCCCCTTGGGGGAGGGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                  <  2:792625/139‑1 (MQ=255)
                                                    tttgCGCAGGAGATAAGGAAACCCGATCGCAACATTGAGCGCGATCGGCCCCCTCCCCCCTCGGGGGAGGGGGTTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                  <  2:1701539/139‑1 (MQ=255)
                                                    gttgGGAAGGAGAAAAGGAAACCGGTTCGCAAAATTGAGCGCGATCGGCCCCCCCGCCCTTTGGGGGAGGGGTTTGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                  <  2:1432776/140‑1 (MQ=255)
                                                    gttgCGGAGGAAAAAAGGAAAACCGACCGCAACATGGAGCGCGACCGGCCCCCCCCCCCCTGGGGGGAGAGGGTTAGGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                  <  2:1305950/140‑1 (MQ=255)
                                                    gttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                  >  1:1303093/1‑140 (MQ=255)
                                                    gttgCGCAGGAGATAAGGCAAACCGATCGAAACATTGAGCGCGATCGGCCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                  <  2:298024/140‑1 (MQ=255)
                                                     ttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                 >  1:1063991/1‑140 (MQ=255)
                                                        cgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACg                                                                              >  2:996936/1‑140 (MQ=255)
                                                           aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                           >  1:1701539/1‑140 (MQ=255)
                                                            ggaaaaaaGGAAACCCGACCGAAACATGGGGCGCGCTCGGCCCCCCCCCCCCTGGGGGGAGGGGTTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                          <  2:1164702/134‑1 (MQ=255)
                                                              agaTAAGGCAACCCGATCGAAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                        <  2:1063991/140‑1 (MQ=255)
                                                              aaaTAAGGCAAACCGATCGCAACATTGAGCGGGATCGGCCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                        <  2:1368249/138‑1 (MQ=39)
                                                                       aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                               >  2:1028982/1‑140 (MQ=255)
                                                                       aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                               >  1:1262399/1‑140 (MQ=255)
                                                                          ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                            >  1:1825727/1‑140 (MQ=255)
                                                                          ccGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                            >  1:970856/1‑140 (MQ=255)
                                                                            gaccGCAAAATGGAGCGGGATGGGCCCCCCCGCCCCTGGGGGGAGGGGGTTAGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                          <  1:845963/137‑1 (MQ=255)
                                                                              tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                        >  1:809169/1‑140 (MQ=255)
                                                                               cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                       >  1:549084/1‑140 (MQ=255)
                                                                               cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGCGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                       >  1:2118029/1‑140 (MQ=255)
                                                                                 cAACATTGAGCGGGTTCGGCCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                     <  1:2010857/140‑1 (MQ=255)
                                                                                    cATTGGGCGCGATCGGCCCCCTCGCCCCTCGGGGGAGGGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                  <  2:549084/140‑1 (MQ=255)
                                                                                    atttGGGCGCGATCGCCCCCCCCCCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                                  <  2:687172/138‑1 (MQ=255)
                                                                                     aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                                 >  2:2384301/1‑140 (MQ=255)
                                                                                         aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                             >  1:186278/1‑140 (MQ=255)
                                                                                          gcggGATCGGCCCCCCCCCCCCTGGGGGGGGGGGGTTGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                            <  2:695019/140‑1 (MQ=255)
                                                                                           cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                           >  1:749344/1‑140 (MQ=255)
                                                                                            gggaccGGCCCCCTCGCCCCTCGGGGGAGGGGGTTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                          <  2:970856/135‑1 (MQ=255)
                                                                                            gcgACCGGCCCCCTCGCCCCTCGGGGGGGGGGGTTGGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                          <  2:1825727/140‑1 (MQ=255)
                                                                                               ttCGGCCCCCTCGCCCCTGGGGGGAGAGGTTTAGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                       <  2:128891/139‑1 (MQ=255)
                                                                                               gcggacccccTCGCCCCTTGGGGGGGGGGGTTGGGTGAGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                       <  2:2440803/134‑1 (MQ=255)
                                                                                               ctcgccccccccccccCTGGGGGGAGGGGGTTGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                       <  2:749344/134‑1 (MQ=255)
                                                                                               aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                       >  1:695019/1‑140 (MQ=255)
                                                                                               aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                      >  1:1072383/1‑140 (MQ=255)
                                                                                                tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                      >  1:1022117/1‑140 (MQ=255)
                                                                                                tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                      >  1:1177982/1‑140 (MQ=255)
                                                                                                tCGGTCCCCTCGCCCCTCTGGGGAGAGGGCTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGGCCTTCCGTCTCTTTAGCAACCAGTTCgtgt                                      >  1:658614/1‑140 (MQ=255)
                                                                                                       cccccccccTTTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                               <  2:2419661/135‑1 (MQ=255)
                                                                                                       ccTGGCCCCTCGGGGGGGGGGGTAGGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                               <  1:1095787/140‑1 (MQ=255)
                                                                                                         tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                             >  1:1035052/1‑140 (MQ=255)
                                                                                                            ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAggct                          >  1:2440803/1‑138 (MQ=255)
                                                                                                             cccTCGGGGGAGGGGGTTGGGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                         <  1:1022868/140‑1 (MQ=255)
                                                                                                              ccccgggggAGAGGGTTGGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                        <  2:2118029/135‑1 (MQ=255)
                                                                                                              ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                        >  1:1124094/1‑140 (MQ=255)
                                                                                                              ccTCGGGGGAGGGGTTGGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                        <  2:1303093/140‑1 (MQ=255)
                                                                                                                tctGGGGAGAGGGTTAGGGTGAGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                       >  1:2480785/1‑140 (MQ=255)
                                                                                                                tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                       >  1:2564573/1‑139 (MQ=255)
                                                                                                                tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                       <  2:2564573/139‑1 (MQ=255)
                                                                                                                tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                        >  1:1629302/1‑138 (MQ=255)
                                                                                                                tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                        <  2:1629302/138‑1 (MQ=255)
                                                                                                                      gggggggTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                <  2:1262399/136‑1 (MQ=255)
                                                                                                                         agGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa             >  2:2494658/1‑140 (MQ=255)
                                                                                                                           ggTTGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata           <  2:186278/140‑1 (MQ=255)
                                                                                                                           ggTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                       <  2:1944685/128‑1 (MQ=255)
                                                                                                                           ggTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                       >  1:1944685/1‑128 (MQ=255)
                                                                                                                           ggTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata           >  1:2543391/1‑140 (MQ=255)
                                                                                                                            gTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                        >  2:1978940/1‑126 (MQ=255)
                                                                                                                            gTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                        <  1:1978940/126‑1 (MQ=255)
                                                                                                                              tAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg        <  2:1726019/140‑1 (MQ=255)
                                                                                                                              tAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg        <  1:1930282/140‑1 (MQ=255)
                                                                                                                              tAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg        >  2:1930282/1‑140 (MQ=255)
                                                                                                                              tAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg        >  1:1726019/1‑140 (MQ=255)
                                                                                                                              tAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg        <  1:2384301/140‑1 (MQ=255)
                                                                                                                              tAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                       >  1:1406693/1‑124 (MQ=255)
                                                                                                                              tAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                       <  2:1406693/124‑1 (MQ=255)
                                                                                                                               aGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                   >  1:1021931/1‑128 (MQ=255)
                                                                                                                               aGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                   <  2:1021931/128‑1 (MQ=255)
                                                                                                                                   gggggggaAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCggg    <  2:1035052/136‑1 (MQ=255)
                                                                                                                                    gggggggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGt   <  1:605658/137‑1 (MQ=255)
                                                                                                                                    gggggggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTa  <  2:1072383/137‑1 (MQ=255)
                                                                                                                                    gAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTa  <  2:1022117/140‑1 (MQ=255)
                                                                                                                                    gAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTa  <  1:2494658/140‑1 (MQ=255)
                                                                                                                                     aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa              <  2:791832/127‑1 (MQ=255)
                                                                                                                                     aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa              >  1:791832/1‑127 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTA  >  NC_000913/3608582‑3608853

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