Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,772,236 A→T 15.5% pseudogene (831/900 nt) yfjV ← pseudogene, CP4‑57 prophage;Phage or Prophage Related; putative pump protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,772,2360AT15.5% 480.4 / 69.1 180pseudogene (831/900 nt)yfjVpseudogene, CP4‑57 prophage;Phage or Prophage Related; putative pump protein
Reads supporting (aligned to +/- strand):  ref base A (79/73);  new base T (15/13);  total (94/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 180 total aligned reads displayed.

TGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAT  >  NC_000913/2772139‑2772334
                                                                                                 |                                                                                                  
 ggTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATg                                                                                                 >  2:2780198/1‑100 (MQ=255)
 ggTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATg                                                                                                 <  2:2802757/100‑1 (MQ=255)
  gTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGc                                                                                                >  1:2949031/1‑100 (MQ=255)
  gTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGc                                                                                                <  2:3100577/100‑1 (MQ=255)
   tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGcc                                                                                               <  1:1432714/100‑1 (MQ=255)
   tAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGcc                                                                                               <  1:2571896/100‑1 (MQ=255)
      tGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGt                                                                                            <  2:806320/100‑1 (MQ=25)
      tGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGt                                                                                            >  1:806324/1‑100 (MQ=25)
      tGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAATCATGACAATGCCGGt                                                                                            >  2:2306220/1‑100 (MQ=255)
        ttGCGCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGACAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTa                                                                                           <  2:2973420/99‑1 (MQ=255)
         tGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgct                                                                                         <  1:2818992/100‑4 (MQ=25)
          gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctg                                                                                        <  1:3196259/100‑5 (MQ=25)
          gCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctg                                                                                        <  2:3681717/100‑5 (MQ=25)
           cTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGa                                                                                       >  1:2414690/1‑100 (MQ=255)
            tCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctgtt                                                                                      <  1:322478/100‑7 (MQ=37)
            tCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGaa                                                                                      >  2:17607/1‑100 (MQ=255)
            tCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGaa                                                                                      >  1:2648332/1‑100 (MQ=255)
             cATATCAGTATCTAATTACGAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTCACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAg                                                                                     <  2:659996/100‑1 (MQ=255)
              atatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACTCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctgtttg                                                                                    <  2:2196479/100‑9 (MQ=25)
               tatCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTa                                                                                   <  2:1513231/100‑1 (MQ=255)
                atCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCATTGTCATGACTGTGCCGGtgctgtttgtc                                                                                  <  2:3944283/100‑11 (MQ=25)
                 tCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAat                                                                                 >  2:598705/1‑100 (MQ=255)
                  cAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAata                                                                                >  1:838402/1‑100 (MQ=255)
                    gTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATc                                                                              <  2:1384723/100‑1 (MQ=255)
                      aTCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATccc                                                                            >  2:3357509/1‑100 (MQ=255)
                        cTAAGTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCa                                                                          >  2:177180/1‑100 (MQ=255)
                           aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGt                                                                       <  2:2954148/100‑1 (MQ=255)
                           aTTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGACGGTCAGGAAGTACTATCCCCAAGt                                                                       >  1:2701077/1‑100 (MQ=255)
                             tACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGa                                                                     >  1:908604/1‑100 (MQ=255)
                              aCAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAt                                                                    >  1:3808079/1‑100 (MQ=255)
                               cAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGAtt                                                                   >  1:1215936/1‑100 (MQ=255)
                                 aaGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGt                                                                 <  2:966260/100‑1 (MQ=255)
                                   gTGACAGAGAGCTGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCa                                                               <  2:2846552/100‑1 (MQ=255)
                                   gTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCa                                                               <  1:3859169/100‑1 (MQ=255)
                                   gTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCa                                                               <  1:1040753/100‑1 (MQ=255)
                                     gACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCata                                                             <  2:3735654/100‑1 (MQ=255)
                                       cAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtt                                                           >  1:3499171/1‑100 (MQ=255)
                                       cAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtt                                                           >  2:3408460/1‑100 (MQ=255)
                                        agagagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttt                                                          >  1:1581963/1‑100 (MQ=255)
                                         gagagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAtttt                                                         <  1:662720/100‑1 (MQ=255)
                                          agagCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATAttttt                                                        >  1:3650740/1‑100 (MQ=255)
                                             gccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGCTTGTCATATTTTTCTg                                                     >  2:2130835/1‑100 (MQ=255)
                                             gccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                     <  2:253840/100‑1 (MQ=255)
                                             gccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                     >  2:1491899/1‑100 (MQ=255)
                                             gccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTg                                                     >  1:476127/1‑100 (MQ=255)
                                             gccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGCAATATCCCCAAGTGATTGTCATATTTTTCTg                                                     >  2:1293032/1‑100 (MQ=255)
                                              ccgccACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTgg                                                    <  1:2257117/100‑1 (MQ=255)
                                               cgccacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGCAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                   >  1:3098241/1‑100 (MQ=255)
                                               cgccacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTggg                                                   <  1:2279242/100‑1 (MQ=255)
                                                gccacgccaGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGc                                                  <  1:2030697/100‑1 (MQ=255)
                                                       caGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGcaca                                           <  1:3946455/100‑1 (MQ=255)
                                                       caGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGcaca                                           >  2:2094842/1‑100 (MQ=255)
                                                       caGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGcaca                                           <  2:3808051/100‑1 (MQ=255)
                                                         gcgcGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                         <  2:767844/100‑1 (MQ=255)
                                                         gcgcGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATg                                         <  2:2553132/100‑1 (MQ=255)
                                                           gcgGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCa                                       <  1:1534944/100‑1 (MQ=255)
                                                            cgGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAg                                      >  1:1623618/1‑100 (MQ=255)
                                                            cgGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAg                                      >  2:1136032/1‑100 (MQ=255)
                                                              gCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGcc                                    >  2:106765/1‑100 (MQ=255)
                                                              gCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGcc                                    <  2:273570/100‑1 (MQ=255)
                                                                cAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCaagatc                                  >  1:3835313/1‑95 (MQ=255)
                                                                cAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcac                                  >  2:44014/1‑100 (MQ=255)
                                                                cAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcac                                  >  1:3659334/1‑100 (MQ=255)
                                                                cAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCcac                                  <  1:1520192/100‑1 (MQ=255)
                                                                  gagTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcag                                >  2:3200571/1‑100 (MQ=255)
                                                                    gTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagca                              >  1:2050313/1‑100 (MQ=255)
                                                                    gTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTAAGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagca                              >  1:3910582/1‑100 (MQ=255)
                                                                    gTGACAAACAGCACCGGCACAGTCATGACAATGCCGATACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCAcagca                              >  1:2781192/1‑100 (MQ=255)
                                                                      gACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCaga                            >  2:1522198/1‑100 (MQ=255)
                                                                           tccgtaCCGGCATTGTCATGACTGTGCCGGTGCGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGc                       >  1:1707127/6‑100 (MQ=25)
                                                                           tccgtaCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGc                       >  1:343037/6‑100 (MQ=37)
                                                                            cAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGcc                      <  1:2284093/100‑1 (MQ=255)
                                                                            cAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGcc                      >  1:232276/1‑100 (MQ=255)
                                                                             cgtaCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                     >  1:1786944/4‑100 (MQ=37)
                                                                             aGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCa                     <  2:3271997/100‑1 (MQ=255)
                                                                              gCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAg                    <  2:476125/100‑1 (MQ=255)
                                                                              gCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAg                    >  1:3057071/1‑100 (MQ=255)
                                                                              gCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAg                    <  1:17607/100‑1 (MQ=255)
                                                                               taCCGGCATTGTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGa                   >  2:330726/2‑100 (MQ=38)
                                                                                 ccggcattgTCATGACTGTGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTACCAGACt                 >  1:37080/9‑100 (MQ=38)
                                                                                  cGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTg                <  1:3339207/100‑1 (MQ=255)
                                                                                   ggCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCGGCCACAGCAGAGTTGCCAGACTGc               <  2:989041/100‑1 (MQ=255)
                                                                                   ggCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTGTCTGGGCAAGTACATGCAGCCACAGCAGTGTTGCCAGACttc               <  1:2661737/100‑3 (MQ=38)
                                                                                    gCACAGTCATTACAATGCCGGTACGCAAGTAATTTCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCACCCACAGCAGAGTTGCCATACTg                <  2:3536254/98‑1 (MQ=255)
                                                                                    gCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc              <  1:3448734/100‑1 (MQ=255)
                                                                                    gCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc              <  1:1519422/100‑1 (MQ=255)
                                                                                    gCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGcc              >  1:1161009/1‑100 (MQ=255)
                                                                                       caGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAAt           >  2:557625/1‑100 (MQ=255)
                                                                                        ngTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATc          >  1:1023733/2‑100 (MQ=255)
                                                                                         gTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCg         >  2:1603815/1‑100 (MQ=255)
                                                                                           cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCggg       >  1:689986/1‑100 (MQ=255)
                                                                                           cATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCgg        >  2:1818342/1‑99 (MQ=255)
                                                                                            acgACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCgggg      <  2:880133/98‑1 (MQ=255)
                                                                                            aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTTTTTTCTTTTACATGAAGCCACAGCAGAGTTGCCAGACTTCCAATCggg       >  2:3674211/1‑99 (MQ=255)
                                                                                            aTGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCgggg      <  2:3124800/100‑1 (MQ=255)
                                                                                             tGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCgggg      >  2:2175229/1‑99 (MQ=255)
                                                                                             tGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATACAGCCACAGCAGAGTTGCCAGACTGCCAATCggggc     >  2:2260423/1‑99 (MQ=255)
                                                                                              gACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTg    <  1:76814/100‑1 (MQ=255)
                                                                                              gACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTg    <  1:3033690/100‑1 (MQ=255)
                                                                                                cAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAt  <  1:2450307/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGCGGCCAGAGTGACAAACAGCACCGGCACAGTCATGACAATGCCGGTACGGAAGTAATATCCCCAAGTGATTGTCATATTTTTCTGGGCAAGCACATGCAGCCACAGCAGAGTTGCCAGACTGCCAATCGGGGTGAT  >  NC_000913/2772139‑2772334

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