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

TTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATC  >  NC_000913/1764636‑1764935
                                                                                                                                                         |                                                                                                                                                      
ttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt                                                                                                                                                        >  1:1152402/1‑150 (MQ=255)
 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                       <  1:1216800/150‑1 (MQ=255)
     agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                  <  2:1146330/151‑1 (MQ=255)
      gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTa                                                                                                                                                 <  1:824591/151‑1 (MQ=255)
        tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTaa                                                                                                                                                <  1:859780/150‑1 (MQ=255)
        tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                                                                                                                               >  1:1162501/1‑151 (MQ=255)
         aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAcc                                                                                                                                              <  2:245799/151‑1 (MQ=255)
             aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                                          >  2:794029/1‑151 (MQ=255)
             aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAACCCATCAGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                                          >  2:535680/1‑151 (MQ=255)
                aaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtt                                                                                                                                       >  1:1082923/1‑151 (MQ=255)
                   gTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtttt                                                                                                                                     >  1:950460/1‑150 (MQ=255)
                    tGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTa                                                                                                                                    >  1:918883/1‑150 (MQ=255)
                    tGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTa                                                                                                                                    >  1:325702/1‑150 (MQ=255)
                         cGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAAcaca                                                                                                                               >  1:1157858/1‑150 (MQ=255)
                             cAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTg                                                                                                                          >  1:122461/1‑151 (MQ=255)
                              aaaaaaaaaaaaCCCCTAAAGTCTGTGCGACACAGGTTAAAGGTTTTCTCCCCCACCCGGCTTTTATCCCCGCCCCTCTCGGGGGCTTAGCTTTTGAAGGGGCAATAAGAAAAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGc                                                                                                                         <  2:625683/149‑1 (MQ=255)
                              aGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                         >  2:358055/1‑151 (MQ=255)
                                  aaaaaaaacCCCTTAAGTCTGTCCGACGGGGGGTAAGGGCTTTCTCCCCAACCGGGACCTAACCTCCGCGTCCCTCGGTGTCTACGCTATTGAGGGGGAAATAAAAATAAACTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                     <  1:864232/146‑1 (MQ=255)
                                  aTAAAAAACCCCTAAAGTCTGTGCGACACGGGCTAAGGGTTTTCTCCCCCACCCGGCGCTTATCTCCGGCCCTCGCAGTGGCTTAGCTTTTGAAGGGGCGAAAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                     <  1:58537/151‑1 (MQ=255)
                                   tAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                     >  1:1123841/1‑150 (MQ=255)
                                    aaaaaaaCCCTTAATTTTGGGCCACCCAGGCTTAAGGGTTTCTCCCCCACCCGGCCTTTATCCCCGGCCCTCCCAGTGTCTTTCCTCTTAGAGGGGCAATAAGAATAACCTCAAAAAACTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                     <  2:909651/146‑1 (MQ=255)
                                    aaGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGGGAt                                                                                                                      >  2:352924/1‑148 (MQ=255)
                                       aaaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtat                                                                                                                >  1:80615/1‑151 (MQ=255)
                                        aaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatt                                                                                                               >  1:1212372/1‑151 (MQ=255)
                                        aaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGCTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatt                                                                                                               >  2:349225/1‑151 (MQ=255)
                                         accctttAAGTCTGTGCCACACAGGTTTTAGGGTTTCTCCCCCACCCGGCGTTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAGGGGGCGATAAGAATAATCTAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatta                                                                                                              <  1:856143/146‑1 (MQ=255)
                                            ccTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                                                                           >  1:1239146/1‑151 (MQ=255)
                                               tAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                        >  1:430544/1‑151 (MQ=255)
                                                aggTCTGTGCGAAACAGGTTAAGGGTTTTCTACCCCACCCGGCGCTTATCCCCGCCACTCTCAGTGGCTTACCTTTGGAGGGGGCAATAAAAATAATTTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt                                                                                                         <  2:1143832/147‑1 (MQ=255)
                                                 aGTCTGTGCGACACAGGCTAAAGGGTTTCCCCCCCACCCGGCTCTTATCCCCGGCACTCTCAGGGGTTTAGCTTTTGAGGGGGAGATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                      <  1:920709/151‑1 (MQ=255)
                                                   tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCACTCGCAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAATCTGCGCTTAGTATTATTTTTGAAc                                                                                                     >  2:4112/1‑150 (MQ=255)
                                                    cTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                                                                    >  1:97251/1‑150 (MQ=255)
                                                      gtgCGACACGGGCTTAGGGGTTTCTCCCCCACCCGGCTTTTATCTCCGGCACTCTCAGTGGTTTAGCTTTTGAAGGGGAAAAAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt                                                                                                         <  1:874658/143‑1 (MQ=255)
                                                           acacAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAgg                                                                                            >  2:201705/1‑151 (MQ=255)
                                                              caGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg                                                                                         >  1:187683/1‑151 (MQ=255)
                                                                   ttAAGGGTTTCCCCCCCATCCGGCGTTTATCCCCGCCACCCTCAGTGGCTTAGCTTTTGAGGGGGCAATAAGAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                                        <  1:406324/147‑1 (MQ=255)
                                                                      aagcttttcacccccAGCCGCCGCTTATCCCCGCCCCCCTCTGTGGTTTACTTTTGGAAGGGGCAAAAAAAATAATTTTAAAATTCCCCCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGa                                                                                     <  2:1027169/137‑1 (MQ=255)
                                                                        ggTTTCCTCCCCCTCCGGGCGCTTATCCCCGGCACTCTAAGTGTTTTAGCTTTTGAAGGGGCAATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGc                                                                                <  2:549562/151‑1 (MQ=255)
                                                                                 ccccACCCGGCGTTTATCTCCGGCCCTCGCAGTGGTTTAGCTCTTGAGGGGGCGATAAGAATAATCTCAAAAAACTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt                                                                          <  2:933945/148‑1 (MQ=255)
                                                                                  cccACCCGGCTCTTATCCCCGGCACTCTCAGTGGTTTAGCTCTGAAGGGGCGGATAAAAATACTCTAAAAAAGTTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTa                                                                            <  1:677839/145‑1 (MQ=255)
                                                                                  cccACCCGGCGCTCATCCCCGGCCCTCCCGGTGGCTTTGCTCTGGAAGGGGCAAAAAGAATAATCTCATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga                                                                      <  2:288525/151‑1 (MQ=255)
                                                                                     aTCCGCCGCTTATCCCCGCCACTCTCAGTGGTTTGGCTTTTGAAGGGGCAAAAAAAAAAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaag                                                                    <  2:1221472/150‑1 (MQ=255)
                                                                                       ccGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                                 >  2:266104/1‑151 (MQ=255)
                                                                                       ccGGCGCTTATCCCCGGCCCTCTCGGTGGTTTAGCTTTTGAAGGGGCAAAAAAAAAAACCAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAgg                                                                   <  1:299628/149‑1 (MQ=255)
                                                                                           gcgcTTATCCCCGCCACCCTCGGTGGTTTACCTCTGGAAGGGGCAATAAAAATAATCCAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                                 <  1:866686/148‑1 (MQ=255)
                                                                                               ttATCTCCGGCCCTCTAAGTGGTTTACTTCTTAAAGGGGCGAAAAGAAAAATCAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGaa                                                           <  2:830652/150‑1 (MQ=255)
                                                                                                tATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAAcc                                                         >  1:200536/1‑151 (MQ=255)
                                                                                                     ccGCCACTCTCAGTGTTTTAGCTTTGAAGGGGGAAATAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                        <  2:459800/147‑1 (MQ=255)
                                                                                                     ccGACACTCTGAGTGGTTTAGCTCTTGAGGGGGCAATAAGAATAATCAAAAAAATTAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                        <  1:238292/147‑1 (MQ=255)
                                                                                                               cAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAAt                                          >  2:56081/1‑151 (MQ=255)
                                                                                                                  tGGTTTAGCTTTTGAGGGGGCTATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                        <  1:658807/150‑1 (MQ=255)
                                                                                                                   gttttATCTCTTGAAGGGGAGAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCa                                       <  1:1027628/147‑1 (MQ=255)
                                                                                                                   ggCTTAGCTCTTGAGGGGGCGAAAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                        <  1:602857/149‑1 (MQ=255)
                                                                                                                   ggCTTAGATCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                        >  2:141157/1‑149 (MQ=255)
                                                                                                                         gCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAt                                >  2:280251/1‑151 (MQ=255)
                                                                                                                             ttGAAGGGGCGATAAAAAAAATCTAATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAttt                              <  1:401812/149‑1 (MQ=255)
                                                                                                                              tGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGc                           >  2:468625/1‑151 (MQ=255)
                                                                                                                                 agggggCAATAAAAATAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                                  <  1:702773/140‑1 (MQ=255)
                                                                                                                                    gggCGATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACt                      <  1:500215/151‑1 (MQ=255)
                                                                                                                                          aaaaaaataatCTCAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGc                   <  1:963513/144‑1 (MQ=255)
                                                                                                                                            aaaaataatCTCAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTg                    <  2:1008736/143‑1 (MQ=255)
                                                                                                                                            aaaaataatCTAAAAACGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGc                   <  2:1117894/144‑1 (MQ=255)
                                                                                                                                            aaaaataacCTAAAAAATCTAGCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTg                    <  2:557001/143‑1 (MQ=255)
                                                                                                                                             aGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCAt              >  1:823455/1‑151 (MQ=255)
                                                                                                                                               aatattCACATAAAGTTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCAtt             <  1:1103043/150‑1 (MQ=255)
                                                                                                                                               aataatCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTa            <  2:1152381/151‑1 (MQ=255)
                                                                                                                                                         tAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATc  >  2:1149004/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                      
TTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATC  >  NC_000913/1764636‑1764935

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