Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 823,777 (A)5→6 intergenic (‑57/‑76) ybhP ← / → ybhQ hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913823,7721.A99.3% 492.4 / ‑3.6 136intergenic (‑52/‑81)ybhP/ybhQhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (55/80);  total (55/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

Only 100 of 157 total aligned reads displayed.

GTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGT  >  NC_000913/823680‑823859
                                                                                              |                                                                                       
gTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcac                                                                                          <  1:82575/93‑1 (MQ=255)
 tGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcaca                                                                                        <  2:3997398/93‑1 (MQ=255)
 tGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcaca                                                                                        <  2:3442784/93‑1 (MQ=255)
   aGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaa                                                                                      >  2:906088/1‑93 (MQ=255)
   aGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaa                                                                                      >  1:52129/1‑93 (MQ=255)
    gCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaa                                                                                     >  2:231823/1‑93 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaa                                                                                   >  1:4555089/1‑92 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaa                                                                                   >  1:4613528/1‑92 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaa                                                                                   <  2:2155383/93‑2 (MQ=255)
        cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaagg                                                                                 <  2:2545586/93‑4 (MQ=255)
        cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaagg                                                                                 <  1:3081931/93‑4 (MQ=255)
        cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaagg                                                                                 <  1:4203806/93‑4 (MQ=255)
        cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaagg                                                                                 <  2:2056541/93‑4 (MQ=255)
         ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaaggt                                                                                <  2:224226/93‑5 (MQ=255)
           gAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGc                                                                               >  2:4526682/1‑93 (MQ=255)
           gAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGc                                                                               >  2:123688/1‑93 (MQ=255)
            aaCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCg                                                                              >  2:1387184/1‑93 (MQ=255)
             aCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCgg                                                                             >  1:930396/1‑93 (MQ=255)
               gAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTc                                                                           <  2:1944797/93‑1 (MQ=255)
               gAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTc                                                                           <  2:158672/93‑1 (MQ=255)
                aaaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCt                                                                          <  2:4456688/93‑1 (MQ=255)
                aaaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCt                                                                          <  2:2810501/93‑1 (MQ=255)
                 aaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCtt                                                                         <  2:636954/93‑1 (MQ=255)
                  aaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTa                                                                        <  2:3670609/93‑1 (MQ=255)
                   aTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTAc                                                                       >  1:883819/1‑93 (MQ=255)
                   aTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTAc                                                                       <  1:1976364/93‑1 (MQ=255)
                    ttgttgTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACg                                                                      >  1:3797641/1‑93 (MQ=255)
                       ttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGaa                                                                   >  1:759220/1‑93 (MQ=255)
                        tgtgTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAt                                                                  >  2:2714841/1‑93 (MQ=255)
                        tgtgTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAt                                                                  <  1:588316/93‑1 (MQ=255)
                             ttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTcc                                                             <  2:319082/93‑1 (MQ=255)
                             ttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTcc                                                             >  1:194984/1‑93 (MQ=255)
                              tGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCg                                                            <  1:4675519/93‑1 (MQ=255)
                              tGATCGGGCACAATTTTCCTGTCTTTGCCTCTAATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCg                                                            >  2:3249996/1‑93 (MQ=255)
                                  cGGGCATTATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                        <  2:4335831/93‑1 (MQ=255)
                                   gggCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAg                                                       >  1:2647022/1‑93 (MQ=255)
                                     gCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGtt                                                     <  1:2313170/93‑1 (MQ=255)
                                      cATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTc                                                    <  2:1774954/93‑1 (MQ=255)
                                      cATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTc                                                    <  1:2209736/93‑1 (MQ=255)
                                      cATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTc                                                    <  1:4082686/93‑1 (MQ=255)
                                      cATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTc                                                    >  1:2403552/1‑93 (MQ=255)
                                       aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                   <  1:2270686/93‑1 (MQ=255)
                                       aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                   >  2:3616132/1‑93 (MQ=255)
                                        tAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCgg                                                  <  2:2680191/93‑1 (MQ=255)
                                        tAATTTGCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAGAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCgg                                                  <  2:3505041/93‑1 (MQ=255)
                                          aTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTa                                                >  1:3002529/1‑93 (MQ=255)
                                           ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               <  1:1060762/93‑1 (MQ=255)
                                           ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               >  1:1283440/1‑93 (MQ=255)
                                             ttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActc                                             <  1:311047/93‑1 (MQ=255)
                                               ccTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTg                                           <  2:4201702/93‑1 (MQ=255)
                                               ccTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTg                                           <  1:3339061/93‑1 (MQ=255)
                                                 tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                         >  1:1737589/1‑93 (MQ=255)
                                                 tGTCTTTGCCTCTTATCTCATTGAAATAGCGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGCAATTTTCCGTAAAGTTCGGTACTCTGag                                         >  1:1597115/1‑93 (MQ=255)
                                                  gTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagt                                        >  1:3467274/1‑93 (MQ=255)
                                                  gTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGGGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagt                                        >  1:1386349/1‑93 (MQ=255)
                                                     tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaag                                     >  2:2903640/1‑93 (MQ=255)
                                                     tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaag                                     <  2:1906503/93‑1 (MQ=255)
                                                       tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACTGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   <  2:3112249/93‑1 (MQ=255)
                                                       tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   <  2:1416637/93‑1 (MQ=255)
                                                       tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   <  2:1233662/93‑1 (MQ=255)
                                                          ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTagag                                >  1:1121609/1‑93 (MQ=255)
                                                          ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTagag                                >  1:60264/1‑93 (MQ=255)
                                                           tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAgaga                               <  2:2439002/93‑1 (MQ=255)
                                                              tATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATaa                            >  2:2805373/1‑93 (MQ=255)
                                                              tATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATaa                            <  1:649938/93‑1 (MQ=255)
                                                               aTCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATaaa                           <  1:2419381/93‑1 (MQ=255)
                                                                tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAt                          <  2:782566/93‑1 (MQ=255)
                                                                 ctcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAtt                         <  2:2194513/93‑1 (MQ=255)
                                                                  tcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttc                        <  2:4071179/93‑1 (MQ=255)
                                                                   cATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttct                       <  1:3598661/93‑1 (MQ=255)
                                                                   cATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttct                       >  2:4556777/1‑93 (MQ=255)
                                                                    aTTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                      >  1:1314274/1‑93 (MQ=255)
                                                                    aTTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                      <  1:4663439/93‑1 (MQ=255)
                                                                    aTTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                      <  1:1853517/93‑1 (MQ=255)
                                                                    aTTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                      >  1:2201545/1‑93 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     >  2:2729380/1‑93 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     <  2:3002529/93‑1 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     >  1:1422958/1‑93 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     >  1:2137545/1‑93 (MQ=255)
                                                                      tGACATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                    <  1:1678662/93‑1 (MQ=255)
                                                                      tGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                    <  2:1787546/93‑1 (MQ=255)
                                                                      tGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                    <  1:1380386/93‑1 (MQ=255)
                                                                      tGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                    <  1:2339272/93‑1 (MQ=255)
                                                                       gAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                   <  2:2099194/93‑1 (MQ=255)
                                                                       gAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                   >  1:1355975/1‑93 (MQ=255)
                                                                       gAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                   <  2:607822/93‑1 (MQ=255)
                                                                         aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  1:1026574/93‑1 (MQ=255)
                                                                         aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  2:407117/93‑1 (MQ=255)
                                                                         aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  1:3374322/93‑1 (MQ=255)
                                                                         aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  1:495372/93‑1 (MQ=255)
                                                                              tgtAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGaa            >  1:1276499/1‑93 (MQ=255)
                                                                                tAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGc          >  1:2683789/1‑93 (MQ=255)
                                                                                 aGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGcc         <  1:669650/93‑1 (MQ=255)
                                                                                  gTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa        >  2:3900371/1‑93 (MQ=255)
                                                                                   tCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCAt       <  1:1102059/93‑1 (MQ=255)
                                                                                   tCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCAt       >  2:2198356/1‑93 (MQ=255)
                                                                                    cGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATg      >  2:632021/1‑93 (MQ=255)
                                                                                     ggCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa     <  1:2083165/93‑1 (MQ=255)
                                                                                      gCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGaa    <  2:3425051/93‑1 (MQ=255)
                                                                                       cGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAg   <  1:161197/93‑1 (MQ=255)
                                                                                              |                                                                                       
GTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGT  >  NC_000913/823680‑823859

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