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

TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTC  >  NC_000913/1764637‑1764929
                                                                                                                                                     |                                                                                                                                                
tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                 >  2:783253/1‑150 (MQ=255)
tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                 >  2:787755/1‑150 (MQ=255)
 aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                >  2:640352/1‑150 (MQ=255)
 aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                <  1:1011224/150‑1 (MQ=255)
   gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                                                                                                             <  2:303557/151‑1 (MQ=255)
     gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTa                                                                                                                                           <  2:377094/151‑1 (MQ=255)
                   tGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTa                                                                                                                              >  2:97843/1‑150 (MQ=255)
                      aGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAacac                                                                                                                          >  2:730841/1‑151 (MQ=255)
                          ttCAGAAATAAGAAAACCCTTCAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAAc                                                                                                                      >  1:1226623/1‑151 (MQ=255)
                            aaaaaaaaaGAAACCCCTTAAGTCTGTGCGACTCAGGCTAAGGGTTTT‑CTCCCCCACCCGGCGCTTATCGCCGGCACTCTGGGTGGTTTAGCTCTTGAAGGGGCGAAAGGAAAAATCTCAAAAGGCTAACCCGCCGTTTTAACACAAACTg                                                                                                                    <  1:1086473/148‑1 (MQ=255)
                             aGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                   >  2:90289/1‑151 (MQ=255)
                             aGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                   >  2:874407/1‑151 (MQ=255)
                             aGAAAAAAAAAACCCCTTAATTCTGTGCGACACAGGTTTAGGGTTTT‑CTCCCCCACCCGGCGCTTATCCCCGCCACTCTCAGGGGTTTAGCTTTTGAGGGGGCAAAAAAAAAAATCTCAAAAACCTACCCCGCCGTTTTAACACAAACTg                                                                                                                    <  1:946269/150‑1 (MQ=255)
                                  aaaGAAACCCCTAAAGTCTGTGCGACACGGGTTTAGGGTTTT‑CTCCCCCACCCGGCTTTTATCCCCGCCACTCTCGGTGGTTTAGCTTTTGAAGGGGCAATAAAAATAATCTAATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTa                                                                                                              <  1:1106622/150‑1 (MQ=255)
                                   aaaaaacCCCTTAAGTCTGTGCACCCCAGTTTTAAGGTTTT‑CTCCCCCACCCGCCCCTTCTCCCCGCCCCTCCCAGGGTTTTAGCTTTTGAGGGGGGGATAAAAAAAATCTCAAAAGCCTACCCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                             <  1:154805/148‑1 (MQ=255)
                                   aaaaaacCCCTAAAGTCTGGGCGACACAGGC‑TAAGGGTTTCCTCCCCCACCCGGCGCTTATCTCCGCCCCTCACAGTGGCTTCCCTTTTAAGGGGGCAAAAAGAATAATCAAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                             <  1:624704/148‑1 (MQ=255)
                                   aaaaaaaCCTTTAAGTTTGTGCGACACGGGCTTAAGGTTTT‑CTCCCCCACCCGGCGTTTATCCCCGGCCCCCTCAGTGGCTTCCCTCTGGAAGGGGCGATAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTa                                                                                                              <  2:212601/147‑1 (MQ=255)
                                   aaGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                               >  1:924694/1‑149 (MQ=255)
                                          cccTTAAGTCTGTGCGACACAGGCTTAAGGTTTT‑CTCCCCCATCCGGCCCTTATCTCCGCCCCTCTCAGTGGCTTAGCTTTTGAAGGGGCGATAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGtatta                                                                                                        <  1:903113/149‑1 (MQ=255)
                                          cccTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTattatt                                                                                                      >  1:880297/1‑151 (MQ=255)
                                            cTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                    >  1:182773/1‑151 (MQ=255)
                                               aaGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                  >  1:1228678/1‑150 (MQ=255)
                                               aaGTCTGGGCGACACAGTCTTAAGGTTTT‑CCCCCCCACCCGGCCCTTATCCCCGGCACTCTCAGTGGCTTAGCTTTTGAAGGGGAAAAAAAAATAATCTAATAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGa                                                                                                 <  1:899077/151‑1 (MQ=255)
                                                  tctggccaccacaGGCTTAAGGGTTT‑CTCCCCCATCCGGCCCTTATCCCCGCCCCTCTCAGTGGTTTCGCTTTTGAAGGGGCAAAAAAAATAATCTTAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                  <  1:540327/138‑1 (MQ=255)
                                                     gtgCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAat                                                                                            >  1:22327/1‑150 (MQ=255)
                                                         gACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                         >  2:1210319/1‑149 (MQ=255)
                                                           cacaGGTTTAAGGTTTT‑CTCCCCCACCCGGCGTTTATCCCCGCCACTCTAAGTGTTTTAGCTTTTGAAGGGGCAATAAAAATAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACaa                                                                                             <  1:1083355/143‑1 (MQ=255)
                                                            acaGGTTTAGGGGTTC‑CTCCCCCACCCGGCGTTTATCCCCGGCACTCTCAGTGGTTTAGCTCTTGAGGGGGCGAAAAAAAAAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                     <  2:605377/150‑1 (MQ=255)
                                                               ggtttAGGGTTTT‑CTCCCCCACCCGGCGTTTATCTCCGCCCCTCTCGGGGGTTTAGCTTTTGAGGGGGCAATAAAAATAATCTAAAAAATCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                     <  1:267545/144‑1 (MQ=255)
                                                                      ggtttt‑CTACCCCACCCGGCGTTTATCTCCGGCCCTCTCAGTGGCTTAGCTTTGGAAGGGGCAAAAAAAATAATTAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAg                                                                           <  1:172366/147‑1 (MQ=255)
                                                                       ggTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGc                                                                          >  1:296672/1‑150 (MQ=255)
                                                                          tttcccccccACCGGGCGCTTACCTCCGGCACTCTCAGTGGTTTAGCTTTTGAGGGGGCAATAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCa                                                                         <  2:660019/143‑1 (MQ=255)
                                                                          tt‑CTCCCCCACCCGGCCTTTACCCCCGGCCCTCTCAGTGGTTTAGCTTTTGAAGGGGCGAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt                                                                       <  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)
                                                                                                                                                     |                                                                                                                                                
TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTC  >  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: ‑