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

CATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCA  >  NC_000913/751632‑751828
                                                                                                    |                                                                                                  
cATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTaga                                                                                                     >  1:325047/1‑100 (MQ=255)
 agtCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTagag                                                                                                   >  1:1768422/3‑100 (MQ=255)
 aTTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTagag                                                                                                   >  1:2132769/1‑100 (MQ=255)
    cAGGTTAGTTCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggg                                                                                                <  2:2149266/100‑1 (MQ=255)
    cAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggg                                                                                                <  2:1129971/100‑1 (MQ=255)
    cAGGTTAGGCCTGCTGCTGCCGTCCCGTGGGATGTAGCTGCTTGCAAACAAGTTGTAAGCCAGCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggg                                                                                                <  2:1097795/100‑1 (MQ=255)
      ggTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggggg                                                                                              <  2:1968860/100‑1 (MQ=255)
      ggTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggggg                                                                                              <  2:1578844/100‑1 (MQ=255)
        tttgTGCTGCTGCTGCCGTCCTGTGTGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAg                                                                                           <  2:1736984/97‑1 (MQ=255)
           gTGCTGCTGCTGCCGTCTTGTGGGCGGTAGCTGCTGGCAAACGGGTTGTAATCCATCAGGACACTGGCAATACCTTCTTTCAATGTAGAggggggg                                                                                             >  2:270212/1‑95 (MQ=255)
            tgctgctgctgcCGTCCTGTCGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGCAGACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGccc                                                                                        <  2:2489101/100‑1 (MQ=255)
             gctgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCa                                                                                       >  1:1479578/1‑100 (MQ=255)
             gctgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAgtctgg                                                                                      <  2:1966215/100‑6 (MQ=255)
             gctgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAg                                                                                      >  1:2449591/1‑100 (MQ=255)
               tgctgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAg                                                                                      <  2:2120985/99‑1 (MQ=255)
                  tgctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGtttt                                                                                  <  2:239746/100‑1 (MQ=255)
                    ctgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCt                                                                                >  1:1415818/1‑100 (MQ=255)
                     tgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCtn                                                                               <  1:2123829/100‑2 (MQ=255)
                     tgcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGa                                                                              >  1:842/1‑100 (MQ=255)
                      gcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAg                                                                             >  2:1272332/1‑100 (MQ=255)
                      gcCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAg                                                                             >  1:1543074/1‑100 (MQ=255)
                       ccGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAg                                                                             >  1:670748/1‑100 (MQ=255)
                         gTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAAAACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTg                                                                           >  2:1665987/1‑100 (MQ=255)
                          tCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGcc                                                                         >  1:2027542/1‑100 (MQ=255)
                           ccTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGcc                                                                         >  1:2313569/1‑100 (MQ=255)
                             tgtgGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGACCAGTTTTCTGAGTGCCAc                                                                       >  1:2024590/1‑100 (MQ=255)
                              gtgGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCaag                                                                      >  2:2228495/1‑98 (MQ=255)
                                 ggCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACgcca                                                                   <  1:93608/100‑1 (MQ=255)
                                  gCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACgccag                                                                  <  2:2174377/100‑1 (MQ=255)
                                    ggTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACgccagcc                                                                <  2:1185530/100‑1 (MQ=255)
                                     gTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACgccagcca                                                               <  1:146457/100‑1 (MQ=255)
                                              ggCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCgg                                                      <  1:2133758/100‑1 (MQ=255)
                                                cAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                    >  1:1530857/1‑100 (MQ=255)
                                                cAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAAt                                                   >  1:1602499/1‑100 (MQ=255)
                                                 aaaCAAGTTGTAATCCATCAGGACACCTGCAACACCCTCTTTCCATGTAGAGGGGGGCAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                    <  1:1057960/99‑1 (MQ=255)
                                                 aaaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAAt                                                   >  1:1399873/1‑100 (MQ=255)
                                                  aaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATa                                                  >  2:1254344/1‑100 (MQ=255)
                                                  aaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATa                                                  >  1:73237/1‑100 (MQ=255)
                                                  aaCAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATAc                                                 <  2:1246501/100‑1 (MQ=255)
                                                      aGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTaa                                              >  1:2482115/1‑100 (MQ=255)
                                                      aGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAat                                             >  2:1654439/1‑100 (MQ=255)
                                                       gTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAata                                            >  1:1293626/1‑100 (MQ=255)
                                                        ttGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAata                                            <  1:1958663/100‑1 (MQ=255)
                                                           tAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTc                                         <  2:73237/100‑1 (MQ=255)
                                                            aaTCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAg                                       <  2:818542/100‑1 (MQ=255)
                                                             aTCCATCAGGACACCGGCAACACCTTCTTTCAATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAg                                       >  2:36423/1‑99 (MQ=255)
                                                              tCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGCATACTAATATTCAGtt                                     >  2:292389/1‑100 (MQ=255)
                                                               ccATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  2:2398370/1‑100 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  2:2464647/1‑100 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    <  2:1189987/100‑1 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  1:572689/1‑100 (MQ=255)
                                                                cATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTg                                    >  1:1189992/1‑100 (MQ=255)
                                                                 aTCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCt                                  >  2:1321696/1‑100 (MQ=255)
                                                                  tCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTg                                 <  1:545825/100‑1 (MQ=255)
                                                                    aGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGa                                >  2:2263384/1‑100 (MQ=255)
                                                                     ggACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAATGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGAt                               >  2:2337420/1‑100 (MQ=255)
                                                                      gACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGAtt                              >  1:315463/1‑100 (MQ=255)
                                                                       acacCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTg                             <  1:2551695/100‑1 (MQ=255)
                                                                       acacCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTg                             <  1:622701/100‑1 (MQ=255)
                                                                        cacCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTgg                            <  1:640432/100‑1 (MQ=255)
                                                                         acCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          <  1:2075910/100‑1 (MQ=255)
                                                                          ccGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          >  1:2220316/1‑100 (MQ=255)
                                                                          ccGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          >  1:884713/1‑100 (MQ=255)
                                                                          ccGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCtt                         <  2:1655118/100‑1 (MQ=255)
                                                                          ccGGAAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCt                          <  2:581223/100‑1 (MQ=255)
                                                                              cAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGAt                      >  1:203662/1‑100 (MQ=255)
                                                                              cAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCAGATTGGCTTGATc                     <  1:1749571/100‑1 (MQ=255)
                                                                               aaCACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATc                     <  2:71571/100‑1 (MQ=255)
                                                                               aaCACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATc                     <  2:93620/100‑1 (MQ=255)
                                                                               aaCACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATc                     <  2:1247835/100‑1 (MQ=255)
                                                                                acacCTTCTTGGCATGTAGAGGGGGGGAGCCCAGTTTTATGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTAATGATTGGCTTGATCg                    >  2:477845/1‑100 (MQ=255)
                                                                                 cacCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaa                  >  2:68926/1‑100 (MQ=255)
                                                                                 cacCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaa                  >  2:1143704/1‑100 (MQ=255)
                                                                                  acCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaaa                 <  1:1404726/100‑1 (MQ=255)
                                                                                   ccTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGTGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGaaa                 >  2:805376/1‑100 (MQ=255)
                                                                                    cttcttTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAAtt               >  1:1577518/1‑100 (MQ=255)
                                                                                    cttcttTCCAAGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTCGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAAt                >  2:1831683/1‑100 (MQ=255)
                                                                                      tcttTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTg              >  1:664958/1‑100 (MQ=255)
                                                                                       cttTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGa             <  2:2220304/100‑1 (MQ=255)
                                                                                        tttCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAag           >  1:2464368/1‑100 (MQ=255)
                                                                                        tttCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAag           <  2:2449578/100‑1 (MQ=255)
                                                                                         ttCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAag           >  2:79193/1‑100 (MQ=255)
                                                                                          tCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAaga          <  2:2217625/100‑1 (MQ=255)
                                                                                            cATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       >  1:1526517/1‑100 (MQ=255)
                                                                                            cATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       <  2:2027531/100‑1 (MQ=255)
                                                                                             aTGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       >  1:1779395/1‑100 (MQ=255)
                                                                                             aTGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCa       <  2:1937167/100‑1 (MQ=255)
                                                                                             aTGTAGAGGGGGGGAGCCCAGTTTTCGGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGc        <  1:2604888/99‑1 (MQ=255)
                                                                                              tGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCAt      >  1:2227565/1‑100 (MQ=255)
                                                                                                tAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTc   <  1:1272338/100‑1 (MQ=255)
                                                                                                 agagGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAAAATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCa  <  1:1112671/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
CATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCA  >  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: