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)
                         aaaaCGATCGTCGCCGCTCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTc                                                                           <  1:700679/100‑1 (MQ=255)
                         aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                            >  2:2644997/1‑100 (MQ=255)
                          aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCg                                                                          >  2:1292341/1‑100 (MQ=255)
                          aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCg                                                                          <  2:1614209/100‑1 (MQ=255)
                          aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCg                                                                          >  2:1160708/1‑100 (MQ=255)
                           aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGt                                                                         >  2:1450046/1‑100 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                        <  2:1721632/100‑1 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGa                                                                       >  2:1932000/1‑100 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGa                                                                       >  1:1285206/1‑100 (MQ=255)
                              gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAg                                                                      <  2:2783230/100‑1 (MQ=255)
                               aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgc                                                                     <  2:116553/100‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGATGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgc                                                                     <  1:452544/99‑1 (MQ=255)
                                 cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgcgc                                                                   >  2:742894/1‑100 (MQ=255)
                                  gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TCTTTTTTGCCTGTTGTTAGCGca                                                                  >  1:502274/1‑100 (MQ=255)
                                    cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCAAGCGATGGCAGG‑TTTTTTTTGCCTGTTGTGATGGcac                                                                 >  2:2099871/1‑99 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCacac                                                               <  2:2217286/100‑1 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCacac                                                               <  1:1662510/100‑1 (MQ=255)
                                      ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACAcc                                                              <  2:2828840/100‑1 (MQ=255)
                                       cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCa                                                             >  2:1850076/1‑100 (MQ=255)
                                          nccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                          >  1:1964571/2‑100 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                        <  1:1235514/100‑1 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                        <  2:2686297/100‑1 (MQ=255)
                                             aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                       >  1:211807/1‑100 (MQ=255)
                                             aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTCAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCGCACGAATGCGa                                                       <  1:1522763/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      <  2:582375/100‑1 (MQ=255)
                                               ccGCGCACCTAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                       >  1:2743750/1‑98 (MQ=255)
                                                cgcgCACCGATAGACGATTAATTTCGTTTTTTAGTCATAAAGAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      <  1:1729628/99‑1 (MQ=255)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                  <  2:1133277/100‑1 (MQ=255)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   >  2:516731/1‑99 (MQ=255)
                                                     aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                               <  2:208448/100‑1 (MQ=255)
                                                      ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                              <  2:847279/100‑1 (MQ=255)
                                                        gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTT‑CCTGTCGTGCGCGAACACCAAAGCGAGAAtc                                                  >  1:1825900/1‑94 (MQ=255)
                                                               gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                     <  2:1786690/100‑1 (MQ=255)
                                                                   aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                 <  2:2720441/100‑1 (MQ=255)
                                                                    aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                <  2:1383471/100‑1 (MQ=255)
                                                                      ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                              >  2:478137/1‑100 (MQ=255)
                                                                      ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTCAGCGCACACGAATGCGAGTATCTCCCTTATAAAACTGTGGt                               >  1:1892120/1‑99 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCTAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACg                           <  1:586836/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                            >  1:809244/1‑99 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                           >  1:2611836/1‑100 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                           <  2:1005238/100‑1 (MQ=255)
                                                                           tttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         <  2:752900/100‑1 (MQ=255)
                                                                           tttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         <  1:2717231/100‑1 (MQ=255)
                                                                           tttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         <  2:1297779/100‑1 (MQ=255)
                                                                            ttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                        <  1:1794010/100‑1 (MQ=255)
                                                                            ttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                        <  2:2611847/100‑1 (MQ=255)
                                                                              ttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                      <  2:1476602/100‑1 (MQ=255)
                                                                               tAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                     >  1:2832264/1‑100 (MQ=255)
                                                                               caggcaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTATGACTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                     >  1:829250/8‑100 (MQ=21)
                                                                                 ggcaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATACGTCTGTGGTTACGTCAATaga                   >  2:352888/6‑100 (MQ=37)
                                                                                 ggcaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                   >  2:2178761/6‑100 (MQ=25)
                                                                                 gTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                   >  1:1044796/1‑100 (MQ=255)
                                                                                  gcaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATGCGTCTGTGGTTAAGTCAATagag                  >  1:228310/5‑100 (MQ=25)
                                                                                   caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  1:38021/4‑100 (MQ=38)
                                                                                   caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  1:1567045/4‑100 (MQ=38)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 <  2:1212154/100‑1 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  1:1469622/1‑100 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  2:2309409/1‑100 (MQ=255)
                                                                                    aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                >  1:1224106/1‑100 (MQ=255)
                                                                                     tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCATTAGAGAGc               >  1:199873/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGACGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAAAGCttt            >  1:2690981/1‑100 (MQ=38)
                                                                                         aaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg           >  2:456050/1‑100 (MQ=39)
                                                                                         aaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg           >  1:1436782/1‑100 (MQ=39)
                                                                                           aCCTGCCAGCGATGGCAGGTTTTTTTTT‑TTTTTTTTTAACGCACAGCAATGAGAGTATATCCCTTATAGGTCTGTGGTTACGTc                         >  1:220426/1‑84 (MQ=255)
                                                                                           aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGATTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt          <  1:1508846/99‑1 (MQ=255)
                                                                                           aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt          <  2:1224108/99‑1 (MQ=255)
                                                                                           aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGGCTGTGGTTACGTCAATAGAGAGCTTTg           <  2:827967/98‑1 (MQ=255)
                                                                                            ccTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca        <  2:643087/100‑1 (MQ=255)
                                                                                                  aGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa  >  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: ‑