Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 15.2% intergenic (+31/‑186) 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,2610GT15.2% 378.7 / 45.6 136intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (47/68);  new base T (13/8);  total (60/76)
Fisher's exact test for biased strand distribution p-value = 9.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 137 total aligned reads displayed.

CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4060162‑4060360
                                                                                                   |                                                                                                   
cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg                                                                                                     >  2:847854/1‑100 (MQ=255)
  tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCg                                                                                                   >  2:1662515/1‑100 (MQ=255)
   aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                  >  1:1305853/1‑100 (MQ=255)
     tCGTAAACGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                 <  1:2542280/99‑1 (MQ=255)
     gcGTAAACGTCACCTGACGGAAAACGATCGACGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                <  1:1290213/99‑1 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGATTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                               >  2:1662105/1‑100 (MQ=255)
       gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                              <  2:662531/100‑1 (MQ=255)
       gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                              <  2:606759/100‑1 (MQ=255)
          aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                           <  2:2290586/100‑1 (MQ=255)
          aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATCAATTTCGCTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                            >  1:1110108/1‑99 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                          <  2:1781375/100‑1 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                          <  2:1208505/100‑1 (MQ=255)
            cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                         >  1:16687/1‑100 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTTGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                    <  2:2054986/100‑1 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                    <  1:408898/100‑2 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                    >  2:2779958/1‑100 (MQ=255)
                  cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                   <  2:662063/100‑1 (MQ=255)
                   ttactgAAACCGATCATCGTCGCGACAACGGCGCACCGAAAGACTATTATTTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                   <  1:2478317/94‑3 (MQ=21)
                   ctaCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                  <  1:2363809/98‑4 (MQ=38)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAACAAACCTGCCAGCGATGGCAGGtttttttt                                                                                   >  1:1896586/1‑98 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                 <  1:258724/100‑3 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGAACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                 <  1:2817937/100‑3 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgag                                                                                >  2:2363820/1‑95 (MQ=38)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgac                                                                                <  2:2500061/100‑6 (MQ=38)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                <  1:1197965/100‑4 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                <  2:331906/100‑4 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                <  1:558017/100‑4 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                <  1:2475088/100‑4 (MQ=255)
                      cGGAAAACTATCGTCGCCGCGCCAACCGCGCACCGCAAGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGttttttgtgcct                                                                               <  1:1005218/100‑7 (MQ=255)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                               <  2:2704382/100‑5 (MQ=255)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                               <  1:1938515/100‑1 (MQ=255)
                         aaaaCGATCGTCGCCGCTCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTc                                                                           <  1:700679/100‑1 (MQ=255)
                         aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                            >  2:2644997/1‑100 (MQ=255)
                          aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                          >  2:1292341/1‑100 (MQ=255)
                          aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                          <  2:1614209/100‑1 (MQ=255)
                          aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                          >  2:1160708/1‑100 (MQ=255)
                           aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGt                                                                         >  2:1450046/1‑100 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                        <  2:1721632/100‑1 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGa                                                                       >  2:1932000/1‑100 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGa                                                                       >  1:1285206/1‑100 (MQ=255)
                              gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAg                                                                      <  2:2783230/100‑1 (MQ=255)
                               aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                     <  2:116553/100‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGATGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                     <  1:452544/99‑1 (MQ=255)
                                 cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                   >  2:742894/1‑100 (MQ=255)
                                  gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTCTTTTTTGCCTGTTGTTAGCGca                                                                  >  1:502274/1‑100 (MQ=255)
                                    cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCAAGCGATGGCAGGTTTTTTTTGCCTGTTGTGATGGcac                                                                 >  2:2099871/1‑99 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                               <  2:2217286/100‑1 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                               <  1:1662510/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                                                        <  1:1235514/100‑1 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                        <  2:2686297/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)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                  <  2:1133277/100‑1 (MQ=255)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   >  2:516731/1‑99 (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)
                                                                   aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                 <  2:2720441/100‑1 (MQ=255)
                                                                    aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                <  2:1383471/100‑1 (MQ=255)
                                                                      ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                              >  2:478137/1‑100 (MQ=255)
                                                                      ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTCAGCGCACACGAATGCGAGTATCTCCCTTATAAAACTGTGGt                               >  1:1892120/1‑99 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCTAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACg                           <  1:586836/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                            >  1:809244/1‑99 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                           >  1:2611836/1‑100 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                           <  2:1005238/100‑1 (MQ=255)
                                                                           tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         <  2:752900/100‑1 (MQ=255)
                                                                           tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         <  1:2717231/100‑1 (MQ=255)
                                                                           tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         <  2:1297779/100‑1 (MQ=255)
                                                                            ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                        <  1:1794010/100‑1 (MQ=255)
                                                                            ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                        <  2:2611847/100‑1 (MQ=255)
                                                                              ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                      <  2:1476602/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)
                                                                                  gcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATGCGTCTGTGGTTAAGTCAATagag                  >  1:228310/5‑100 (MQ=25)
                                                                                   caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  1:38021/4‑100 (MQ=38)
                                                                                   caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  1:1567045/4‑100 (MQ=38)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 <  2:1212154/100‑1 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  1:1469622/1‑100 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  2:2309409/1‑100 (MQ=255)
                                                                                    aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                >  1:1224106/1‑100 (MQ=255)
                                                                                     tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCATTAGAGAGc               >  1:199873/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGACGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAAAGCttt            >  1:2690981/1‑100 (MQ=38)
                                                                                         aaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg           >  2:456050/1‑100 (MQ=39)
                                                                                         aaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg           >  1:1436782/1‑100 (MQ=39)
                                                                                           aCCTGCCAGCGATGGCAGGTTTTTTTTTTTTTTTTTTAACGCACAGCAATGAGAGTATATCCCTTATAGGTCTGTGGTTACGTc                         >  1:220426/1‑84 (MQ=255)
                                                                                           aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGATTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt          <  1:1508846/99‑1 (MQ=255)
                                                                                           aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt          <  2:1224108/99‑1 (MQ=255)
                                                                                           aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGGCTGTGGTTACGTCAATAGAGAGCTTTg           <  2:827967/98‑1 (MQ=255)
                                                                                            ccTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca        <  2:643087/100‑1 (MQ=255)
                                                                                                  aGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa  >  1:891012/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4060162‑4060360

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: