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.91.4% 385.9 / 31.1 105intergenic (+42/‑175)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base . (49/47);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 132 total aligned reads displayed.

ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAG  >  NC_000913/4060173‑4060369
                                                                                                   |                                                                                                 
aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                                   <  2:1208505/100‑1 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                  >  1:16687/1‑100 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTTGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                             <  2:2054986/100‑1 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                             <  2:1781828/100‑2 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                             >  2:2779958/1‑100 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                             <  2:2721146/100‑1 (MQ=255)
        ttactgAAACCGATCATCGTCGCGACAACGGCGCACCGAAAGACTATTATTTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                            <  1:2478317/94‑3 (MQ=21)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                           <  2:2380963/100‑4 (MQ=255)
        ctaCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                           <  1:2363809/98‑4 (MQ=38)
         gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                          <  1:258724/100‑3 (MQ=255)
         gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                          <  1:1472104/100‑3 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCTCACCGAAACACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCATCGCTGGCAGGttttttatg                                                                                           <  1:1630034/98‑4 (MQ=37)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgac                                                                                         <  2:2500061/100‑6 (MQ=38)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                         <  2:331906/100‑4 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                         <  1:1197965/100‑4 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                         <  1:2475088/100‑4 (MQ=255)
           cGGAAAACTATCGTCGCCGCGCCAACCGCGCACCGCAAGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGttttttgtgcct                                                                                        <  1:1005218/100‑7 (MQ=255)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                                        >  1:491849/1‑96 (MQ=255)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                                        <  2:2704382/100‑5 (MQ=255)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                        <  1:1938515/100‑1 (MQ=255)
            ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCACAAAAAACCTGCAAGCGATGGCAGGttttttttgcc                                                                                         >  1:1955610/1‑95 (MQ=255)
             gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAGCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                                      >  2:2223124/1‑99 (MQ=255)
              aaaaCGATCGTCGCCGCTCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTc                                                                                    <  1:700679/100‑1 (MQ=255)
               aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                                   <  2:1614209/100‑1 (MQ=255)
               aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                                   >  1:2783220/1‑100 (MQ=255)
               aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                                   >  2:1292341/1‑100 (MQ=255)
               aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                                   >  2:1160708/1‑100 (MQ=255)
                aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGt                                                                                  >  2:1450046/1‑100 (MQ=255)
                  cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                                 <  2:1721632/100‑1 (MQ=255)
                  cGATCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGa                                                                                >  1:1285206/1‑100 (MQ=255)
                   gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAg                                                                               <  2:2783230/100‑1 (MQ=255)
                   gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAg                                                                               <  1:2203830/100‑1 (MQ=255)
                    aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                              >  2:1980974/1‑100 (MQ=255)
                     tcgtcgCCGCGCCAACCGCGCACCGAAAGATGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                              <  1:452544/99‑1 (MQ=255)
                      cgtcgCCGCGCCAGCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAGCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                            <  1:1276262/100‑1 (MQ=255)
                      cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                            >  2:742894/1‑100 (MQ=255)
                       gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGca                                                                           >  2:1283781/1‑100 (MQ=255)
                       gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTCTTTTTTGCCTGTTGTTAGCGca                                                                           >  1:502274/1‑100 (MQ=255)
                        tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcac                                                                          <  1:1172840/100‑1 (MQ=255)
                           ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                       <  2:2828840/100‑1 (MQ=255)
                            cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCa                                                                      >  2:1850076/1‑100 (MQ=255)
                               nccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                   >  1:1964571/2‑100 (MQ=255)
                                 cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                                 <  2:2686297/100‑1 (MQ=255)
                                 cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                                 <  2:627190/100‑1 (MQ=255)
                                 cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                                 <  1:1235514/100‑1 (MQ=255)
                                  aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                >  1:211807/1‑100 (MQ=255)
                                  aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTCAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCGCACGAATGCGa                                                                <  1:1522763/100‑1 (MQ=255)
                                   aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                               <  2:582375/100‑1 (MQ=255)
                                    ccGCGCACCTAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                >  1:2743750/1‑98 (MQ=255)
                                     cgcgCACCGATAGACGATTAATTTCGTTTTTTAGTCATAAAGAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                               <  1:1729628/99‑1 (MQ=255)
                                     cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                              <  1:2644986/100‑1 (MQ=255)
                                       cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                           <  2:1133277/100‑1 (MQ=255)
                                       cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                            >  2:516731/1‑99 (MQ=255)
                                          aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                        <  2:211809/100‑1 (MQ=255)
                                          aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                        <  2:208448/100‑1 (MQ=255)
                                           ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                       <  2:847279/100‑1 (MQ=255)
                                             gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTCCTGTCGTGCGCGAACACCAAAGCGAGAAtc                                                           >  1:1825900/1‑94 (MQ=255)
                                                    gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                              <  2:1786690/100‑1 (MQ=255)
                                                         aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                         <  2:1383471/100‑1 (MQ=255)
                                                          tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                          <  2:1073750/98‑1 (MQ=255)
                                                           ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                                       >  2:478137/1‑100 (MQ=255)
                                                           ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                                       >  2:701886/1‑100 (MQ=255)
                                                           ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTCAGCGCACACGAATGCGAGTATCTCCCTTATAAAACTGTGGt                                        >  1:1892120/1‑99 (MQ=255)
                                                              gTTTTTTAGTCATAAAAAACCTGCTAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACg                                    <  1:586836/100‑1 (MQ=255)
                                                              gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                    >  1:2611836/1‑100 (MQ=255)
                                                              gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                    <  2:1005238/100‑1 (MQ=255)
                                                                tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                  <  2:1297779/100‑1 (MQ=255)
                                                                tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                  <  2:1720205/100‑1 (MQ=255)
                                                                tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                  <  1:2717231/100‑1 (MQ=255)
                                                                tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                  <  2:752900/100‑1 (MQ=255)
                                                                 ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                                 <  2:2611847/100‑1 (MQ=255)
                                                                  tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                                >  2:836093/1‑100 (MQ=255)
                                                                   ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                               <  2:1476602/100‑1 (MQ=255)
                                                                    tAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                              <  1:2474021/100‑1 (MQ=255)
                                                                    tAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                              >  1:2832264/1‑100 (MQ=255)
                                                                    caggcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTATGACTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                              >  1:829250/8‑100 (MQ=21)
                                                                      ggcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATACGTCTGTGGTTACGTCAATaga                            >  2:352888/6‑100 (MQ=37)
                                                                      ggcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                            >  2:2178761/6‑100 (MQ=25)
                                                                      gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                            >  1:1044796/1‑100 (MQ=255)
                                                                        caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                          >  1:38021/4‑100 (MQ=38)
                                                                        caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                          >  1:1567045/4‑100 (MQ=38)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                          >  1:1469622/1‑100 (MQ=255)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                          >  2:2309409/1‑100 (MQ=255)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                          <  2:1212154/100‑1 (MQ=255)
                                                                          aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                        >  2:1203147/2‑100 (MQ=38)
                                                                          aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                        >  1:297440/2‑100 (MQ=38)
                                                                           aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                        >  1:1886362/1‑99 (MQ=255)
                                                                           aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                       >  1:655413/1‑100 (MQ=255)
                                                                              aaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                    >  2:456050/1‑100 (MQ=39)
                                                                              aaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                    >  1:1436782/1‑100 (MQ=39)
                                                                                aCCTGCCAGCGATGGCAGGTTTTTTTTTTTTTTTTTTAACGCACAGCAATGAGAGTATATCCCTTATAGGTCTGTGGTTACGTc                                  >  1:220426/1‑84 (MQ=255)
                                                                                aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                   <  2:1224108/99‑1 (MQ=255)
                                                                                 ccTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                  <  1:173664/100‑1 (MQ=255)
                                                                                 ccTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                 <  2:643087/100‑1 (MQ=255)
                                                                                     ccAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt             <  2:2670550/100‑1 (MQ=255)
                                                                                       aGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa           >  1:891012/1‑100 (MQ=255)
                                                                                           aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg       <  2:907899/100‑1 (MQ=255)
                                                                                            tGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc      <  2:2157132/100‑1 (MQ=255)
                                                                                             ggCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa     >  2:2222105/1‑100 (MQ=255)
                                                                                              gCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaa    >  2:2514381/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAG  >  NC_000913/4060173‑4060369

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: