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.4% 654.7 / ‑3.8 174intergenic (‑52/‑81)ybhP/ybhQhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (87/86);  total (88/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 194 total aligned reads displayed.

TGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGT  >  NC_000913/823681‑823859
                                                                                            |                                                                                        
 gAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaa                                                                                        >  2:5816337/1‑93 (MQ=255)
     aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaa                                                                                   >  1:5866454/1‑92 (MQ=255)
     aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaa                                                                                   <  2:1489649/93‑2 (MQ=255)
      ccTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaag                                                                                  >  2:3333561/1‑91 (MQ=255)
        ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaaggt                                                                                <  2:4043330/93‑5 (MQ=255)
        ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaaggt                                                                                >  1:1648051/1‑89 (MQ=255)
        ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaaggg                                                                                >  1:1281245/1‑89 (MQ=255)
        ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaaggg                                                                                >  1:3520214/1‑89 (MQ=255)
        ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACaaaaaaggc                                                                                >  1:2981211/1‑89 (MQ=255)
          gAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGc                                                                               <  2:5436020/93‑1 (MQ=255)
           aaCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCg                                                                              <  2:6226801/93‑1 (MQ=255)
             cGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGt                                                                            >  2:3726451/1‑93 (MQ=255)
              gAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTc                                                                           <  1:1028988/93‑1 (MQ=255)
              gAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTc                                                                           >  2:2970458/1‑93 (MQ=255)
               aaaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCt                                                                          <  2:3213331/93‑1 (MQ=255)
                aaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAATTAGTGTAGTCGGCGTCACAAACAAGGTGCGGTCtt                                                                         >  1:38882/1‑93 (MQ=255)
                aaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCtt                                                                         <  1:938932/93‑1 (MQ=255)
                 aaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTa                                                                        >  2:1312561/1‑93 (MQ=255)
                  aTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTAc                                                                       <  1:4731186/93‑1 (MQ=255)
                        gtgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAtt                                                                 <  1:1528663/93‑1 (MQ=255)
                         tgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAttt                                                                <  1:2469543/93‑1 (MQ=255)
                           tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTc                                                              >  2:5767678/1‑93 (MQ=255)
                           tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTc                                                              <  2:692261/93‑1 (MQ=255)
                            ttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTcc                                                             <  1:3438041/93‑1 (MQ=255)
                              gATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGt                                                           <  1:5409839/93‑1 (MQ=255)
                              gATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGt                                                           >  2:5629333/1‑93 (MQ=255)
                              gATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGt                                                           >  2:3297459/1‑93 (MQ=255)
                              gATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGt                                                           >  1:5142988/1‑93 (MQ=255)
                                   ggCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGt                                                      >  2:5911544/1‑93 (MQ=255)
                                      aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                   >  1:6224034/1‑93 (MQ=255)
                                       tAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCgg                                                  >  2:3947067/1‑93 (MQ=255)
                                        aaTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGt                                                 >  2:1884187/1‑93 (MQ=255)
                                         aTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTa                                                >  1:1794428/1‑93 (MQ=255)
                                          ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               <  2:2732627/93‑1 (MQ=255)
                                          ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               >  2:2540387/1‑93 (MQ=255)
                                            ttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTGCGTAAAGTTCGGTActc                                             >  2:3480429/1‑93 (MQ=255)
                                            ttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActc                                             >  2:652098/1‑93 (MQ=255)
                                              ccTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTg                                           >  2:3784469/1‑93 (MQ=255)
                                               cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                          >  2:2059336/1‑93 (MQ=255)
                                               cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                          <  1:5623206/93‑1 (MQ=255)
                                               cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCAGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTACAGTTCGGTACTCTGa                                          >  2:921583/1‑93 (MQ=255)
                                                tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                         <  1:2850178/93‑1 (MQ=255)
                                                tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                         >  1:1407544/1‑93 (MQ=255)
                                                   cTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaa                                      <  1:633007/93‑1 (MQ=255)
                                                   cTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaa                                      <  1:5283310/93‑1 (MQ=255)
                                                   cTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaa                                      <  2:4936847/93‑1 (MQ=255)
                                                    tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaag                                     >  1:2888225/1‑93 (MQ=255)
                                                    tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaag                                     >  1:4578233/1‑93 (MQ=255)
                                                    tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaag                                     >  2:5282398/1‑93 (MQ=255)
                                                      tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   >  2:4123651/1‑93 (MQ=255)
                                                      tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   >  2:4315219/1‑93 (MQ=255)
                                                       gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTag                                  <  1:3333561/93‑1 (MQ=255)
                                                       gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTag                                  >  2:1376935/1‑93 (MQ=255)
                                                       gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGAGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTag                                  >  2:4904575/1‑93 (MQ=255)
                                                        ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTaga                                 <  2:1235499/93‑1 (MQ=255)
                                                        ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTaga                                 <  2:4378367/93‑1 (MQ=255)
                                                        ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTaga                                 <  2:4731916/93‑1 (MQ=255)
                                                         ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTagag                                >  2:5062959/1‑93 (MQ=255)
                                                         ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTagag                                >  1:2432058/1‑93 (MQ=255)
                                                          tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAgaga                               <  2:2751445/93‑1 (MQ=255)
                                                          tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAgaga                               >  2:4160333/1‑93 (MQ=255)
                                                          tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGGGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAgaga                               >  2:1158254/1‑93 (MQ=255)
                                                           ctcaTCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGAt                              >  2:4690409/4‑93 (MQ=255)
                                                            ttATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATa                             >  1:5229725/1‑93 (MQ=255)
                                                            ttATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATa                             <  2:1036174/93‑1 (MQ=255)
                                                             tATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATaa                            <  2:602224/93‑1 (MQ=255)
                                                              aTCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATaaa                           <  1:614211/93‑1 (MQ=255)
                                                               tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAt                          <  1:4511203/93‑1 (MQ=255)
                                                               tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAt                          >  1:4134670/1‑93 (MQ=255)
                                                                 tcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttc                        <  1:532440/93‑1 (MQ=255)
                                                                   aTTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                      <  1:1011952/93‑1 (MQ=255)
                                                                   aTTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                      <  1:5069732/93‑1 (MQ=255)
                                                                     tGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                    <  1:2099764/93‑1 (MQ=255)
                                                                     tGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                    >  2:5314600/1‑93 (MQ=255)
                                                                        aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  1:2651236/93‑1 (MQ=255)
                                                                        aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 >  2:6309209/1‑93 (MQ=255)
                                                                        aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  2:145762/93‑1 (MQ=255)
                                                                          tAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgaga               >  2:4998974/1‑93 (MQ=255)
                                                                          tAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgaga               <  2:1560145/93‑1 (MQ=255)
                                                                           aGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagag              >  1:823246/1‑93 (MQ=255)
                                                                           aGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagag              <  2:749738/93‑1 (MQ=255)
                                                                           aGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagag              <  1:2502780/93‑1 (MQ=255)
                                                                           aGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagag              <  2:1270380/93‑1 (MQ=255)
                                                                            gtgtAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagaga             <  1:116650/93‑1 (MQ=255)
                                                                              gtAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAg           >  1:1368953/1‑93 (MQ=255)
                                                                              gtAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAg           <  1:5539884/93‑1 (MQ=255)
                                                                              gtAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAg           <  2:4013102/93‑1 (MQ=255)
                                                                               tAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGc          >  2:4089593/1‑93 (MQ=255)
                                                                               tAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGc          >  2:6107703/1‑93 (MQ=255)
                                                                               tAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGc          >  1:1889502/1‑93 (MQ=255)
                                                                                aGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGcc         <  1:4806850/93‑1 (MQ=255)
                                                                                aGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGcc         <  2:1917208/93‑1 (MQ=255)
                                                                                 gTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa        >  1:3660614/1‑93 (MQ=255)
                                                                                  tCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCAt       <  2:1279766/93‑1 (MQ=255)
                                                                                    ggCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa     >  1:2300222/1‑93 (MQ=255)
                                                                                    ggCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa     <  2:4464450/93‑1 (MQ=255)
                                                                                     gCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGaa    >  2:3407/1‑93 (MQ=255)
                                                                                     gCGTCACAAAAAAGGTACGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGaa    <  2:3545694/93‑1 (MQ=255)
                                                                                       gTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGt  >  2:5642370/1‑93 (MQ=255)
                                                                                       gTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGt  <  2:5193800/93‑1 (MQ=255)
                                                                                            |                                                                                        
TGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGT  >  NC_000913/823681‑823859

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