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

TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTC  >  NC_000913/1764637‑1764929
                                                                                                                                                     |                                                                                                                                                
tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                 >  2:783253/1‑150 (MQ=255)
tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                 >  2:787755/1‑150 (MQ=255)
 aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                >  2:640352/1‑150 (MQ=255)
 aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                <  1:1011224/150‑1 (MQ=255)
   gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                                                                                                             <  2:303557/151‑1 (MQ=255)
     gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTa                                                                                                                                           <  2:377094/151‑1 (MQ=255)
                   tGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTa                                                                                                                              >  2:97843/1‑150 (MQ=255)
                      aGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAacac                                                                                                                          >  2:730841/1‑151 (MQ=255)
                          ttCAGAAATAAGAAAACCCTTCAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAAc                                                                                                                      >  1:1226623/1‑151 (MQ=255)
                            aaaaaaaaaGAAACCCCTTAAGTCTGTGCGACTCAGGCTAAGGGTTTTCTCCCCCACCCGGCGCTTATCGCCGGCACTCTGGGTGGTTTAGCTCTTGAAGGGGCGAAAGGAAAAATCTCAAAAGGCTAACCCGCCGTTTTAACACAAACTg                                                                                                                    <  1:1086473/148‑1 (MQ=255)
                             aGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                   >  2:90289/1‑151 (MQ=255)
                             aGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                   >  2:874407/1‑151 (MQ=255)
                             aGAAAAAAAAAACCCCTTAATTCTGTGCGACACAGGTTTAGGGTTTTCTCCCCCACCCGGCGCTTATCCCCGCCACTCTCAGGGGTTTAGCTTTTGAGGGGGCAAAAAAAAAAATCTCAAAAACCTACCCCGCCGTTTTAACACAAACTg                                                                                                                    <  1:946269/150‑1 (MQ=255)
                                  aaaGAAACCCCTAAAGTCTGTGCGACACGGGTTTAGGGTTTTCTCCCCCACCCGGCTTTTATCCCCGCCACTCTCGGTGGTTTAGCTTTTGAAGGGGCAATAAAAATAATCTAATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTa                                                                                                              <  1:1106622/150‑1 (MQ=255)
                                   aaaaaacCCCTTAAGTCTGTGCACCCCAGTTTTAAGGTTTTCTCCCCCACCCGCCCCTTCTCCCCGCCCCTCCCAGGGTTTTAGCTTTTGAGGGGGGGATAAAAAAAATCTCAAAAGCCTACCCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                             <  1:154805/148‑1 (MQ=255)
                                   aaaaaacCCCTAAAGTCTGGGCGACACAGGCTAAGGGTTTCCTCCCCCACCCGGCGCTTATCTCCGCCCCTCACAGTGGCTTCCCTTTTAAGGGGGCAAAAAGAATAATCAAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                             <  1:624704/148‑1 (MQ=255)
                                   aaaaaaaCCTTTAAGTTTGTGCGACACGGGCTTAAGGTTTTCTCCCCCACCCGGCGTTTATCCCCGGCCCCCTCAGTGGCTTCCCTCTGGAAGGGGCGATAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTa                                                                                                              <  2:212601/147‑1 (MQ=255)
                                   aaGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                               >  1:924694/1‑149 (MQ=255)
                                          cccTTAAGTCTGTGCGACACAGGCTTAAGGTTTTCTCCCCCATCCGGCCCTTATCTCCGCCCCTCTCAGTGGCTTAGCTTTTGAAGGGGCGATAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGtatta                                                                                                        <  1:903113/149‑1 (MQ=255)
                                          cccTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTattatt                                                                                                      >  1:880297/1‑151 (MQ=255)
                                            cTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                    >  1:182773/1‑151 (MQ=255)
                                               aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                  >  1:1228678/1‑150 (MQ=255)
                                               aaGTCTGGGCGACACAGTCTTAAGGTTTTCCCCCCCACCCGGCCCTTATCCCCGGCACTCTCAGTGGCTTAGCTTTTGAAGGGGAAAAAAAAATAATCTAATAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGa                                                                                                 <  1:899077/151‑1 (MQ=255)
                                                  tctggccaccacaGGCTTAAGGGTTTCTCCCCCATCCGGCCCTTATCCCCGCCCCTCTCAGTGGTTTCGCTTTTGAAGGGGCAAAAAAAATAATCTTAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                  <  1:540327/138‑1 (MQ=255)
                                                     gtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAat                                                                                            >  1:22327/1‑150 (MQ=255)
                                                         gACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                         >  2:1210319/1‑149 (MQ=255)
                                                           cacaGGTTTAAGGTTTTCTCCCCCACCCGGCGTTTATCCCCGCCACTCTAAGTGTTTTAGCTTTTGAAGGGGCAATAAAAATAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACaa                                                                                             <  1:1083355/143‑1 (MQ=255)
                                                            acaGGTTTAGGGGTTCCTCCCCCACCCGGCGTTTATCCCCGGCACTCTCAGTGGTTTAGCTCTTGAGGGGGCGAAAAAAAAAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                     <  2:605377/150‑1 (MQ=255)
                                                               ggtttAGGGTTTTCTCCCCCACCCGGCGTTTATCTCCGCCCCTCTCGGGGGTTTAGCTTTTGAGGGGGCAATAAAAATAATCTAAAAAATCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                     <  1:267545/144‑1 (MQ=255)
                                                                      ggttttCTACCCCACCCGGCGTTTATCTCCGGCCCTCTCAGTGGCTTAGCTTTGGAAGGGGCAAAAAAAATAATTAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAg                                                                           <  1:172366/147‑1 (MQ=255)
                                                                       ggTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGc                                                                          >  1:296672/1‑150 (MQ=255)
                                                                          tttcccccccACCGGGCGCTTACCTCCGGCACTCTCAGTGGTTTAGCTTTTGAGGGGGCAATAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCa                                                                         <  2:660019/143‑1 (MQ=255)
                                                                          ttCTCCCCCACCCGGCCTTTACCCCCGGCCCTCTCAGTGGTTTAGCTTTTGAAGGGGCGAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt                                                                       <  2:727623/150‑1 (MQ=255)
                                                                             cTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt                                                                    >  2:141639/1‑151 (MQ=255)
                                                                                 cccATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa                                                                   >  1:48960/1‑148 (MQ=255)
                                                                                 cccACCCGGCTCTTATCCCCGGCCCTCTCAGTGGTTTAGCTTTTGAAGGGGCAATAAAAATAATCTAAAAAACCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga                                                                <  1:853490/151‑1 (MQ=255)
                                                                                 cccACCCGGCGTTTACCCCCGGCCCTCTCAGTGGCTTACCTCTTGAGGGGGCAAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga                                                                <  2:217526/151‑1 (MQ=255)
                                                                                  ccccccGGCGCTTATCCCCGCCCCTCTCAGTGGTTTAGCTTTTGAGGGGGAAAAAAAAATAATCTCAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaa                                                               <  1:616543/147‑1 (MQ=255)
                                                                                  ccACCCGGCTCTTATCCCCGCCCCTCTCGGTGGCTTAGCTTTTGAAGGGGCAATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaa                                                               <  2:395645/151‑1 (MQ=255)
                                                                                       cGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                           >  1:529673/1‑150 (MQ=255)
                                                                                        ggCGCTTATCTCCGGCCCTCTCAGTGGTTTAGCTTTTGAGGGGGCAATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCa                                                          <  2:998745/150‑1 (MQ=255)
                                                                                         gcgcTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATc                                                        >  2:375379/1‑151 (MQ=255)
                                                                                                 ccccgcccctctCAGTGGCTTAGCTCTTGAAGGGGCAAAAAAAATAATCAAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAAc                                                    <  2:46427/139‑1 (MQ=255)
                                                                                                 ccccGGCCCCCTCAGTGGTTTAGCTCTTAAAGGGGCAATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                  <  2:856911/147‑1 (MQ=255)
                                                                                                     gcccctctCAGTGGTTTAGCTCTTGAAGGGGCAAAAAAAATAATCTCAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACg                                             <  1:1146606/146‑1 (MQ=255)
                                                                                                       cACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACgg                                            >  1:995763/1‑149 (MQ=255)
                                                                                                        aCTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAgg                                         >  2:222505/1‑151 (MQ=255)
                                                                                                             cAGTGGTTTAGCTCTTGAAGGGGCGAAAAAAATAATCTAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTa                                      <  2:138356/149‑1 (MQ=255)
                                                                                                             cAGTGGCTTAGCTCTTGAAGGGGCAAAAAAAATAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACg                                             <  1:592386/142‑1 (MQ=255)
                                                                                                                 ggtttAGCTTTTGAGGGGGGAATAAAAAAAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCa                                 <  2:944507/147‑1 (MQ=255)
                                                                                                                 ggttaaGCTCTTGAAGGGCAAAAAAAAATAATCTAATAAAGTCACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                  <  1:522812/144‑1 (MQ=255)
                                                                                                                 ggCTTAGTTCTGGAGGGGGAAATAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAg                                <  1:456271/151‑1 (MQ=255)
                                                                                                                    ttAGCTTTTGAGGGGGCAAAAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGt                              <  1:1185728/150‑1 (MQ=255)
                                                                                                                     tAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                            >  2:237034/1‑151 (MQ=255)
                                                                                                                      aGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                            >  2:16756/1‑150 (MQ=255)
                                                                                                                      aGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                            >  2:249984/1‑150 (MQ=255)
                                                                                                                        ctctTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAtt                         >  2:676590/1‑151 (MQ=255)
                                                                                                                         tctTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAttt                        >  2:1228743/1‑151 (MQ=255)
                                                                                                                                ggggaaaaaaaaataatCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAAc                 <  1:1119321/140‑1 (MQ=255)
                                                                                                                                      aaaaGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCg            <  2:1064915/148‑1 (MQ=255)
                                                                                                                                              aatCTAAAAAACCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAg    <  1:770756/150‑1 (MQ=255)
                                                                                                                                                tctcATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTc  >  2:49644/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTC  >  NC_000913/1764637‑1764929

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