Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,667 T→G 15.9% 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,6670TG15.9% 367.7 / 61.0 148intergenic (+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 (61/63);  new base G (16/8);  total (77/71)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 148 total aligned reads displayed.

TGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCT  >  NC_000913/2767569‑2767765
                                                                                                  |                                                                                                  
tGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTc                                                                                                   >  2:1155506/1‑100 (MQ=255)
tGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTc                                                                                                   <  1:1719001/100‑1 (MQ=255)
   tGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgcg                                                                                                <  1:745293/100‑1 (MQ=255)
   tGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgcg                                                                                                <  2:2345200/100‑1 (MQ=255)
    gCTGTGAGTTCAGAAAGATTCTGTAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                               <  1:2941983/100‑1 (MQ=255)
    gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                               <  1:3245407/100‑1 (MQ=255)
     cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                              <  1:1332463/100‑1 (MQ=255)
             tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                      >  1:1170758/1‑100 (MQ=255)
             tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                      <  2:2016980/100‑1 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                     >  2:1730818/1‑100 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                     >  1:3068410/1‑100 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                     >  2:1572585/1‑100 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATGGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                      <  2:760424/99‑1 (MQ=39)
                taaaGAGTCTACAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATAACCGGTATAACCTCGCGGTTCTGTCGtttt                                                                                   <  2:900492/99‑1 (MQ=255)
                gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGTCAGAACCTCGCGGTTCTGCCGtttt                                                                                   >  2:14151/1‑100 (MQ=255)
                gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                   >  1:2116099/1‑100 (MQ=255)
                gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                   >  2:2522524/1‑100 (MQ=255)
                gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                   >  1:2919721/1‑100 (MQ=255)
                gAAAGATTCTGCAATCTATGTGCATTGAAAAAACCTTCAGATAGATATCATTCTCCTGAAAAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGt                                                                                      >  1:3272574/1‑97 (MQ=255)
                  aaGATTCTGCAATCTATGTGCAGTGAAAAAAGCTTCAGAGAGATATCATTGTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTg                                                                                 <  2:3602572/100‑1 (MQ=255)
                   aGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttcatcg                                                                                <  1:2165849/100‑6 (MQ=38)
                     aTTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGttt                                                                              <  2:2735574/100‑1 (MQ=255)
                     aTTCTGCAATCTATGCGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttcatcgg                                                                               <  2:2511162/99‑7 (MQ=37)
                       tCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTa                                                                            >  1:1024825/1‑100 (MQ=255)
                        cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                           >  1:1752264/1‑100 (MQ=255)
                          gCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTGTCTTGAACAATCCGATGAACGGCAGAATCTCGTTGCTCTGCCGTCTTTGTTTTa                                                                            >  2:2908719/1‑97 (MQ=38)
                          gCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                         <  1:2503229/100‑1 (MQ=255)
                             aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                      <  2:3549320/100‑1 (MQ=255)
                             aTCTATGTGCATTAAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                      <  2:1199886/100‑1 (MQ=255)
                              tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                     <  2:2178555/100‑1 (MQ=255)
                               cTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                    <  1:2412370/100‑1 (MQ=255)
                               cTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                    <  1:2033767/100‑1 (MQ=255)
                                 aTGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCTGTTTTTGTTTTACTTTTCagag                                                                  <  2:2116086/100‑1 (MQ=255)
                                 aTGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGATTCTGCCGTTTTTGTTTTACTTTTCagag                                                                  <  1:2945650/100‑1 (MQ=255)
                                  tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                 >  2:846451/1‑100 (MQ=255)
                                  tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                 <  2:1855067/100‑1 (MQ=255)
                                  tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                 >  1:2681807/1‑100 (MQ=255)
                                    tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                               <  2:2636633/100‑1 (MQ=255)
                                    tgCATTAAAAAAATCTTCATATAGATATCATTCTCCTGAACAATGCGATTAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                               <  2:1613210/100‑1 (MQ=255)
                                      aaaaaaaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                 <  2:368263/91‑1 (MQ=255)
                                       aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                            <  1:1561628/100‑1 (MQ=255)
                                       aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                            >  2:2762106/1‑100 (MQ=255)
                                       aTTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGACGTTTTTGTTTTACTTTTCAGCGGCgtg                                                             >  2:3147098/1‑99 (MQ=255)
                                        ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                           <  1:3597743/100‑1 (MQ=255)
                                        ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                           >  2:2235286/1‑100 (MQ=255)
                                          gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggtg                                                          >  1:544568/1‑99 (MQ=255)
                                          gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                         >  1:1158888/1‑100 (MQ=255)
                                          gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                         >  1:871741/1‑100 (MQ=255)
                                           aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                        >  1:2588852/1‑100 (MQ=255)
                                            aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                       <  1:1910278/100‑1 (MQ=255)
                                            aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                       <  2:3341813/100‑1 (MQ=255)
                                             aaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatga                                                      <  2:1170749/100‑1 (MQ=255)
                                              aaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaa                                                     <  1:1890971/100‑1 (MQ=255)
                                              aaaTCTTCAGATAGATATCATTATGCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCCGAGGAgtggtggtgca                                                     >  2:1376649/1‑98 (MQ=255)
                                                     cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                              <  1:1405908/100‑1 (MQ=255)
                                                      agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                             <  1:872596/100‑1 (MQ=255)
                                                       gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                            >  1:970867/1‑100 (MQ=255)
                                                       gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                            >  1:796577/1‑100 (MQ=255)
                                                       gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                            >  1:3147979/1‑100 (MQ=255)
                                                          agataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGt                                         <  1:2814325/100‑1 (MQ=255)
                                                               tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                    <  1:1309362/100‑1 (MQ=255)
                                                                cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                   <  1:1177988/100‑1 (MQ=255)
                                                                 aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                  <  2:1890979/100‑1 (MQ=255)
                                                                   tctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGcc                                <  2:1763861/100‑1 (MQ=255)
                                                                   tctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGcc                                <  1:2471334/100‑1 (MQ=255)
                                                                      ccTGAACAGTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGTAAAATTTTGAAGTTCTGCAGCCACt                             >  2:503059/1‑100 (MQ=255)
                                                                      ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                             >  1:1448390/1‑100 (MQ=255)
                                                                      ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                             >  1:2382078/1‑100 (MQ=255)
                                                                      ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                             >  1:1448352/1‑100 (MQ=255)
                                                                      ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                             <  1:2503243/100‑1 (MQ=255)
                                                                          aaCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAg                         >  1:1486838/1‑100 (MQ=255)
                                                                           aCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGa                        <  1:1397180/100‑1 (MQ=255)
                                                                           aCAATCCGATGAACGGCAGAACCTAGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGa                        <  2:1657736/100‑1 (MQ=255)
                                                                            cAATCCGATGAACGGCAGAACCGCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGTGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGaa                       >  2:2117547/1‑100 (MQ=255)
                                                                                ccGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTc                    >  2:223207/1‑99 (MQ=255)
                                                                                 cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                  <  1:2990183/100‑1 (MQ=255)
                                                                                  caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTATTCAGCCACTACAGAACTc                    >  2:3308489/5‑97 (MQ=37)
                                                                                    tggaCGGCAGAACCTCGCGGTTCTGACGTTTTTGTTGTACTTTTCATAGGAGTTGTGATGATAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt               <  2:704808/97‑1 (MQ=255)
                                                                                     gAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTc              <  1:1277351/100‑1 (MQ=255)
                                                                                       aCGGCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTTATAGGAGTGGTGATGAAAAATTTTTATTTTTTTCAGTCACTTCAGAACTCCCttt               >  2:3463071/1‑97 (MQ=16)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg           <  2:1202439/100‑1 (MQ=255)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg           <  1:523984/100‑1 (MQ=255)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg           >  1:2986548/1‑100 (MQ=255)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTTCTTTTTAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg           >  2:2638354/1‑100 (MQ=39)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGCGTGGTGATGAAAAACTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg            >  2:3120521/1‑99 (MQ=38)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg           >  1:3405950/1‑100 (MQ=255)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGGGGTGATGAAAAATTTTGACGTTCTGCAGCCACTACAGAACTCCCTTTCCgg           >  2:571607/1‑100 (MQ=39)
                                                                                         ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt          >  1:1296647/1‑100 (MQ=255)
                                                                                         ggCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTCAGAGGAGTGGTGATGAAAAATTTTGTTTTTCTGCAGCCACTACAGAACTCCCTTTCCGGt          >  2:97205/1‑100 (MQ=255)
                                                                                         ggCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTTTTTTCAGAGGAGTGGTGATGAAAAATTTTGATGTTTTTCAGCCACTAAAGAACTTCCTTTCAGGt          >  2:828089/1‑100 (MQ=255)
                                                                                         ggCAGAACCGCGAGGTTCTACCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt          >  1:1334183/1‑100 (MQ=255)
                                                                                          gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc         >  2:2659910/1‑100 (MQ=255)
                                                                                          gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc         >  2:403372/1‑100 (MQ=255)
                                                                                          gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTCCCGCCACTACAGAACTCCCTTTCCg            >  2:3601408/1‑97 (MQ=255)
                                                                                           cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt        <  1:2027201/100‑1 (MQ=255)
                                                                                           cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt        <  2:1453209/100‑1 (MQ=255)
                                                                                              aaCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTAcc     <  1:3553593/100‑1 (MQ=255)
                                                                                                ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc   >  2:2321054/1‑100 (MQ=255)
                                                                                                 cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAGCAATTTTGACGTTCTTCAGACACTACAGAACTACATTTCCGGTCTACCGc   >  2:2927352/1‑99 (MQ=255)
                                                                                                 cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCt  >  1:3522851/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCT  >  NC_000913/2767569‑2767765

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: