Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7850TA40.7% 75.1 / 30.9 59intergenic (‑399/+149)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (28/7);  new base A (0/24);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.23e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGAT  >  NC_000913/1764637‑1764934
                                                                                                                                                       |                                                                                                                                                      
tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATa                                                                                                                                                       <  2:769027/150‑1 (MQ=255)
 aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATaaa                                                                                                                                                     <  1:871981/151‑1 (MQ=255)
  agagATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATaa                                                                                                                                                      >  2:283499/1‑149 (MQ=255)
    agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                  >  1:1234932/1‑151 (MQ=255)
       tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTaa                                                                                                                                                >  1:1007274/1‑150 (MQ=255)
       tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑c                                                                                                                                              >  2:1221346/1‑151 (MQ=255)
         ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑ccc                                                                                                                                            >  2:834403/1‑151 (MQ=255)
          gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑Cccg                                                                                                                                           >  1:939293/1‑151 (MQ=255)
                  gTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGCA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGgttt                                                                                                                                    >  1:828263/1‑150 (MQ=255)
                         gTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACaa                                                                                                                             >  1:224318/1‑150 (MQ=255)
                               aaaaaaaaaaCCCT‑TTAAGTCGGTGCGACACAGGCTA‑AAGGGTTTCTCCCCCACCCGGCCCTTATCCCCGGCCCTCTCAGTGTTTTACCTTTGGAG‑GGGGCAAAAAAAATAATCTAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGa                                                                                                                      <  1:1048269/151‑1 (MQ=255)
                               aaaaaaGAAACCCC‑TTAGGTGTGGGCGACAC‑GGCTTAAAGTGTTCCTCCCCCATCCGTCCTTTATCTCCGCCACTCTCAGTGGATTGGCTCTTGAA‑GGGGCAAAAAAAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTg                                                                                                                         <  1:462145/148‑1 (MQ=255)
                                aaaaaaaaacCCC‑TTAAGTTTGTGCGACACGGGTTA‑AAGGTTTTCTCCCCCACCCGGCGCTTATCCCCGCCACTCTCAGGGGTTTAGCTCTGGAA‑GGGGCGAAAAAAATAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGAt                                                                                                                     <  1:748845/145‑1 (MQ=255)
                                    aaaaaaCCCTTTAAGTCTGTGCGACACAGGCTT‑AAGGTTTTCTACCCCACCCGGCGTTTATCTCCGGCACTCTCAGGGGCTTAGCTCTTGAG‑GGGGCAAAAAAAATAATCTCAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGATTa                                                                                                                   <  1:346579/148‑1 (MQ=255)
                                     aaaacCCC‑TTATTTCTGCGCGACACAGGCTT‑AGGGTTTTCCCCCCCTCCCGGCTTTTATCACCGCCACTCTCGGTGTCTTAGCTTTTAAG‑GGGGCAATAAAAATAATCAAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGta                                                                                                                <  2:318235/150‑1 (MQ=255)
                                      aaccccc‑TTAATTCGGTGCGCCACAGGCTT‑AGGGTTTCCCCCCCCATCCGGCGTTTACCCCCCCCCCTCTGAGGGGCTTAGCTCTTGAG‑GGGGCAAAAAGAATAACCTCAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGtat                                                                                                               <  2:356305/147‑1 (MQ=255)
                                          ccc‑TAAATTCTGTGCGACACAGGCTT‑AGGGTTTTCCCCCCCACCCGGCGTTTATCCCCGGCATTCTCAGTGGCTTAGCTTTTGAG‑GGGGCAAAAAAAATAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGtattat                                                                                                            <  2:647974/150‑1 (MQ=255)
                                          ccc‑TAAAGTTTGGGCAACACAGGCTT‑AAGGTTTCCTCCCCCACCCGGCGTTTATCCCCGCCCCTCTGGGTGTTTTAGCTTTTGAG‑GGGGCAAAAAAAATAATCTAAAAAACCTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGtatta                                                                                                             <  2:3551/149‑1 (MQ=255)
                                           cc‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGCA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                                                                          >  1:903065/1‑151 (MQ=255)
                                            c‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                         >  1:1254548/1‑151 (MQ=255)
                                                 gttcggtgAGACACCGCCTA‑AAGGTTTTCTCCCCCAACCGTTTTTAATCTCCGCCTCTCTCAGTGGTTTAGCTCTGGAG‑GGGGCAAAAAAAATATTCACATAAAGCAAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGa                                                                                                      <  1:169686/145‑1 (MQ=255)
                                                    cTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                                                                   >  1:386986/1‑150 (MQ=255)
                                                    cTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAAATGCGATTAGTATTATTTTTGAACa                                                                                                   >  1:520142/1‑150 (MQ=255)
                                                      gtgCGACACGGGTTA‑AAGGGTTCCTCCCCCACCCGGCGCTTACCCCCGCCCCTCCCAGTGCTTTAGCTGTGGAG‑GGGGCAAAAAAAAAAATCTAAAAAAGCAAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                     <  1:342109/146‑1 (MQ=255)
                                                      gtgCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAat                                                                                                 >  1:512356/1‑150 (MQ=255)
                                                             acaGGCTT‑AGGGTTTTCTCCCCCACCCGCCCCTTATCCCCGCCACTCTAAGTGTCTTAGTTTTTGAA‑GGGGAAAAAAAAATAATCTAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                          <  2:703327/150‑1 (MQ=255)
                                                              caGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg                                                                                        >  2:17418/1‑151 (MQ=255)
                                                                ggttt‑AGGGGTTCCTCCCCCACCCGGCCTTTATCCCCGCCCCCCCCGGTGGTTTAGTTCTTAAG‑GGGGAGATAAGAATAATCTAAAAAAGTAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTa                                                                                      <  1:1039619/148‑1 (MQ=255)
                                                                   tt‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                                       >  1:8934/1‑147 (MQ=255)
                                                                      aaGGGTTTCTCCCCCACCCGGCGTTTATCCCCGGCACTCTGAGTGGCTTAGCTTTTGAG‑GGGGCGATAAAAATATTCTCAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATa                                                                                  <  1:856077/150‑1 (MQ=255)
                                                                      aaGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATa                                                                                  >  2:473140/1‑150 (MQ=255)
                                                                        gtttttCCACCCCACCCGGCGTTTATCCCCGCCCCTCCAAGGGGCTTAGCTCTGAAA‑GGGGAAATAAAAATAATCTAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGc                                                                               <  2:1255505/148‑1 (MQ=255)
                                                                           ttcccccccctcccgccgCTTTTCCCCGGCACCCACGGTGGATTATTTCTGGAA‑GGGGCTAAAAGAATAACCTCAAAAAGACAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAg                                                                                <  1:828679/135‑1 (MQ=255)
                                                                              cccccccACCGGGCTTTTACCTCCGGCACTCCTGGGGGTTTAGCTCTTGAG‑GGGGCAAAAAAAAAAATCAAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt                                                                         <  2:684429/148‑1 (MQ=255)
                                                                                  cccACCCGGCTCTTTCCCCCGCCCCTCTCGGTGGTTTACTTTTTGAA‑GGGGCAATAAAAAAAATCT‑AAAAAACTAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga                                                                     <  2:964081/151‑1 (MQ=255)
                                                                                     aTCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAgg                                                                  >  1:152021/1‑151 (MQ=255)
                                                                                     aTCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAgg                                                                  >  1:374300/1‑151 (MQ=255)
                                                                                       ccGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                                >  1:364274/1‑151 (MQ=255)
                                                                                        cGGCGCTTATCCCCGCCCCTCCCGGTGGTTTTTTTCTTGAG‑GGGGAAAAAAAAATAATCTAAAAAAGCAAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaag                                                                   <  2:365571/147‑1 (MQ=255)
                                                                                         gcctcTTATCCCCGCCACTCACAGTGTTTTTGCTTTTGGG‑GGGGTGAAAAAAATAATCAAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                                <  2:758000/145‑1 (MQ=255)
                                                                                           cgcTTATCTCGGGCACTCTAGGTGCTTTAGCTTTTGAG‑GGGGCAAAAAAAATAACCTAAAAAACCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                                <  1:753727/147‑1 (MQ=255)
                                                                                            gcTTATCTCCGGCACTCCCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCg                                                            >  2:884059/1‑150 (MQ=255)
                                                                                                  ctcCGGCTCGCTCATTGTTTTAGCT‑TTGAAGGGGGCAAAAAAAATAATCTCAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCg                                                            <  1:970327/144‑1 (MQ=255)
                                                                                                   tcCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTgg                                                     >  1:103211/1‑150 (MQ=255)
                                                                                                   cccgcccccctctGTGTCTTAGCTTTAGAG‑GGGGCAAAAAGAATAATCTAAAAAACCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGa                                                    <  2:663855/142‑1 (MQ=255)
                                                                                                    ccgcccctctCAGGGTTTTAGTTCTTGAA‑GGGGAAAAAAAAATAATCTAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAAc                                                         <  2:817803/139‑1 (MQ=255)
                                                                                                    ccGGCCCTCTCAGTGGCTTAGCTTTTGAA‑GGGGCAATAAAAATAATCTAAAAAAACTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTg                                                      <  2:683411/148‑1 (MQ=255)
                                                                                                       cccctctCAGTGTTTTAGCTTTTGAA‑GGGTTAATAAGAATATTCTCAAAAGGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACgg                                                 <  2:158369/147‑1 (MQ=255)
                                                                                                                gTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                       >  1:361375/1‑151 (MQ=255)
                                                                                                                gTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGGCATCGAACCTGGACGGAGGTTAAt                                         >  1:391080/1‑149 (MQ=255)
                                                                                                                 tGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                       >  2:263728/1‑150 (MQ=255)
                                                                                                                  ggtttAGCTCTTAAG‑GGGGCAAAAAAAATAATCTCAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCa                                      <  1:202296/147‑1 (MQ=255)
                                                                                                                     ttAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGCTAATCCAgg                                    >  2:833002/1‑149 (MQ=255)
                                                                                                                                  ggggCAATAAAAATAATCTAATAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGa                        <  2:917695/149‑1 (MQ=255)
                                                                                                                                        aTAAAAATAATTTAATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGc                  <  2:512334/149‑1 (MQ=255)
                                                                                                                                          aaaaataatCTAATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCg                 <  1:976387/145‑1 (MQ=255)
                                                                                                                                               taatCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTa           >  2:110615/1‑149 (MQ=255)
                                                                                                                                                actctcATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAg         >  2:1127450/3‑150 (MQ=255)
                                                                                                                                                      aTAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGAt  >  2:997922/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCC‑TTAAGTCTGTGCGACACAGGCTT‑AAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGAT  >  NC_000913/1764637‑1764934

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