Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 751,737:1 (G)6→7 100% coding (433/2448 nt) ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913751,7311.G63.0% 171.4 / 123.9 81coding (439/2448 nt)ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (17/13);  new base G (28/23);  total (45/36)
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.99e-01

CATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCA  >  NC_000913/751632‑751828
                                                                                                    |                                                                                                  
cATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTaga                                                                                                     >  1:325047/1‑100 (MQ=255)
 agtCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTaga‑g                                                                                                   >  1:1768422/3‑100 (MQ=255)
 aTTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTaga‑g                                                                                                   >  1:2132769/1‑100 (MQ=255)
    cAGGTTAGTTCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑gggg                                                                                                <  2:2149266/100‑1 (MQ=255)
    cAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑gggg                                                                                                <  2:1129971/100‑1 (MQ=255)
    cAGGTTAGGCCTGCTGCTGCCGTCCCGTGGGATGTAGCTGCTTGCAAACAAGTTGTAAGCCAGCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑gggg                                                                                                <  2:1097795/100‑1 (MQ=255)
      ggTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑gggggg                                                                                              <  2:1968860/100‑1 (MQ=255)
      ggTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑gggggg                                                                                              <  2:1578844/100‑1 (MQ=255)
        tttgTGCTGCTGCTGCCGTCCTGTGTGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑Ag                                                                                           <  2:1736984/97‑1 (MQ=255)
           gTGCTGCTGCTGCCGTCTTGTGGGCGGTAGCTGCTGGCAAACGGGTTGTAATCCATCAGGACACTGGCAATACCTTCTTTCAATGTAGA‑ggggggg                                                                                             >  2:270212/1‑95 (MQ=255)
            tgctgctgctgcCGTCCTGTCGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGCAGACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGccc                                                                                        <  2:2489101/100‑1 (MQ=255)
             gctgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCa                                                                                       >  1:1479578/1‑100 (MQ=255)
             gctgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑Agtctgg                                                                                      <  2:1966215/100‑6 (MQ=255)
             gctgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAg                                                                                      >  1:2449591/1‑100 (MQ=255)
               tgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAg                                                                                      <  2:2120985/99‑1 (MQ=255)
                  tgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGtttt                                                                                  <  2:239746/100‑1 (MQ=255)
                    ctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCt                                                                                >  1:1415818/1‑100 (MQ=255)
                     tgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCtn                                                                               <  1:2123829/100‑2 (MQ=255)
                     tgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGa                                                                              >  1:842/1‑100 (MQ=255)
                      gcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAg                                                                             >  2:1272332/1‑100 (MQ=255)
                      gcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAg                                                                             >  1:1543074/1‑100 (MQ=255)
                       ccGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAg                                                                             >  1:670748/1‑100 (MQ=255)
                         gTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAAAACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTg                                                                           >  2:1665987/1‑100 (MQ=255)
                          tCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGcc                                                                         >  1:2027542/1‑100 (MQ=255)
                           ccTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGcc                                                                         >  1:2313569/1‑100 (MQ=255)
                             tgtgGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGACCAGTTTTCTGAGTGCCAc                                                                       >  1:2024590/1‑100 (MQ=255)
                              gtgGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCaag                                                                      >  2:2228495/1‑98 (MQ=255)
                                 ggCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACgcca                                                                   <  1:93608/100‑1 (MQ=255)
                                  gCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACgccag                                                                  <  2:2174377/100‑1 (MQ=255)
                                    ggTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACgccagcc                                                                <  2:1185530/100‑1 (MQ=255)
                                     gTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACgccagcca                                                               <  1:146457/100‑1 (MQ=255)
                                              ggCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCgg                                                      <  1:2133758/100‑1 (MQ=255)
                                                cAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                    >  1:1530857/1‑100 (MQ=255)
                                                cAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAAt                                                   >  1:1602499/1‑100 (MQ=255)
                                                 aaaCAAGTTGTAATCCATCAGGACACCTGCAACACCCTCTTTCCATGTAGA‑GGGGGGCAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                    <  1:1057960/99‑1 (MQ=255)
                                                 aaaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAAt                                                   >  1:1399873/1‑100 (MQ=255)
                                                  aaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATa                                                  >  2:1254344/1‑100 (MQ=255)
                                                  aaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATa                                                  >  1:73237/1‑100 (MQ=255)
                                                  aaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATAc                                                 <  2:1246501/100‑1 (MQ=255)
                                                      aGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTaa                                              >  1:2482115/1‑100 (MQ=255)
                                                      aGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAat                                             >  2:1654439/1‑100 (MQ=255)
                                                       gTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAata                                            >  1:1293626/1‑100 (MQ=255)
                                                        ttGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAata                                            <  1:1958663/100‑1 (MQ=255)
                                                           tAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTc                                         <  2:73237/100‑1 (MQ=255)
                                                            aaTCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAg                                       <  2:818542/100‑1 (MQ=255)
                                                             aTCCATCAGGACACCGGCAACACCTTCTTTCAATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAg                                       >  2:36423/1‑99 (MQ=255)
                                                              tCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGCATACTAATATTCAGtt                                     >  2:292389/1‑100 (MQ=255)
                                                               ccATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  2:2398370/1‑100 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  2:2464647/1‑100 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    <  2:1189987/100‑1 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  1:572689/1‑100 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  1:1189992/1‑100 (MQ=255)
                                                                 aTCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCt                                  >  2:1321696/1‑100 (MQ=255)
                                                                  tCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTg                                 <  1:545825/100‑1 (MQ=255)
                                                                    aGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGa                                >  2:2263384/1‑100 (MQ=255)
                                                                     ggACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAATGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGAt                               >  2:2337420/1‑100 (MQ=255)
                                                                      gACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGAtt                              >  1:315463/1‑100 (MQ=255)
                                                                       acacCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTg                             <  1:2551695/100‑1 (MQ=255)
                                                                       acacCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTg                             <  1:622701/100‑1 (MQ=255)
                                                                        cacCGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTgg                            <  1:640432/100‑1 (MQ=255)
                                                                         acCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          <  1:2075910/100‑1 (MQ=255)
                                                                          ccGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          >  1:2220316/1‑100 (MQ=255)
                                                                          ccGGCAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          >  1:884713/1‑100 (MQ=255)
                                                                          ccGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCtt                         <  2:1655118/100‑1 (MQ=255)
                                                                          ccGGAAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          <  2:581223/100‑1 (MQ=255)
                                                                              cAACACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGAt                      >  1:203662/1‑100 (MQ=255)
                                                                              cAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCAGATTGGCTTGATc                     <  1:1749571/100‑1 (MQ=255)
                                                                               aaCACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATc                     <  2:71571/100‑1 (MQ=255)
                                                                               aaCACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATc                     <  2:93620/100‑1 (MQ=255)
                                                                               aaCACCTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATc                     <  2:1247835/100‑1 (MQ=255)
                                                                                acacCTTCTTGGCATGTAGAGGGGGGG‑AGCCCAGTTTTATGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTAATGATTGGCTTGATCg                    >  2:477845/1‑100 (MQ=255)
                                                                                 cacCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaa                  >  2:68926/1‑100 (MQ=255)
                                                                                 cacCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaa                  >  2:1143704/1‑100 (MQ=255)
                                                                                  acCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaaa                 <  1:1404726/100‑1 (MQ=255)
                                                                                   ccTTCTTTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGTGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaaa                 >  2:805376/1‑100 (MQ=255)
                                                                                    cttcttTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAAtt               >  1:1577518/1‑100 (MQ=255)
                                                                                    cttcttTCCAAGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTCGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAAt                >  2:1831683/1‑100 (MQ=255)
                                                                                      tcttTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTg              >  1:664958/1‑100 (MQ=255)
                                                                                       cttTCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGa             <  2:2220304/100‑1 (MQ=255)
                                                                                        tttCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAag           >  1:2464368/1‑100 (MQ=255)
                                                                                        tttCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAag           <  2:2449578/100‑1 (MQ=255)
                                                                                         ttCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAag           >  2:79193/1‑100 (MQ=255)
                                                                                          tCCATGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAaga          <  2:2217625/100‑1 (MQ=255)
                                                                                            cATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       >  1:1526517/1‑100 (MQ=255)
                                                                                            cATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       <  2:2027531/100‑1 (MQ=255)
                                                                                             aTGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       >  1:1779395/1‑100 (MQ=255)
                                                                                             aTGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       <  2:1937167/100‑1 (MQ=255)
                                                                                             aTGTAGAGGGGGGG‑AGCCCAGTTTTCGGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGc        <  1:2604888/99‑1 (MQ=255)
                                                                                              tGTAGAGGGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCAt      >  1:2227565/1‑100 (MQ=255)
                                                                                                tAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTc   <  1:1272338/100‑1 (MQ=255)
                                                                                                 aga‑gGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAAAATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCa  <  1:1112671/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
CATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGA‑GGGGGG‑AGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCA  >  NC_000913/751632‑751828

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