Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 18.8% intergenic (+33/‑50) ygeI → / → pbl uncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,8890TA18.8% 357.1 / 75.0 154intergenic (+33/‑50)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (69/55);  new base A (15/14);  total (84/70)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

Only 100 of 159 total aligned reads displayed.

GCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATA  >  NC_000913/2993790‑2993986
                                                                                                   |                                                                                                 
gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                                                                                                   >  2:1535105/1‑100 (MQ=255)
gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                                                                                                   <  1:122304/100‑1 (MQ=255)
 ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                  <  1:1006772/100‑1 (MQ=255)
 ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                  <  1:1773383/100‑1 (MQ=255)
 ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                  <  1:1761706/100‑1 (MQ=255)
     ttCAAAATGAAAGTTATTCAGAAGATACTCCTATCTCACCTTACTATAAAATTAAGTTGTATACGATTTACCAAATAAATCATAACGAGATACTTTta                                                                                                <  2:2551913/98‑1 (MQ=255)
      tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataa                                                                                             <  1:1735053/100‑1 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                                                                                              >  2:1332400/1‑98 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                            >  2:1767205/1‑100 (MQ=255)
        aaaaTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataata                                                                                           >  1:231611/1‑100 (MQ=255)
          aaTGATATTTATTCAGAAGATAATCCTATCTCACCTTACTATATAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt                                                                                         <  2:1796221/100‑1 (MQ=255)
           aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTa                                                                                        <  1:952653/100‑1 (MQ=255)
           aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTa                                                                                        <  1:2005322/100‑1 (MQ=255)
                atGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCg                                                                                   <  1:511698/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGATTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                  <  2:344162/100‑1 (MQ=39)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                  >  2:443286/1‑100 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                  <  1:1449780/100‑1 (MQ=255)
                  gTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt                                                                                 <  2:3075214/100‑1 (MQ=255)
                    aTTCAGAATATAATCCTATCTCACATGATTATAAAATGGAATGGAATACTATTGACGAAAGAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                               <  2:376150/100‑1 (MQ=255)
                    aTTCAGAAGATAATGCTATCTCACCTTACTAGAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                               <  2:1634187/100‑1 (MQ=255)
                    aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                               <  2:2409769/100‑3 (MQ=255)
                     ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCt                                                                              >  1:438510/1‑100 (MQ=255)
                     ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                              <  1:3533159/100‑4 (MQ=38)
                     ttCAGAAGATAATCCTATCTCACCTTACTATAAAAATGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                               <  2:1820181/99‑1 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                            >  2:816384/1‑100 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                            <  2:320756/100‑1 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                            >  1:3268013/1‑100 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                            >  1:3257500/1‑100 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatgat                                                                            <  1:1016041/100‑6 (MQ=38)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                           >  1:2961495/1‑100 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                           <  1:1604479/100‑1 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatgatt                                                                           <  1:3150501/100‑7 (MQ=37)
                          cgatctaTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                         <  2:641513/94‑1 (MQ=255)
                          aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                         <  2:3375576/100‑1 (MQ=255)
                            gaTAATCCTATCTCAGCTTACTATTAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattc                                                                       <  2:1571802/100‑1 (MQ=255)
                            cttAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattc                                                                       <  2:2516798/98‑1 (MQ=255)
                             aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                      <  1:2653815/100‑1 (MQ=255)
                             aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                      <  1:1133548/100‑1 (MQ=255)
                              tAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattcag                                                                     >  1:2516814/1‑99 (MQ=255)
                              tAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCaa                                                                     >  2:251630/1‑100 (MQ=255)
                                aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                   >  2:3321094/1‑100 (MQ=255)
                                 tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                  <  1:2092892/100‑1 (MQ=255)
                                   cTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTc                                                                >  2:2303498/1‑100 (MQ=255)
                                      tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcaca                                                             <  1:1802040/100‑1 (MQ=255)
                                        tcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacac                                                            >  1:695451/1‑99 (MQ=255)
                                         cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACAt                                                          >  1:3324343/1‑100 (MQ=255)
                                           ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                        >  2:1725985/1‑100 (MQ=255)
                                            cTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                       <  2:2961472/100‑1 (MQ=255)
                                                cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaaa                                                   >  1:575163/1‑100 (MQ=255)
                                                cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaaa                                                   >  1:2250218/1‑100 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGATCAAATAAATCATAACGAGATACTTTGATAATAGTATCTCGTTCCTTTCATTCAATGGTCACACAGGaacaaa                                                   >  2:3081209/1‑98 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAat                                                  >  2:2170801/1‑99 (MQ=255)
                                                    aaaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATg                                               >  1:569827/1‑100 (MQ=255)
                                                    aaaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATg                                               <  1:401715/100‑1 (MQ=255)
                                                       aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcgc                                            >  1:3578670/1‑100 (MQ=255)
                                                            gTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAcc                                       >  2:3136333/1‑100 (MQ=255)
                                                              ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                     >  2:2153042/1‑100 (MQ=255)
                                                              ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                     >  1:505/1‑100 (MQ=255)
                                                              ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTAACCaa                                     >  1:1619143/1‑100 (MQ=255)
                                                               gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaa                                    >  2:642524/1‑100 (MQ=255)
                                                                tAGACGATTGACCAAATAAATAATAACAAGATACTTTTATAATAGTATCTCATTCCTTACATTCAAAACTCACACATGAAAAAATATGCGCTCACCcaaa                                   >  2:1261031/1‑100 (MQ=255)
                                                                 aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                  >  2:1631969/1‑100 (MQ=255)
                                                                 aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                  <  2:1016517/100‑1 (MQ=255)
                                                                  gACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGt                                 >  1:859280/1‑100 (MQ=255)
                                                                     gATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaa                              >  1:2601057/1‑100 (MQ=255)
                                                                     gATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaa                              >  2:2007810/1‑100 (MQ=255)
                                                                      aTTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaaa                             <  2:1828818/100‑1 (MQ=255)
                                                                         gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                          >  1:2396997/1‑100 (MQ=255)
                                                                         gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAAACCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                          >  2:2563062/1‑100 (MQ=255)
                                                                         gACCAAATAAATCATAACAAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTAACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                          >  2:2053745/1‑100 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                         <  1:103290/100‑1 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                         <  2:2348183/100‑1 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                         <  1:3346138/100‑1 (MQ=255)
                                                                           ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                        >  2:2044974/1‑100 (MQ=255)
                                                                            cAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtaa                       <  1:2451507/100‑3 (MQ=255)
                                                                            cAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                       <  2:3461312/100‑1 (MQ=255)
                                                                              gaatgaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                     >  2:2062043/11‑100 (MQ=25)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                     >  1:615524/1‑100 (MQ=255)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                     >  1:1673601/1‑100 (MQ=255)
                                                                                 tgaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCAt                  >  1:3364628/8‑100 (MQ=37)
                                                                                  aatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                 <  1:1418930/100‑1 (MQ=255)
                                                                                   atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCt                <  1:1729253/100‑1 (MQ=255)
                                                                                   atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCt                <  1:1378240/100‑1 (MQ=255)
                                                                                   aaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCt                >  2:2691922/6‑100 (MQ=38)
                                                                                   aaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCt                >  1:2166540/6‑100 (MQ=38)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa             >  1:265180/3‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc            >  2:443882/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc            >  1:1322382/2‑100 (MQ=255)
                                                                                       aaaCGAGATACTATTACAAAAGTATCTAGTTCCTTTCATTCAATCCTAACACATGACAAAATATGCGCTGCCCAAAAGTTAAATTGTTAACAATCttaa             >  2:3103915/2‑97 (MQ=255)
                                                                                        aaCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg           <  2:2736639/100‑1 (MQ=255)
                                                                                         aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt          >  2:1412392/1‑100 (MQ=255)
                                                                                         aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt          >  1:738716/1‑100 (MQ=255)
                                                                                          cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt         >  2:1584933/1‑100 (MQ=255)
                                                                                          cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt         >  2:1974191/1‑100 (MQ=255)
                                                                                           gagaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt        <  1:189503/100‑1 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg      >  1:2520382/1‑100 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg      <  1:2892194/100‑1 (MQ=255)
                                                                                               tACTTTTATAATACTGTCGCGTTCCTTTCATTCAATCCCCACACATGAAAAATTATTCTCTCACTAAAAGTTAAATTCTTAACCATATTAACgctct       >  2:2249977/1‑93 (MQ=255)
                                                                                               tACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta    <  1:2094721/100‑1 (MQ=255)
                                                                                                 cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCAGACATGAAAAAATATGCGCTCACCAAACGTTAAATTCTTAACCATCTTCACGTTTTGtat   >  1:791227/1‑99 (MQ=255)
                                                                                                   |                                                                                                 
GCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATA  >  NC_000913/2993790‑2993986

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: