Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,036,9060TA47.4% 64.9 / 33.6 57L53* (TTA→TAA) hyaFhydrogenase‑1 protein nickel incorporation factor
Reads supporting (aligned to +/- strand):  ref base T (26/4);  new base A (1/26);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 6.40e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGC  >  NC_000913/1036756‑1037053
                                                                                                                                                        |                                                                                                                                                   
aaaCTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGtt                                                                                                                                                     <  1:12507/151‑1 (MQ=255)
 aaCTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGtt                                                                                                                                                     >  2:112731/1‑150 (MQ=255)
               tgctgGGGCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTacaac                                                                                                                                       >  2:650230/1‑150 (MQ=255)
                gctgGGGCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTacaaca                                                                                                                                      >  2:1039698/1‑150 (MQ=255)
                     ggCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACt                                                                                                                                <  1:513380/151‑1 (MQ=255)
                      gCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GCGTATTGCGGCCCTGGAGCAA‑TGTGATCACACCCCGCAAGGGATTTCGCTGTAAAACGGGTTacaac                                                                                                                                       >  2:106645/1‑143 (MQ=255)
                                cAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGcc                                                                                                                     >  1:436313/1‑151 (MQ=255)
                                 aaCCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGcc                                                                                                                     >  1:281031/1‑150 (MQ=255)
                                   ccGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCTACTAAGCGAACGCCAAc                                                                                                                  <  2:638776/151‑1 (MQ=255)
                                     gAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAAccg                                                                                                                >  1:1036813/1‑151 (MQ=255)
                                         gATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAAccgccg                                                                                                             >  2:198825/1‑150 (MQ=255)
                                         gATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGt                                                                                                            >  2:831501/1‑151 (MQ=255)
                                                    cAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTg                                                                                                 >  2:207600/1‑151 (MQ=255)
                                                       cATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTg                                                                                              >  1:528196/1‑151 (MQ=255)
                                                        ataactccccggCCGCTCCCCTGTCAGGTAAATGATGACCC‑GAGTTTGGGGGCCGTGGAGAAG‑TGTCCCCACAGCCCGGGGGTGATGGCGTTGTTAAACGGGTTACAACATAAACTAAGCGAACCCCAACCGCCGTTGGGCGAGGTGCTg                                                                                              <  1:184931/139‑1 (MQ=255)
                                                          gAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGc                                                                                            >  1:566477/1‑150 (MQ=255)
                                                           aacccccggCCGATCCCGTGTCGGGTAAATGATAAACC‑GAGTAGGGCGCCCGGGGAGAAA‑TGTCCCCCCAGCCCGGGGGTGATGGCGTTGTAAAACGAGTAACAACAAAAACAAAGCGAACGCAAACCGCCGTTGGGCGAGGTGCTgg                                                                                             <  2:714993/140‑1 (MQ=255)
                                                           aacccccTGCCGTCCCCCTGTCGGGGAAATGATGAACC‑GAGTAGGGGGCCCGTGGAGAAA‑GGTCCTCACACCCCGCAGGTGATTGCCTTGTAAAACGGTTAAAAACATAAACAAAGCAACCGCAAACCGCCGTTGGGCGAGGTGCTGGc                                                                                            <  1:959939/143‑1 (MQ=255)
                                                           aaTCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGc                                                                                            >  2:641101/1‑149 (MQ=255)
                                                              ccccTGCCGATCCCCTGTCAGGTAAATGATAACCC‑GAGTATGGC‑GCCCCGGAGGAATTGTGCTCACAGCCCGCGGGTGTTTGCGCTGTAAAACGATTAACAACACCAACAAAGCAAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGt                                                                                         <  1:958066/146‑1 (MQ=255)
                                                               cACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCg                                                                                       >  2:922808/1‑150 (MQ=255)
                                                                 cTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGAGCTGGCAGTCGATc                                                                                    >  2:1005134/1‑151 (MQ=255)
                                                                    ccaaTCCCCTGTCAGGTAAATGATGACCC‑GAGTAGGGCGGCCCGGGAGAAA‑TGTCCTCACAGCCCGCAGGTGATGGCGCTGTAAAACGAGTTACAACATCAACAAAGCGAACGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTg                                                                                  <  1:1058377/147‑1 (MQ=255)
                                                                    ccGATCCCCTGTCGGGTGAATGATGAACG‑GAGTAGGGCGGCCGTGGAGAAA‑TGTGCTCACAGCCCGCAGGTGTTTGCGTTTTAAAACGAGTTAAAACATCAACAAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATc                                                                                    <  1:552953/148‑1 (MQ=255)
                                                                      gATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGtt                                                                                >  2:519334/1‑150 (MQ=255)
                                                                      gATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGtt                                                                                >  1:301188/1‑150 (MQ=255)
                                                                         ccccTGTCGGGTAAATGATGAACC‑GAGTAGGGCGCCCCTGGAGAAA‑TGTCCTCACCCCCCGGGGGTGTTTGCGTTGTAAAAGGGGTTACAACAACAAGAAAGCAACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGAt                                                                                     <  2:352507/140‑1 (MQ=255)
                                                                         cACCTGCAAGGTAAATGA‑GAACCGGAGTAGGGGGGCCCGGGAAAAA‑TGTGCCCACACCCCGCGGGTGATGGCGTTTTAAAACGATTAACAACATCAACAAAGCAAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTaaa                                                                             <  1:810771/150‑1 (MQ=255)
                                                                          cccTGCCAGGTAAAGGAGAACCC‑GAGTTGGGCGCCCTTGGACCAA‑TGTCCCCACCGCCCGCGGGTTATGGCGTTGTAAACGGAGTTACAACAACAACAAGGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTaa                                                                              <  2:73131/147‑1 (MQ=255)
                                                                          aCCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAAtc                                                                           >  1:734257/1‑151 (MQ=255)
                                                                           ccggtgaaggaaaTGGTGACCG‑GAGTAGGGGGGCCCGGGAGCAA‑GGGCCTCGCGCCCGGCAGGTGATTGCTTTTTTAAGCGAGTTAAAACAACAAGAAGGGGAGCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGtt                                                                                <  2:381012/134‑1 (MQ=255)
                                                                           ccTGTCAGGAAAATGATGACCC‑GAGTATGGGGCCCCTGGGGCAA‑TGCGCTCCCGCCCCGCGGGTGATTGCGCTTTAAAACGATTTACAACATCAACTAGGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGtt                                                                                <  1:298976/145‑1 (MQ=255)
                                                                              gtcgggaaattaatgaCCG‑GAGTAGGGCGCCCCTGGAGAAT‑TGCCCCCACGGCCCGCAGGTTTTGGCGCTGTAAAACGAGTAACAACACCAACAAGGCGACCCCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAAtctc                                                                         <  1:438643/137‑1 (MQ=255)
                                                                               tCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAAc                                                                      >  1:683723/1‑151 (MQ=255)
                                                                                  ggTGAATGATGAACG‑GAGTGGGGCGGCCCGGGAGCAA‑TGTCCCCACGCCCCGCAGGTGATGGCGTTGTAAAACGAGTTACAACACCAACTAAGCAACCCCCACCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGcc                                                                   <  2:300530/151‑1 (MQ=255)
                                                                                  ggTGAATGATGAACC‑GAGTAGGGCGCCCCTGGAGAAA‑TGTGCTCCCGGCCCGCAGGTGTTTGCGCTGTAAAACGATTTACAACATCAACAAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGc                                                                    <  1:590532/150‑1 (MQ=255)
                                                                                  ggTGAATGATAACCC‑AAGTAGGGCGGCCCGGGAGAAT‑TGTGCCCACCCCCCGGGGGTGGTTGCCCGTTAAAAGGATTTACAACAACAACTAGGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCaa                                                                       <  1:792932/147‑1 (MQ=255)
                                                                                           tgaACG‑GAGTGGGGCGCCCTGGGAAAAT‑TGTGCTCACAGCCGGAGGGTTATGGCGTTTTAAAACGTTTTAAAACATCAACTAAGCGAACCCCATCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGAt                                                             <  1:582898/148‑1 (MQ=255)
                                                                                               cc‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACt                                                      >  2:911485/1‑151 (MQ=255)
                                                                                                cggAGTAGGGCGGCCCGGGAGAAA‑TGTCCCCACGGCCCGCGGGTGTTTGCGTTGTAAAACGAGTTACAACATAAACAAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATc                                                            <  1:615982/143‑1 (MQ=255)
                                                                                                      aTGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATc                                                 >  2:628276/1‑150 (MQ=255)
                                                                                                        ggcggcCCGGGAGCAA‑TGTGCCCACACCCCGCAGGTGATGGCGTTGTAAAACGAGTAAAAACAACAACAAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCa                                                <  1:11071/149‑1 (MQ=255)
                                                                                                        ggcgccCCGGGAGAAA‑TGTCCCCCCACCCCGGAGGTGTTTGCGTTTTAAAAGGAGTTACACCATCAACAACGCGACCCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAAt                                              <  2:880628/151‑1 (MQ=255)
                                                                                                               cTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTtc                                       >  2:825852/1‑151 (MQ=255)
                                                                                                                 ggAGAAA‑TGTGCCCACAGCCCGCAGGTGTTGGCGTTGTAAAACGAGTAACAACATAAACAAAGCGACCGCCAGCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTtc                                       <  1:104287/149‑1 (MQ=255)
                                                                                                                          gtgCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaagg                             >  2:830309/1‑151 (MQ=255)
                                                                                                                          gtccccacaCCCCGGAGTTGTTTCCGTTTTAAACCGATTTACAACACCAACAAACCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaa                               <  1:1054378/144‑1 (MQ=255)
                                                                                                                           tgCTCCCACCCCGCGGGTGATTGCGTTGTAAAACGAGTTACAACATCAACTAAGCAAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCg                                    <  2:748486/143‑1 (MQ=255)
                                                                                                                             cccacaGCCCGCGGGTGATTGCGTTTTAAAACGAGTTACAACATCAACAAAGCAAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaagggg                           <  1:367903/148‑1 (MQ=255)
                                                                                                                              tCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaagggga                          >  2:676132/1‑150 (MQ=255)
                                                                                                                                    cccGGAGGTATTTGCGCTATAAACGGTGTTACAACACCAAAAAGGCCAACGCCACCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaaggggaa                         <  1:900418/145‑1 (MQ=255)
                                                                                                                                      ccggggTGATCGCGCTTTAAAACGAGTAACACCATCAACAACGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAg                  <  2:1078664/146‑1 (MQ=255)
                                                                                                                                          ggTGATGGCGTTGTAAAACGATTAACAACATAAAATAAGCGACCCCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTg                <  2:142330/148‑1 (MQ=255)
                                                                                                                                           gTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTgcgc             >  2:22685/1‑150 (MQ=255)
                                                                                                                                              aTTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTc         >  1:639898/1‑151 (MQ=255)
                                                                                                                                                   gcTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTcagc      >  1:710908/1‑149 (MQ=255)
                                                                                                                                                     tGTAAACCGAGTAAAAACACCAACAAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTcagcag    <  1:667936/149‑1 (MQ=255)
                                                                                                                                                     tGTAAAACGAGTTACAACATCAACTAAGCGACCGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTca        <  1:604308/145‑1 (MQ=255)
                                                                                                                                                      gTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGc  >  2:267330/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
AAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGGCGGCCCTGGAGCAA‑TGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGC  >  NC_000913/1036756‑1037053

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