Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 15.3% 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.3% 374.4 / 52.4 145intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (50/72);  new base T (6/16);  total (57/88)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

Only 100 of 148 total aligned reads displayed.

CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAT  >  NC_000913/4060162‑4060358
                                                                                                   |                                                                                                  
cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg                                                                                                    >  2:868870/1‑100 (MQ=255)
cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg                                                                                                    >  2:1438751/1‑100 (MQ=255)
 gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGc                                                                                                   <  1:1283039/100‑1 (MQ=255)
 gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGc                                                                                                   <  2:471836/100‑1 (MQ=255)
  tATTCGTAAACGTCACCTGACAGAAAACTATCGTCTCCGCTTCAACCGCGCACCGAAAAACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGc                                                                                                   <  2:2256007/99‑1 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTTATAAAAAACCTGCCAACGATgg                                                                                              >  2:2020706/1‑100 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                              >  1:12110/1‑100 (MQ=255)
        tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                            <  2:2348932/100‑1 (MQ=255)
         aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                           >  2:415792/1‑100 (MQ=255)
          aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGACATAAAAAACCTGCCAGCGATGGCAgg                                                                                          <  2:350255/100‑1 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                        <  1:1104868/100‑1 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                        <  1:2103501/100‑1 (MQ=255)
            cGTCACCTGACGGAAAACGATCGTCGCCGAGCCAACCGCGCACCGAAAGACGATGAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                       <  1:607808/100‑1 (MQ=39)
             gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttt                                                                                      <  1:123857/100‑1 (MQ=255)
              tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                     <  1:327596/100‑1 (MQ=255)
              tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                     <  1:1409471/100‑1 (MQ=255)
              tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                     <  2:1727355/100‑1 (MQ=255)
                 ccTGACGGAAAACGTTCTTCGCCGCGCCAACCGCGCACCTAAAGCCGATTTATTTCGTTTTTTAGTCATAAAAATCCTGCCATCGCTGGCAGGtttttt                                                                                   <  2:723555/99‑1 (MQ=37)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                  <  1:1303530/100‑2 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                  >  2:1115586/1‑100 (MQ=255)
                  cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                   <  1:2145851/98‑1 (MQ=255)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCTCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                 <  2:560832/99‑1 (MQ=255)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                <  1:1727543/100‑4 (MQ=38)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                <  2:2192832/100‑4 (MQ=38)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCATGttttttttt                                                                                >  2:474311/1‑100 (MQ=255)
                   tGACGAAAAACGATCATCGAAGCGACAACCACGCACCGAAAGACAATTAATTTCGTGTTTTAGACATAAAAAACCTGCCAGCACTGGCAGGttttttt                                                                                  >  2:229833/1‑98 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                               <  2:1622479/100‑1 (MQ=255)
                     aCGGAAAACGATCGTCGACGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTCTTTTGc                                                                              <  2:1744330/100‑1 (MQ=255)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                               >  1:343411/1‑98 (MQ=255)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                             >  2:1219178/1‑100 (MQ=255)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                             >  1:2220683/1‑100 (MQ=255)
                       ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                             <  2:881977/99‑1 (MQ=255)
                        gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTg                                                                           >  1:1862717/1‑100 (MQ=255)
                          aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGACTGTc                                                                         >  1:784484/1‑100 (MQ=255)
                            aCGATCTTCGTCGTGCTAACCGCGCACCGAAAGACGATTAGTTTCGTTTTTTAGTCATAAAAATCCTGGCAGCGATGGCAGGTTTTTTTTTGc                                                                              <  2:419466/93‑1 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                      >  1:85177/1‑100 (MQ=255)
                              gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGa                                                                     >  1:893494/1‑100 (MQ=255)
                               aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTTTCGTGAg                                                                    >  1:775416/1‑100 (MQ=255)
                               aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                    <  1:449930/100‑1 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGca                                                                >  2:778904/1‑98 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGaaca                                                              >  2:294742/1‑100 (MQ=255)
                                      ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCacac                                                             <  1:878992/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGATCGCACACCAATGCGa                                                     >  1:2365409/1‑100 (MQ=255)
                                               ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                     <  1:1847038/99‑1 (MQ=255)
                                               ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                    <  1:2013139/100‑1 (MQ=255)
                                                cgcgGACCGAAAGACGATTAATTTCGTTTTTTAGTCATGAAAAACCTGCCGGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                    <  2:1761080/99‑1 (MQ=255)
                                                 gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                    <  1:194818/98‑1 (MQ=255)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                  <  2:523848/99‑1 (MQ=255)
                                                   gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCTATGTGAGTAtc                                                <  2:1252625/100‑1 (MQ=255)
                                                       cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                            <  1:2135794/100‑1 (MQ=255)
                                                        gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                            <  1:1058188/100‑1 (MQ=255)
                                                        gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                            <  2:1362337/100‑1 (MQ=255)
                                                          aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                         >  1:196289/1‑100 (MQ=255)
                                                            gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                       >  1:2246089/1‑100 (MQ=255)
                                                             aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                      <  2:2233911/100‑1 (MQ=255)
                                                             aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                      <  2:803999/100‑1 (MQ=255)
                                                             aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                      <  2:633825/100‑1 (MQ=255)
                                                              cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                     >  1:67726/1‑100 (MQ=255)
                                                                aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                   >  1:2338160/1‑100 (MQ=255)
                                                                aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                   <  1:985860/100‑1 (MQ=255)
                                                                 ttATTTTTGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                         <  1:13011/93‑1 (MQ=255)
                                                                   aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                >  2:1951111/1‑100 (MQ=255)
                                                                     tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                               <  1:248652/99‑1 (MQ=255)
                                                                     tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                              <  2:12109/100‑1 (MQ=255)
                                                                     tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                              <  2:2212198/100‑1 (MQ=255)
                                                                     tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                              <  2:164760/100‑1 (MQ=255)
                                                                      ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTTTTTTGAGCGCACACCAATGCGTGTATCTCTTTTATAAGTCtgt                                >  2:1167701/1‑97 (MQ=255)
                                                                      ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                             >  1:1670590/1‑100 (MQ=255)
                                                                       tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAACTCTGTGGtt                            >  2:1829524/1‑100 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                          >  1:2355460/1‑100 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                          <  1:566052/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGACAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGATTATCTGCCTTATAAGTCTGTAGTTa                           >  2:2220937/1‑99 (MQ=255)
                                                                            ttttAGTGATAGAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAAGGCGAGTATCTACCTTATAAGTATGTGGTTACGTc                       >  2:1255027/1‑100 (MQ=255)
                                                                            ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                       <  1:468203/100‑1 (MQ=255)
                                                                            ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                       <  2:1502376/100‑1 (MQ=255)
                                                                             tttAGTCATAAAAAACCTGCCAGCGAGGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                      <  2:39886/100‑1 (MQ=255)
                                                                               tAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                    <  2:1168565/100‑1 (MQ=255)
                                                                                 gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGATTATCCCCCTTATAAAGCTGTGGTTACGTCAATa                   >  1:350533/1‑99 (MQ=255)
                                                                                 gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                   >  1:2267683/1‑99 (MQ=255)
                                                                                  caaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                 >  1:1430407/5‑100 (MQ=38)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                <  2:2313964/100‑1 (MQ=255)
                                                                                    aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga               <  1:657961/100‑1 (MQ=255)
                                                                                    aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga               <  1:2291129/100‑1 (MQ=255)
                                                                                    aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga               <  2:67725/100‑1 (MQ=255)
                                                                                     tAAAAAACCTGCCAGCGATGGCAGGTTTTTATTCTTTTTTTTTGAGAGCACACCAATGCGAGTATCTTCCTTATAAGTCTGTGGTTACGGCATTaga                 >  2:1577481/1‑97 (MQ=255)
                                                                                      aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc             >  1:1668102/1‑100 (MQ=255)
                                                                                       aaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt            >  1:1054003/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCATCGCTGGCAGGTTTTTTTTTTTTTTTTTTGAACGCACACAAATCCGAGTATCTCAAATATAAGTTTGTGTTTAAGGCAATaga                 >  2:1738674/1‑94 (MQ=9)
                                                                                        aaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTTTCGTGAGCGCACACCAATGGGAGTATCTCCCTTTTAAGTCTGTGGTTACGTCAATag                  >  2:1597351/1‑93 (MQ=37)
                                                                                        aaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           >  2:260495/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCAGTGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           >  1:901430/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCCCACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGGCAATAGAGAGCtt           >  1:1724692/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           <  2:1887675/100‑1 (MQ=255)
                                                                                         aaaCCTGCCATCGCTGGCAGGTTTTTTTTTTTCTGTCTTTAGAGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTAAGTAGAGAACttt          >  2:2075930/1‑100 (MQ=25)
                                                                                          aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  1:1109698/100‑1 (MQ=255)
                                                                                          aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTCATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  2:1008442/100‑1 (MQ=255)
                                                                                           aCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt        >  1:102551/1‑100 (MQ=255)
                                                                                            ccTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc       <  1:1268155/100‑1 (MQ=255)
                                                                                               gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca    >  2:712848/1‑100 (MQ=255)
                                                                                                 cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt  <  2:2354075/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAT  >  NC_000913/4060162‑4060358

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