Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7850TA43.0% 112.8 / 46.4 79intergenic (‑399/+149)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (39/6);  new base A (0/34);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.74e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATC  >  NC_000913/1764636‑1764935
                                                                                                                                                        |                                                                                                                                                      
ttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt                                                                                                                                                        >  1:808194/1‑150 (MQ=255)
ttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt                                                                                                                                                        >  1:523314/1‑150 (MQ=255)
 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                       >  2:548859/1‑150 (MQ=255)
  aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaaa                                                                                                                                                     <  1:1077863/151‑1 (MQ=255)
  aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                      >  2:932487/1‑150 (MQ=255)
   agagATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                      >  2:74637/1‑149 (MQ=255)
   agagATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                      >  1:412213/1‑149 (MQ=255)
     agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                  >  2:67297/1‑151 (MQ=255)
     agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                  >  1:886089/1‑151 (MQ=255)
     agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                  <  2:919462/151‑1 (MQ=255)
      gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CACCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGACGGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                  >  1:348178/1‑150 (MQ=255)
       aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTa                                                                                                                                                 >  2:219443/1‑150 (MQ=255)
        tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                                                                                                                               >  1:1036474/1‑151 (MQ=255)
          ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc                                                                                                                                             <  2:558795/151‑1 (MQ=255)
               aaaaGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGt                                                                                                                                        >  2:697176/1‑151 (MQ=255)
               aaaaGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGt                                                                                                                                        >  2:747711/1‑151 (MQ=255)
               aaaaGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGt                                                                                                                                        >  1:319131/1‑151 (MQ=255)
                aaaGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtt                                                                                                                                       >  2:375495/1‑151 (MQ=255)
                 aaGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCGGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAAAAAGCTAACCCGCCGttt                                                                                                                                      <  1:1092127/151‑1 (MQ=255)
                  aGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtttt                                                                                                                                     >  1:651927/1‑151 (MQ=255)
                   gtccggCGTTCAAAAAAAAGAAACCCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑CCCC‑CCCCACCCGGCTTTTATCCCCGCCCCTCTGAGTGTTTTGCCTTTTGAGGGGGCGAAAAGAATAATCTCAAAAAGCAAACCCGCCGtttt                                                                                                                                     <  1:3245/145‑1 (MQ=255)
                     gCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAac                                                                                                                                  >  1:687886/1‑151 (MQ=255)
                      cAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAac                                                                                                                                  >  2:393920/1‑150 (MQ=255)
                       aGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAacac                                                                                                                                >  2:730908/1‑151 (MQ=255)
                           ttCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAAc                                                                                                                            >  2:347607/1‑151 (MQ=255)
                             cAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTg                                                                                                                          >  1:147045/1‑151 (MQ=255)
                             cAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTg                                                                                                                          >  1:51173/1‑151 (MQ=255)
                             aaaaaaaaaaaaaCCCC‑TTAAGTCGGGGCGACACGGGTTTAAGGGTT‑TCTC‑CCCCACCCGGCGCTTATCCCCGCCACTCTCAGGGGTTTAGCTCTGGAGGGGGCGATAAAAATAATCTCAAAAAGCAACCCCGCCGTTTTAACACAAACTg                                                                                                                          <  2:869249/148‑1 (MQ=255)
                              aaaaaaaaaaaaCCCC‑AAAATTCTGTGCGACAC‑GGTTTAAGGGTTCTCCC‑CCCCATCGGGCGTTTATCCCCGGCACTCTAAGTGGTTTGGCTTTAAAGGGGGCAAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTg                                                                                                                          <  2:393129/148‑1 (MQ=255)
                                   tAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTGTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                     >  1:904556/1‑150 (MQ=255)
                                    aaaaaacCCC‑TTAAGTCTGTGCGACACAGGCTTAAGGTTT‑TCTC‑CCCCACCCGGCGCTTATCCCCGGCACTCTCAGTGGCTTAGCTCTTGAGGGGGCGATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                                   <  2:1013353/148‑1 (MQ=255)
                                    aaGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                     >  2:302574/1‑149 (MQ=255)
                                     aaaaaaCCCTTTAGGTGTGGGCGCCACGGGCTTAGGGGTT‑TCCC‑CCCCCCCCGGCTCTTCTCTCCGCCACTCTCGGTGCTTTCGCTCTTGAAGGGTCAAAAAGAATAATCTGATAAGGTTACCCCGCCGTTTTAACACAAACTGCGAt                                                                                                                      <  2:773331/146‑1 (MQ=255)
                                     aGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                                   >  2:1101915/1‑150 (MQ=255)
                                     aGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGt                                                                                                                  >  1:1053524/1‑151 (MQ=255)
                                       aaacCCC‑TTAAGTCTGTGCGACACAGGCTTAAGGTTT‑TCTC‑CCCCACCCGGCGTTTATCTCCGGCACTCTCAGGGGTTTAGCTTTTGAAGGGGCGAAAAAAATAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGta                                                                                                                 <  2:612145/150‑1 (MQ=255)
                                       aaaaCCCTTTAAGTCTGTGCGACACAGGTTTAAGGTTT‑TCTC‑CCCCACCCGGCTCTTATCCCCGCCACTCTCAGTGGCTTAGCTTTTGAAGGGGAGATAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGta                                                                                                                 <  1:1030414/151‑1 (MQ=255)
                                        aaaCCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatt                                                                                                               >  2:1120963/1‑151 (MQ=255)
                                            ccctaagctCTGGGCGACCGAGACTAAGGGTTT‑TCCC‑CCCCCCCCGGCGCTTACCCCCGGCCCCCTCAGTGGTTTAGCTTTTGAAGGGCGAATAAAAATAATCTAAAAAAGCCACCCCGCCGTTTTAACACAAACTGCGATTAGtattat                                                                                                             <  1:9348/142‑1 (MQ=255)
                                            cc‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                                                                           >  2:294044/1‑151 (MQ=255)
                                               taaaGTCTGTGCGACACAGGCTAAAGGTTT‑TCCC‑CCCCACCCGGCTCTTATCCCCGCCACCCTAAGTGGTTTAGCTTTTGAAGGGGCAAAAAAAATAATCTAAAAAACCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                          <  1:1103925/148‑1 (MQ=255)
                                                 agtccTGTGCGCCACAGGCTAAGGGTTT‑TCTC‑CCCCACCGGGCGTTTATCCCCGCCACTCTCAGTGGTTTAGCTCTTGAGGGGGCAAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                        <  2:5958/146‑1 (MQ=255)
                                                 aattcggtgCCACACGGGTTAAAGGGTT‑TCAC‑CCCCCCCGGGCGTTTATCCCCGCCCCTCTCAGTGCTTTTTCTTTTAAGGGGGAAATAAAAATATTCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                        <  1:983065/144‑1 (MQ=255)
                                                 aaGTCTGTGCGACACAGGTTTAAGGTTT‑TCTC‑CCCCATCCGGCGCTTATCTCCGCCACTCTCAGGGGCTTAGCTCTTGAAGGGGCAATAAGAATAATCTCAAAAATCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGa                                                                                                       <  2:815812/151‑1 (MQ=255)
                                                       gtgCGACCCAGGCTTAGGGGTT‑TCTA‑CCCCATCGGGCGCTTATCCCCGGCACTCTCAGTGTTTTAGCTTTGGAAGGGGCGATAAAAATAATCTAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAata                                                                                                 <  2:586486/151‑1 (MQ=255)
                                                       gtgCGACACAGGCTTAAGGTTT‑TCTC‑CCCCCCCCGGCGCTTATCTCCGGCCCTCTCAGTGGTTTAGCTTTTGAGGGGGCGATAAGAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAat                                                                                                  <  2:151041/150‑1 (MQ=255)
                                                       gtgCGACACAGGCTAAAGGTTT‑TCTC‑CCCCACCCGGCGCTTACCTCCGGCCCCCCCAGTGGTTTGGCTTTTGAAGGGGAAAAAAGAATAATTTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                      <  1:147844/146‑1 (MQ=255)
                                                          cGACACGGCTTTAAGGGTT‑TCTC‑CCCCACCGGGCTTTTACCTCCGCCTCCCCCAGTGGCTAAGTTCTTGAGGGGGCGAAAAAAATAATCTCAAAAAGTTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                               <  2:515107/150‑1 (MQ=255)
                                                          cGACACATGATT‑AGGGTT‑TCCACCCCCCCCCGCCGCTTATCCCCGGCACTCCCAGTGTTTTGGCTTTTAAAGGGGAAATAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCa                                                                                              <  1:76653/151‑1 (MQ=255)
                                                           gACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                               >  2:635310/1‑149 (MQ=255)
                                                           gACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                               >  2:670168/1‑149 (MQ=255)
                                                            acacAGTCTTAAGGTTT‑CCTC‑CCCCACCCGGCTCTTATCCCCGCCCCTCTCAGTGCCTTACCTCTTGAGGGGGCAAAAAAAATAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                                                                    <  2:1023110/143‑1 (MQ=255)
                                                             cacggggtaaaGGGTT‑TCTC‑CCCCATCGGGCTTTTATCCCCGCCTCTTTAAGGGCTTTGGCTTTTAAGGGGTCAAAAAGAATAATCCCAAAAAGTTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACaa                                                                                                   <  2:716043/134‑1 (MQ=255)
                                                              acagtttaaaGGGTT‑TCTC‑CCCCACCCGGCTCTTATCTCCGCCACTCTCAGTGGTTTAGCTTTTGAAGGGGCAAAAAAAATAATCAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAgg                                                                                            <  1:141641/141‑1 (MQ=255)
                                                               caGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg                                                                                         >  2:625038/1‑151 (MQ=255)
                                                                aGGTTTAGGGTTT‑TCTC‑CCCCACCCGGCGCTTATCCCCGGCCCTCTCAGTGGTTTAGCTTTTGAAGGGGCAAAAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCg                                                                                          <  1:558972/149‑1 (MQ=255)
                                                                 ggCTTAAGGGTT‑TCTC‑CCCCACCCGCCGCTTTTCTCCGCCCCTCTCAGTGGTTTAGCTCTTGAAGGGGCGATAAAAAAAATCTAAAAAAGTTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCg                                                                                          <  2:47222/148‑1 (MQ=255)
                                                                    ttAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                                        >  2:307439/1‑147 (MQ=255)
                                                                      aaGGTTT‑TCTC‑CCCCACCCGCCCTTTATCTCCGCCCCTCTCGGTGGTTTAGCTTTGGAAGGGGAAATAAAAATAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATaa                                                                                  <  1:28512/151‑1 (MQ=255)
                                                                        agggttttcccccccCTCCGGCGTTTATCTCCGCCACTCTCAGTGGCTTAGCTTTTGAAGGGGAAAAAAAAATAATCTCATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGa                                                                                     <  2:796364/137‑1 (MQ=255)
                                                                         gttt‑tCTC‑CCCCACCCGGCCCTTATACCCGCCACTCTCGGTGGTTTAGTTTTGGAGGGGGGAAAAAAAAAAATCTAAAAAGGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCa                                                                               <  2:1061853/149‑1 (MQ=255)
                                                                            t‑tCTC‑CCCCACCCGGCCTTTACCCCCGCCCCTCTCAGTGGTTTATCTTTTAAGGGGGCAAAAAAAAAAATTTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt                                                                             <  1:1155623/150‑1 (MQ=255)
                                                                            t‑tCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGAAGt                                                                             >  1:777282/1‑150 (MQ=255)
                                                                                 a‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa                                                                         >  1:1117828/1‑150 (MQ=255)
                                                                                   ccccACCCGGCGCTTCTCCCCGGCCCTCTCATTGGTTTAGCTCTTAAGGGGGCGAAAAGAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt                                                                             <  2:1202615/145‑1 (MQ=255)
                                                                                      cATCCGGCGTTTATCCCCGGCCCTCTCAGTGGTTTAGCTTTTAAAGGGGAGATAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaag                                                                    <  1:263881/151‑1 (MQ=255)
                                                                                        cccGGCGCTTGCCTCGGCCCCTCTCAGTGATTTAGCTTTTGAGGGGGAAAAAAAAAAAATCCAAAAACTCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGt                                                                  <  1:46332/150‑1 (MQ=255)
                                                                                                              ctcAGTGGTTTAGCTCTAAAGGGGGAAATAAAAATAATCAAAAAAAGAAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGt                                              <  1:471022/149‑1 (MQ=255)
                                                                                                                 aGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATc                                         >  1:274062/1‑151 (MQ=255)
                                                                                                                 aGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATc                                         >  1:905437/1‑151 (MQ=255)
                                                                                                                         acgtctTGAAGGGGCGACAAGAATAATCTCATAAAGCTAACGCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGCAGATAATCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                                  >  1:329179/4‑150 (MQ=255)
                                                                                                                         aGCTTTTGAGGGGGAAATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                                  <  2:952596/150‑1 (MQ=255)
                                                                                                                            ttttGAGGGGGCAAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAtt                               <  2:886134/148‑1 (MQ=255)
                                                                                                                               tGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGc                           >  2:422702/1‑151 (MQ=255)
                                                                                                                                   ggggCAAAAAAAATAATCTAAAAAATATAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAAc                       <  1:747681/151‑1 (MQ=255)
                                                                                                                                    gggCGAAAAAAATAACCTAAAAAAAATACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACt                      <  2:1036534/151‑1 (MQ=255)
                                                                                                                                    gggAGATAAGAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACt                      <  1:295125/151‑1 (MQ=255)
                                                                                                                                     ggcaaaaaaaataatCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt                     <  2:891966/142‑1 (MQ=255)
                                                                                                                                        gATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCg                  >  1:722459/1‑151 (MQ=255)
                                                                                                                                                    ctcATAAAGCTAACCCGCCGTTTTAACACAAACTGCGCTTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTc        >  2:643702/1‑149 (MQ=255)
                                                                                                                                                     tcATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAgg     >  1:1004327/1‑151 (MQ=255)
                                                                                                                                                        tAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATc  >  2:181687/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
TTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTTAAGGGTT‑TCTA‑CCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATC  >  NC_000913/1764636‑1764935

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