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.67.2% 179.0 / 212.2 131coding (1277/2448 nt)ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (27/16);  new base . (45/43);  total (72/59)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 140 total aligned reads displayed.

TAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGG  >  NC_000913/750794‑750992
                                                                                                   |                                                                                                   
tAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGc                                                                                                     >  2:2268164/1‑100 (MQ=255)
tAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGc                                                                                                     >  2:2266131/1‑100 (MQ=255)
tAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGc                                                                                                     <  2:1733481/100‑1 (MQ=255)
 aGCTGTGAAAATGGCGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcc                                                                                                    <  1:3097362/100‑1 (MQ=255)
 aGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcc                                                                                                    <  2:683000/100‑1 (MQ=255)
  gCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAGATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGccc                                                                                                   <  2:413072/100‑1 (MQ=255)
           aTGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCct                                                                                         >  1:2604792/1‑100 (MQ=255)
            tGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTCAGCCCCGCTCGTCct                                                                                         >  2:3345008/1‑100 (MQ=255)
            tGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCct                                                                                         >  1:3110751/1‑100 (MQ=255)
            tGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCctg                                                                                        <  2:3293250/100‑1 (MQ=255)
             gACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCctg                                                                                        >  2:1109949/1‑100 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                                   <  1:1650284/100‑1 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                                   <  2:951399/100‑1 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                                   <  1:377492/100‑1 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                                   <  2:2461326/100‑1 (MQ=255)
                 aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAg                                                                                   <  2:1519882/100‑1 (MQ=255)
                  tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAg                                                                                   <  1:194472/100‑1 (MQ=255)
                  tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGt                                                                                  <  2:1687497/100‑1 (MQ=255)
                  tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGt                                                                                  <  2:447834/100‑1 (MQ=255)
                    tGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTAt                                                                                >  2:3127757/1‑100 (MQ=255)
                     gagaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTAt                                                                                >  1:1258705/1‑100 (MQ=255)
                     gagaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATc                                                                               <  1:1218431/100‑1 (MQ=255)
                      agaagaGATAAGCCGCGAGCTAAACCGTGTTGTTCGTGGCATCCACATGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTAGCa                                                                              <  1:1003433/100‑1 (MQ=255)
                      agaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATc                                                                               >  1:2474386/1‑100 (MQ=255)
                      agaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCa                                                                              >  1:1448187/1‑100 (MQ=255)
                      agaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGACATCCACGTGTTTGCTATAATTCAAACGGTAGCCTAAGCCCCGCCCGTCCTGCTGAGTATc                                                                               >  1:901155/1‑100 (MQ=255)
                       gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCaa                                                                             >  2:1529217/1‑100 (MQ=255)
                          gCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCaaaa                                                                           <  2:2932005/100‑1 (MQ=255)
                           cGATAAGCGGCGGGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATg                                                                         <  1:2700063/100‑1 (MQ=255)
                           cGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATg                                                                         <  1:3008406/100‑1 (MQ=255)
                             aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                       >  1:2194483/1‑100 (MQ=255)
                              tAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                       >  2:3258731/1‑100 (MQ=255)
                                 gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTAc                                                                    >  2:981034/1‑100 (MQ=255)
                                  cggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACt                                                                   >  1:1509099/1‑100 (MQ=255)
                                  cggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTAc                                                                    >  1:2135635/1‑98 (MQ=255)
                                   ggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTgg                                                                 <  1:374894/100‑1 (MQ=255)
                                   ggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTgg                                                                 >  1:2668052/1‑100 (MQ=255)
                                   ggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTgg                                                                 <  1:1038576/100‑1 (MQ=255)
                                     cgagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGcc                                                               <  1:1721721/100‑1 (MQ=255)
                                     cgagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGcc                                                               >  2:436867/1‑100 (MQ=255)
                                      gagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGccc                                                              >  2:1582522/1‑100 (MQ=255)
                                        gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCa                                                             >  2:1089250/1‑100 (MQ=255)
                                        gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                            >  2:748456/1‑100 (MQ=255)
                                        gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                            >  1:2401572/1‑100 (MQ=255)
                                         cgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAgg                                                           <  2:507168/100‑1 (MQ=255)
                                           aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAgg                                                           >  1:2686320/1‑99 (MQ=255)
                                            aaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTACGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                            <  1:1654102/97‑1 (MQ=255)
                                            aaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                         <  2:267014/100‑1 (MQ=255)
                                            aaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCCACTGGCCCAggag                                                         <  2:2294843/100‑3 (MQ=255)
                                             aTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGa                                                        <  1:222939/100‑1 (MQ=255)
                                             aTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGAc                                                       <  2:2959582/100‑1 (MQ=255)
                                               cGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                      <  2:1554824/100‑1 (MQ=255)
                                                gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                                    >  2:1859951/1‑100 (MQ=255)
                                                 tGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCa                                                   >  2:2652307/1‑100 (MQ=255)
                                                 tGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCa                                                   <  2:437617/100‑1 (MQ=255)
                                                  gCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCaa                                                  >  2:1717873/1‑100 (MQ=255)
                                                    tGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAg                                                <  2:2753727/100‑1 (MQ=255)
                                                    tGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAg                                                >  1:174972/1‑100 (MQ=255)
                                                    tGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAg                                                <  2:1552198/100‑1 (MQ=255)
                                                     gTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAg                                                <  2:559342/100‑1 (MQ=255)
                                                     gTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAgg                                               <  2:111551/100‑1 (MQ=255)
                                                     gTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAgg                                               <  1:995009/100‑1 (MQ=255)
                                                       tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGa                                              >  1:450796/1‑100 (MQ=255)
                                                       tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                             >  2:1419497/1‑100 (MQ=255)
                                                        ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGCCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                             >  1:1723060/1‑100 (MQ=255)
                                                         gtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAg                                           >  1:1139832/1‑100 (MQ=255)
                                                           ggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGt                                          >  2:2067216/1‑100 (MQ=255)
                                                           ggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTg                                         <  2:2550767/100‑1 (MQ=255)
                                                            gCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTgc                                        <  1:3141771/100‑1 (MQ=255)
                                                                ccACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGcc                                     <  1:3166046/100‑1 (MQ=255)
                                                                    ttGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGc                                 <  1:229802/100‑1 (MQ=255)
                                                                           cTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                         >  1:278709/1‑100 (MQ=255)
                                                                           cTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                         >  2:2026330/1‑100 (MQ=255)
                                                                           cTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                         >  1:3240243/1‑100 (MQ=255)
                                                                           cTGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                         >  2:247746/1‑100 (MQ=255)
                                                                            tGTAATTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCACAGGACAGTGCGCCAAGCCACAGCata                        >  2:644316/1‑100 (MQ=255)
                                                                               aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt                     <  1:2797377/100‑1 (MQ=255)
                                                                               aaTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt                     >  2:2604114/1‑100 (MQ=255)
                                                                                aTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCATGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt                    <  2:1612583/100‑1 (MQ=255)
                                                                                aTTAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt                    >  2:3304077/1‑100 (MQ=255)
                                                                                 ttAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTg                   >  2:2647018/1‑100 (MQ=255)
                                                                                 ttAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTg                   >  2:2519468/1‑100 (MQ=255)
                                                                                  tAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc                  <  2:237149/100‑1 (MQ=255)
                                                                                  tAAAACGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc                  <  1:208487/100‑1 (MQ=255)
                                                                                   aaaaCGGTAGCTTAAGCCCGCTCGTCGTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc                  <  1:406435/99‑1 (MQ=255)
                                                                                   aaaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc                 >  1:962865/1‑100 (MQ=255)
                                                                                    aaaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGccc                >  1:2491174/1‑100 (MQ=255)
                                                                                     aaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCg               <  1:2592548/100‑1 (MQ=255)
                                                                                     aaCGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCg               >  2:1408062/1‑100 (MQ=255)
                                                                                      aCGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCg               >  1:2034248/1‑100 (MQ=255)
                                                                                       cGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGAt             <  2:2744856/100‑1 (MQ=255)
                                                                                       cGGTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGAt             <  2:1242020/100‑1 (MQ=255)
                                                                                        ggTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGAt             <  1:908217/100‑1 (MQ=255)
                                                                                         gTAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAc           <  1:102465/100‑1 (MQ=255)
                                                                                          tAGCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAcc          >  2:2963913/1‑100 (MQ=255)
                                                                                            gCTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCa        >  2:2723666/1‑100 (MQ=255)
                                                                                             cTTAAGCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATACTTTGCCCGATACCCAt       >  2:3099941/1‑100 (MQ=255)
                                                                                              ttAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCAt       >  1:440300/1‑100 (MQ=255)
                                                                                                aaGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTg     <  2:1371978/100‑1 (MQ=255)
                                                                                                   ccccGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCgg  <  2:1275425/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGG  >  NC_000913/750794‑750992

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