Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,280 (T)9→8 100% intergenic (+50/‑167) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2720T.92.5% 512.4 / 34.3 133intergenic (+42/‑175)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base . (68/55);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCA  >  NC_000913/4060173‑4060366
                                                                                                    |                                                                                              
aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTTATTTATTACCTGCAATCGCTGGCAGGt                                                                                                <  2:674957/100‑1 (MQ=255)
aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                                <  2:1953897/100‑1 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  1:441895/1‑100 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  1:1162137/1‑100 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  2:341018/1‑100 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  2:1249188/1‑100 (MQ=255)
  gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                               <  1:1852899/99‑1 (MQ=255)
    cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                            <  1:1994523/100‑1 (MQ=255)
     aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                           <  2:2818476/100‑1 (MQ=255)
     aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                             >  2:56559/1‑98 (MQ=255)
      ccTGACGGAAAACGTTCGTTGTCGCGCTTACCGTGCAACGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                           <  1:171345/99‑1 (MQ=25)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAATACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                          <  2:1549760/100‑1 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                          >  2:2904830/1‑100 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                          <  2:2207932/100‑1 (MQ=255)
       cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGTACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                          <  1:344061/99‑2 (MQ=38)
       cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                         <  1:3695125/100‑1 (MQ=255)
       cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                         >  1:2827237/1‑100 (MQ=255)
        tGACGGAAAACTTTCTTCGCTGCGCTAATCTCGCACCGAAAGACGATTAATTTCGTTTTTTATTCATAAAAAACTTGTTATCGCTGGCAGGtttttta                                                                                          <  2:271046/98‑2 (MQ=255)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCTCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                        <  1:180461/100‑4 (MQ=38)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                        <  1:208462/100‑4 (MQ=255)
        tGACGGAAAACGATCGTCGCCGCGCCAACCACGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                        <  2:2646438/100‑2 (MQ=255)
         gACTGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                        <  2:655983/100‑4 (MQ=37)
         gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                       >  1:2212335/1‑98 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgac                                                                                      <  2:546263/100‑4 (MQ=255)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgact                                                                                     <  2:3768296/100‑7 (MQ=37)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcat                                                                                     >  1:1638558/1‑96 (MQ=255)
            ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                                   >  1:3379183/1‑100 (MQ=255)
            ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                                   >  1:902431/1‑100 (MQ=255)
              aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACTATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTc                                                                                 <  1:2402528/100‑1 (MQ=255)
              aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTACTCATAAAAAACCTGCCAGCCACGGCAGTctttttttcc                                                                                       >  2:3757962/1‑93 (MQ=255)
                aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGt                                                                               <  1:3528242/100‑1 (MQ=255)
                aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTACCAGCGATGGCAGGTTTTTTTTGCCTGTCGt                                                                               >  1:431629/1‑100 (MQ=255)
                   gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAg                                                                            <  2:1699289/100‑1 (MQ=255)
                    cgatctgccgcGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                           <  1:2995110/94‑1 (MQ=255)
                     tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCGGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                            <  1:2753996/99‑1 (MQ=255)
                      cgtcgCCGCGCCAACCGCGCACTGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                         <  1:1485521/100‑1 (MQ=255)
                       gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTAAAAGCAATGGCAGGAATTTTTTGCCTGTCGTGAGCGca                                                                        >  1:2196104/1‑100 (MQ=255)
                        tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcac                                                                       <  1:2257241/100‑1 (MQ=255)
                        tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAGACCTGCCAGCGATCGCACGTTTTTTTTGCCTGTCGTGAGCGcac                                                                       <  2:3461371/100‑1 (MQ=255)
                         cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcaca                                                                      <  1:1103681/100‑1 (MQ=255)
                          gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcagatcg                                                                     >  2:2995083/1‑94 (MQ=255)
                          gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                                     >  2:829087/1‑100 (MQ=255)
                             gcgcCAACCGCGCAACGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                                  >  2:927589/1‑100 (MQ=255)
                                ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCATGTAGTGAGCGCACACCAATGc                                                               >  1:3114709/1‑100 (MQ=255)
                                 cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                              >  2:58965/1‑100 (MQ=255)
                                 cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                              >  1:2096032/1‑100 (MQ=255)
                                  aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                             >  2:1824498/1‑100 (MQ=255)
                                  aaCCGCGCACCGAAAGACGATTAAAATCATTTTTTAGTGATAAAAAACATGACAGCGATGACAGGTTATTTTAACTAGTAGTGAACGCACACCAATgaga                                                             >  2:1016195/1‑100 (MQ=255)
                                   aCCGCGTACCGAAAGACGATTAATTTCGTTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                            <  2:755609/100‑1 (MQ=255)
                                   aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                            <  1:3450391/100‑1 (MQ=255)
                                   aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                            >  2:42099/1‑100 (MQ=255)
                                     cgtgtaCCTAAAGACGATTAATTTCGTTTTTTGGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGGCTGTCGTTAGCGCACACCAATGCGAGTa                                                          <  2:1641151/95‑1 (MQ=255)
                                       cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                        <  2:2028937/100‑1 (MQ=255)
                                       cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                        <  2:73504/100‑1 (MQ=255)
                                         cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                      >  2:3100297/1‑100 (MQ=255)
                                         cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                      >  2:2445727/1‑100 (MQ=255)
                                          aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                     >  1:599722/1‑100 (MQ=255)
                                           ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGGTTTTTTGTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                    <  1:3609338/100‑1 (MQ=255)
                                            cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                   <  2:3205365/100‑1 (MQ=255)
                                             ggaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                  >  1:101655/3‑100 (MQ=255)
                                             gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                  >  2:3807252/1‑100 (MQ=255)
                                                      ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                         <  1:1008336/100‑1 (MQ=255)
                                                        aaTTTCGTTTTTTAGTCATAAAAAACCTGTCAGCGATGGCAGGTTTTTTTT‑‑CCTGTCGTGAGCGCACACTAATGCGAGTATTTCCCTTATAAGTCTGTgg                                      <  2:1673726/100‑1 (MQ=255)
                                                         aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                      >  1:882177/1‑100 (MQ=255)
                                                         aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                      <  1:3445799/100‑1 (MQ=255)
                                                         aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                      <  1:2598641/100‑1 (MQ=255)
                                                          tttGGTTTTTTAGTCATAAAACACCTGGCAGCGAGGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACCCCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                     <  2:3706005/100‑1 (MQ=255)
                                                             cGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTCATAAGTCTGTGATTAc                                  >  2:3000141/1‑100 (MQ=255)
                                                              gTTTTTTAGTTATAAAAATCCTGCCAGTGATGGTAGGTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                  <  1:766726/99‑1 (MQ=255)
                                                              gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                 >  1:2800691/1‑100 (MQ=255)
                                                              gTTTTTTAGTCATAAAAAACCTGCCAGATATGTCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                 <  1:2287015/100‑1 (MQ=255)
                                                                      gTCATAAAAAACCTGCCAGCGAGGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAGCTCCCTTATAAGTCTGTGGTTACGTCAATag                          >  2:3898826/1‑99 (MQ=255)
                                                                        gcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                        >  2:2774968/5‑100 (MQ=37)
                                                                        cATAAAAAACCTGTCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAATTCTGTGGTTAAGTCAATata                         >  2:549861/1‑96 (MQ=255)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       <  2:2479935/100‑1 (MQ=255)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGACAATagaaa                       >  1:3680476/1‑98 (MQ=255)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAATATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       >  1:596486/1‑100 (MQ=255)
                                                                         caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCACTTATAAGTCTGTGGTTACGTCAATagaga                       >  1:1532326/4‑100 (MQ=37)
                                                                         aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      <  1:804374/100‑1 (MQ=255)
                                                                         aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  1:1888897/1‑100 (MQ=255)
                                                                          tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       >  1:2837469/1‑98 (MQ=255)
                                                                          aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  1:3344987/3‑100 (MQ=38)
                                                                          aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  1:2957881/3‑100 (MQ=38)
                                                                           aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                     >  1:2725461/2‑100 (MQ=38)
                                                                           aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTAACGTCAATAGAGAGc                     >  1:3875849/2‑100 (MQ=38)
                                                                            aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTCCTGTCATGAGCGCACACCCATGCAAGTATCTCTGTTATCAGTCTGTGGTTACGTAAATAGACAGCt                    >  2:3237417/1‑100 (MQ=21)
                                                                            aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                    >  1:2324585/1‑100 (MQ=38)
                                                                            aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTTTTTTCTTGAGCGCACAACAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                    >  2:191447/1‑100 (MQ=255)
                                                                            aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                    >  2:191522/1‑100 (MQ=255)
                                                                              aaaaCCTGCCATCGCTGGCAGGTTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                  >  1:1625109/1‑100 (MQ=38)
                                                                                aaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAATCTGTGGTTACGTCAATAGAGAGCTTTGt                >  1:532916/1‑100 (MQ=38)
                                                                                  ccTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca              >  1:1744835/1‑100 (MQ=255)
                                                                                   cTGCCAGCGATGGCAGGTTTTTTTTGCTTGTCGTGAGCGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac             >  2:2556234/1‑100 (MQ=255)
                                                                                    tGCCAGTGATGGCAGGTTTTTTTTGCCTGTCTTGAGCGCATACCAATGCGAGTATCTCCCTAATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca            <  1:583367/100‑1 (MQ=255)
                                                                                     gCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcataa           >  2:2163903/1‑97 (MQ=255)
                                                                                      ccAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt          >  2:3286128/1‑100 (MQ=255)
                                                                                       cAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt         >  2:1757853/1‑100 (MQ=255)
                                                                                         gCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAt       >  1:3748646/1‑100 (MQ=255)
                                                                                             tGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc   <  2:2202317/100‑1 (MQ=255)
                                                                                               gCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa  <  1:659932/99‑1 (MQ=255)
                                                                                                    |                                                                                              
ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCA  >  NC_000913/4060173‑4060366

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