Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 750,896 Δ1 bp 100% coding (1274/2448 nt) ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913750,8930C.99.0% 446.6 / ‑1.8 97coding (1277/2448 nt)ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (49/47);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 104 total aligned reads displayed.

TAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCA  >  NC_000913/750794‑750986
                                                                                                   |                                                                                             
tAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGc                                                                                               <  1:1091628/100‑1 (MQ=255)
 aGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcc                                                                                              <  2:1918294/100‑1 (MQ=255)
 aGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcc                                                                                              <  2:2163507/100‑1 (MQ=255)
   cTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccg                                                                                            >  1:853634/1‑99 (MQ=255)
     gtgAAAATGACGATCTGAGAAGCGATAAGTGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgct                                                                                          <  2:1725465/100‑4 (MQ=255)
     gtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgct                                                                                          >  2:2441256/1‑97 (MQ=255)
        aaaaTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGt                                                                                      >  2:1088414/1‑100 (MQ=255)
        aaaaTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGt                                                                                      >  2:70032/1‑100 (MQ=255)
           aTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCct                                                                                   >  2:2441047/1‑100 (MQ=255)
           aTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCct                                                                                   <  1:716149/100‑1 (MQ=255)
               cGATCTGAGAAGCGATAAGCTGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctgctg                                                                               <  2:1380923/100‑1 (MQ=255)
               cGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctgctg                                                                               <  1:1346021/100‑1 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                             <  1:786400/100‑1 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                             <  2:3004618/100‑1 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                             <  1:260371/100‑1 (MQ=255)
                    tGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTAt                                                                          >  1:89123/1‑100 (MQ=255)
                     gagaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATc                                                                         >  2:2007652/1‑100 (MQ=255)
                      agaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCa                                                                        <  2:11208/100‑1 (MQ=255)
                       gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCACCGCTCGTCCTGCTGAGTATCa                                                                        >  1:124894/1‑100 (MQ=255)
                       gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                       <  1:2214446/100‑1 (MQ=255)
                       gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                       >  1:1440084/1‑100 (MQ=255)
                       gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                       >  1:656014/1‑100 (MQ=255)
                       gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCACGCTCGTCCTGCTGAGTATCaa                                                                       >  2:2783076/1‑100 (MQ=255)
                       gaAGCGATAAGCGGCGAGCAAAATCGTGAAGATAATGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                       >  2:454281/1‑100 (MQ=255)
                         aGCGATAAGCGGCGAGCGAAGTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaaaa                                                                     >  2:1852083/1‑100 (MQ=255)
                          gCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAAt                                                                    <  2:2245486/100‑1 (MQ=255)
                          gCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAAt                                                                    <  1:1336261/100‑1 (MQ=255)
                          gCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGATCGTCCTGCTGAGTATCAACAt                                                                    >  1:1003225/1‑100 (MQ=255)
                           cGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATg                                                                   <  1:2222861/100‑1 (MQ=255)
                             aTAAGCGGCGAGCGAAATTGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                 <  1:66686/100‑1 (MQ=255)
                                aGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTAc                                                              <  1:2034938/100‑1 (MQ=255)
                                 gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACt                                                             >  2:2087726/1‑100 (MQ=255)
                                   ggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTgg                                                           <  2:758932/100‑1 (MQ=255)
                                   ggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTgg                                                           <  2:49569/100‑1 (MQ=255)
                                    gcgagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGc                                                          <  1:708841/100‑1 (MQ=255)
                                    gcgagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGc                                                          <  1:1653839/100‑1 (MQ=255)
                                      gagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGccc                                                        <  2:1377518/100‑1 (MQ=255)
                                      gagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGccc                                                        <  1:2703622/100‑1 (MQ=255)
                                       agcgaAATCGTGCTGTTGGTGGCATACACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCACGCTCGTCATGCTGAGTATAAAAATGGATACTgg                                                           >  2:2569146/1‑96 (MQ=255)
                                         cgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCa                                                       >  1:748791/1‑98 (MQ=255)
                                         cgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAgg                                                     <  1:2263831/100‑1 (MQ=255)
                                          gaatttgTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTAAAGTCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCa                                                       <  2:1751899/91‑1 (MQ=255)
                                           aaaTCGTGCTGTTGGTGGCATCCACTTGTTTTCTGTAATTAAACCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATTGCTACTGGCCCAGGt                                                    <  2:2361366/99‑1 (MQ=255)
                                           aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                   <  1:763788/100‑1 (MQ=255)
                                           aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                   >  2:1342810/1‑100 (MQ=255)
                                           aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                   >  1:191772/1‑100 (MQ=255)
                                           aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTGAAGCCCGCTTGTCCTCCTGAGTATAAAAATGGCTACTGGCCCa                                                       >  2:2048067/1‑96 (MQ=255)
                                            aaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGACCAGGTg                                                   <  2:1930644/99‑1 (MQ=255)
                                             aTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                   <  1:681528/98‑1 (MQ=255)
                                              tCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                <  2:2925035/100‑1 (MQ=255)
                                                gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGTTAGCTTAAGCCCGATCTTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                              <  2:1357892/100‑1 (MQ=255)
                                                 tGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCa                                             >  1:865874/1‑100 (MQ=255)
                                                  gCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaa                                            <  1:2014802/100‑1 (MQ=255)
                                                   cTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaaa                                           >  2:2697371/1‑100 (MQ=255)
                                                   cTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaaa                                           >  2:709810/1‑100 (MQ=255)
                                                    tGTTGGTGGCATCCACTTGTTTGCTGTGCTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTTAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAg                                          <  2:3035465/100‑1 (MQ=255)
                                                    tGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAg                                          <  1:1842114/100‑1 (MQ=255)
                                                     gTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAgg                                         >  1:2777598/1‑100 (MQ=255)
                                                        ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACa                                      >  1:357627/1‑100 (MQ=255)
                                                        ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACa                                      >  1:887662/1‑100 (MQ=255)
                                                        ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACa                                      >  2:2925994/1‑100 (MQ=255)
                                                         gtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAg                                     >  2:2755214/1‑100 (MQ=255)
                                                          tggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGACCCAGGTGACATCAAAGGACAg                                     >  1:728839/1‑99 (MQ=255)
                                                           ggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGACCAGGTGACATCAAAGGACAGTg                                   <  1:644227/100‑1 (MQ=255)
                                                              aTCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTgcgc                                >  1:71126/1‑100 (MQ=255)
                                                                ccACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCa                              >  2:768959/1‑100 (MQ=255)
                                                                 cACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTAGCCCAGGTGACATCAAAGGACAGTGCGCCaa                             >  1:386900/1‑100 (MQ=255)
                                                                   cTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGACTACTGACCCAGGTGACATCAAAGGACAGTGCGCCAAGc                           >  1:2321097/1‑100 (MQ=255)
                                                                    ttGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGTTCGTCTTGTTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGcc                          <  1:1025792/100‑1 (MQ=255)
                                                                    ttGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGcc                          <  2:382051/100‑1 (MQ=255)
                                                                      gTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcac                        <  2:292693/100‑1 (MQ=255)
                                                                      gTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAAGTAACATCAAAGGACAGTGCGCCAAGCcac                        >  1:2408890/1‑100 (MQ=255)
                                                                      gTTTGCTCTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcac                        >  1:1583698/1‑100 (MQ=255)
                                                                      gTTTGCTATAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcac                        >  1:2457968/1‑100 (MQ=255)
                                                                        ttGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAg                      <  1:1146328/100‑1 (MQ=255)
                                                                           cTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                   >  2:488895/1‑100 (MQ=255)
                                                                            tGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTAGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCata                  >  1:802389/1‑100 (MQ=255)
                                                                             gTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCatat                 >  2:2500200/1‑100 (MQ=255)
                                                                             gTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCatat                 >  1:72766/1‑100 (MQ=255)
                                                                              tAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtt                >  2:1938400/1‑100 (MQ=255)
                                                                               aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt               <  1:70032/100‑1 (MQ=255)
                                                                               aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt               >  2:79221/1‑100 (MQ=255)
                                                                               aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt               <  1:2360047/100‑1 (MQ=255)
                                                                               aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt               >  1:3003608/1‑100 (MQ=255)
                                                                               aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt               >  2:1210300/1‑100 (MQ=255)
                                                                                aTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttgt              >  2:1836625/1‑98 (MQ=255)
                                                                                 ttAAAACGGTAGCGTAAGCCCGCTCGTCCTGCTGAGTAACAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTg             >  2:1872805/1‑100 (MQ=255)
                                                                                  tAAAACGTTAGCTTAAGCCCGCTCGTCTTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc            <  1:735580/100‑1 (MQ=255)
                                                                                  tAAAACGGTAGCTTAAGCTCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc            <  1:666361/100‑1 (MQ=255)
                                                                                  tAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc            <  2:2437290/100‑1 (MQ=255)
                                                                                  tAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc            >  1:1796420/1‑100 (MQ=255)
                                                                                  tAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGGGCGCCAAGCCACAGCATATTTTGc            >  2:2958613/1‑100 (MQ=255)
                                                                                   aaaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc           <  2:1227560/100‑1 (MQ=255)
                                                                                   aaaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc           <  2:1440070/100‑1 (MQ=255)
                                                                                   aaaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc           <  2:2221340/100‑1 (MQ=255)
                                                                                    aaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGccc          <  2:1692020/100‑1 (MQ=255)
                                                                                    aaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGACCAGGTGACATCAAAGGACAGTACGCCAAGCCACAGCATATTTTGccc          >  1:1382467/1‑100 (MQ=255)
                                                                                      aCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGa        >  2:1381236/1‑100 (MQ=255)
                                                                                          tAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAcc    >  1:195771/1‑100 (MQ=255)
                                                                                            gCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCa  <  2:2631604/100‑1 (MQ=255)
                                                                                                   |                                                                                             
TAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCA  >  NC_000913/750794‑750986

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