Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,667 T→G 16.8% intergenic (+312/‑43) rnlB → / → yfjP CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,6670TG16.8% 365.9 / 72.4 160intergenic (+312/‑43)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (72/60);  new base G (17/10);  total (90/70)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 162 total aligned reads displayed.

TGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTA  >  NC_000913/2767569‑2767766
                                                                                                  |                                                                                                   
 ggTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcg                                                                                                   >  2:2704882/1‑100 (MQ=255)
 ggTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcg                                                                                                   >  1:2846746/1‑100 (MQ=255)
  gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                  <  1:2496825/100‑1 (MQ=255)
     cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGt                                                                                               <  1:244709/100‑1 (MQ=255)
      tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                              >  1:388498/1‑100 (MQ=255)
        tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                             >  2:3371284/1‑99 (MQ=255)
        tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                            >  1:2061689/1‑100 (MQ=255)
         gAGTTCAGAAAGATTCTGCACTCTATGTGCATTGAAAAAATCTTCAGATAGCTATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                            >  2:670592/1‑99 (MQ=255)
         gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                           >  2:1144706/1‑100 (MQ=255)
         gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTg                                                                                           <  1:580895/100‑1 (MQ=255)
          aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATGCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                          >  2:1674168/1‑100 (MQ=255)
          aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                          <  1:3764112/100‑1 (MQ=255)
             tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                       >  2:332664/1‑100 (MQ=255)
              ctgAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGtt                                                                                      <  1:2089570/98‑1 (MQ=255)
               aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                     >  2:1286385/1‑100 (MQ=255)
               aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttc                                                                                     <  2:2418038/100‑2 (MQ=255)
                gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACATCGCGGTTCTGCCGtt                                                                                      >  2:2381076/1‑98 (MQ=255)
                gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGACAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttag                                                                                    >  2:2089558/1‑98 (MQ=39)
                  aaGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTg                                                                                  >  1:1740002/1‑100 (MQ=255)
                   aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGt                                                                                 <  1:857935/100‑1 (MQ=255)
                    gATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtt                                                                                <  1:543282/100‑1 (MQ=255)
                        cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                            <  1:3216964/100‑1 (MQ=255)
                          gCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                          >  2:1826332/1‑100 (MQ=255)
                             aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                       <  1:2970343/100‑1 (MQ=255)
                             aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                       <  2:575718/100‑1 (MQ=255)
                              tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                      <  1:1842891/100‑1 (MQ=255)
                              tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                      <  1:2310113/100‑1 (MQ=255)
                                 aTGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCTATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCagag                                                                   <  2:1022631/100‑1 (MQ=255)
                                  tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                  >  2:1318073/1‑100 (MQ=255)
                                   gtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagga                                                                 >  2:1978172/1‑100 (MQ=255)
                                    tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                >  2:194917/1‑100 (MQ=255)
                                    tgCATAGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                  <  2:153059/98‑1 (MQ=255)
                                     gCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                >  1:3153345/1‑99 (MQ=255)
                                     gCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgt                                                               >  2:1250166/1‑100 (MQ=255)
                                      cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                              >  2:1188736/1‑100 (MQ=255)
                                          gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                          >  1:38073/1‑100 (MQ=255)
                                           aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                         >  2:2746010/1‑100 (MQ=255)
                                               aaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaaa                                                     <  2:3684400/100‑1 (MQ=255)
                                                  cTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAt                                                  >  2:249990/1‑100 (MQ=255)
                                                   ttCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTTCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtt                                                 >  2:2677272/1‑100 (MQ=255)
                                                    tCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAttt                                                <  1:3318581/100‑1 (MQ=255)
                                                     cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                               <  1:2529993/100‑1 (MQ=255)
                                                     cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                               <  1:2600977/100‑1 (MQ=255)
                                                     cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                               <  2:2435147/100‑1 (MQ=255)
                                                      agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                              >  1:359400/1‑100 (MQ=255)
                                                      agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                              >  2:2099458/1‑100 (MQ=255)
                                                      agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                              <  1:621857/100‑1 (MQ=255)
                                                      agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                              <  2:716232/100‑1 (MQ=255)
                                                       gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                             <  2:185347/100‑1 (MQ=255)
                                                        atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                            <  2:2194802/100‑1 (MQ=255)
                                                         tagataTCATTCTCCTGAACAATCCGATGAACGGAAGAACCCCGCGGTTCTGCCGTTTTTGGTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                           >  2:2080365/1‑100 (MQ=255)
                                                           gataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGtt                                         >  1:3666622/1‑100 (MQ=255)
                                                           gataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGATGATGAAAAATTTTGAAGtt                                         >  1:1567236/1‑100 (MQ=255)
                                                             tatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCt                                       >  1:1362666/1‑100 (MQ=255)
                                                               tCATTCTCCTGAACTATCCGATGTACGGGAGAACCTCGCGGTTCTGACTTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGATAATTTTTCAAGTTCTGc                                     <  2:2241340/100‑1 (MQ=255)
                                                               tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                     <  2:834778/100‑1 (MQ=255)
                                                                cATTCTCCTGAGCATTCCGATGAGCGGCAGAGCCTCGCGGTTGTGCCGTTTTTGTTGTACTTGTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                    <  2:3100764/100‑1 (MQ=255)
                                                                cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                    >  1:473398/1‑100 (MQ=255)
                                                                 aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                   >  2:1371287/1‑100 (MQ=255)
                                                                   tctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGcc                                 <  1:3248817/100‑1 (MQ=255)
                                                                      ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                              >  1:264295/1‑100 (MQ=255)
                                                                       cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTa                             >  1:335970/1‑100 (MQ=255)
                                                                         gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                            >  2:3718054/1‑99 (MQ=255)
                                                                           aCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTCTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGa                         <  2:2463699/100‑1 (MQ=255)
                                                                             aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                       <  2:1473053/100‑1 (MQ=255)
                                                                             aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                       <  1:1813518/100‑1 (MQ=255)
                                                                              aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                      >  2:2592831/1‑100 (MQ=255)
                                                                               tCCGATGAGCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTc                     <  1:187987/100‑1 (MQ=255)
                                                                               aaacaaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                      >  2:1441078/8‑99 (MQ=37)
                                                                                ccGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTcc                    >  1:1529076/1‑100 (MQ=255)
                                                                                  caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                  >  1:2771809/5‑100 (MQ=38)
                                                                                  caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                  >  1:3259139/5‑100 (MQ=38)
                                                                                    tGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt                <  2:736358/100‑1 (MQ=255)
                                                                                    aaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt                <  1:1529451/98‑1 (MQ=255)
                                                                                     gAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTAAGAGAAGTGGTGATAAAAAATTTTGAAGTTCTGCAGCCAATACAGAACTACCTTTc               >  2:3415102/1‑100 (MQ=255)
                                                                                     aaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTc               >  1:3581863/2‑100 (MQ=255)
                                                                                      aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTGCCTTTcc              >  2:2048620/1‑100 (MQ=255)
                                                                                      aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc              <  1:3462617/100‑1 (MQ=255)
                                                                                      aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc              >  1:2346080/1‑100 (MQ=255)
                                                                                      aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc              <  2:948118/100‑1 (MQ=255)
                                                                                       aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg             >  2:990873/1‑100 (MQ=255)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGCAGTTCTGCAGCCACTACAGAACTCCCTTTCCgn            <  2:2291704/100‑2 (MQ=255)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            <  2:890183/100‑1 (MQ=255)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTTTTTTTTTTTGTAGTGGTGATGAAAAATTTTGATGTTCTTCAGCCACTACAGAACTCCCTTTCCgt            >  2:2710786/1‑99 (MQ=25)
                                                                                         ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt           >  2:3022918/1‑100 (MQ=255)
                                                                                         ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt           >  1:52050/1‑100 (MQ=255)
                                                                                         ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt           <  1:808094/100‑1 (MQ=255)
                                                                                          ncAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  1:3045901/2‑100 (MQ=255)
                                                                                          gCAGACCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTAATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            <  2:2111623/98‑1 (MQ=255)
                                                                                          gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  2:1861460/1‑100 (MQ=255)
                                                                                          gCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  2:1515184/1‑100 (MQ=38)
                                                                                          gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  1:3404872/1‑100 (MQ=255)
                                                                                          gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  1:1061813/1‑100 (MQ=255)
                                                                                           cAGAACCTCGCTGTTCTGCCGTTGTTGTTTTACTGTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt         <  2:3301580/100‑1 (MQ=255)
                                                                                           cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt         >  2:2830366/1‑100 (MQ=255)
                                                                                           cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt         >  1:1570731/1‑100 (MQ=255)
                                                                                            aGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt           <  2:1656908/97‑1 (MQ=255)
                                                                                             gAACCGCGAGGTTCTGCCGTCTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTAc       >  2:2330168/1‑100 (MQ=255)
                                                                                               aCCTCGCGGTTCTGCCGTTTTTGTTTGACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCg     <  2:3284153/100‑1 (MQ=255)
                                                                                                 cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAAAAAAAAAAAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCt   <  2:3790215/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTA  >  NC_000913/2767569‑2767766

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