Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 15.8% intergenic (+34/‑183) 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,2640AC15.8% 345.6 / 57.3 143intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (50/70);  new base C (7/16);  total (57/86)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

Only 100 of 143 total aligned reads displayed.

CGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAT  >  NC_000913/4060168‑4060358
                                                                                                |                                                                                               
cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                              >  1:12110/1‑100 (MQ=255)
 gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                             >  2:703868/1‑100 (MQ=255)
  tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                            <  2:2348932/100‑1 (MQ=255)
   aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                           >  2:415792/1‑100 (MQ=255)
    aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCATCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                          <  2:243380/100‑1 (MQ=255)
    aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGACATAAAAAACCTGCCAGCGATGGCAgg                                                                                          <  2:350255/100‑1 (MQ=255)
     aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                        <  1:1104868/100‑1 (MQ=255)
     aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                        <  1:2103501/100‑1 (MQ=255)
       gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttt                                                                                      <  1:123857/100‑1 (MQ=255)
        tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                     <  1:327596/100‑1 (MQ=255)
        tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                     <  2:1727355/100‑1 (MQ=255)
        tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                     <  1:1409471/100‑1 (MQ=255)
         cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                    <  1:2341504/100‑1 (MQ=255)
          acttGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                   <  1:742081/97‑1 (MQ=255)
          aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                   >  2:1268142/1‑100 (MQ=255)
           ccTGACGTAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                   <  1:159576/99‑1 (MQ=255)
            cTGCCGGAAACCGATCGTCGCCGCGCCAACCGCGCACCGACAGACGATTAATTTCGTTTTTTAGTCCTAACAGCTCTGCCATCGCTGGCAGGtctttt                                                                                   <  2:1645086/98‑1 (MQ=255)
            cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                 <  1:34171/100‑3 (MQ=255)
            cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                   <  1:2145851/98‑1 (MQ=255)
             tGACGGAAAACGATCGTCGCCGCGCCAACCTCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                 <  2:560832/99‑1 (MQ=255)
             tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                <  1:1727543/100‑4 (MQ=38)
             tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCATGttttttttt                                                                                >  2:474311/1‑100 (MQ=255)
             tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                <  2:628595/100‑1 (MQ=255)
             tGACGAAAAACGATCATCGAAGCGACAACCACGCACCGAAAGACAATTAATTTCGTGTTTTAGACATAAAAAACCTGCCAGCACTGGCAGGttttttt                                                                                  >  2:229833/1‑98 (MQ=255)
              gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCTAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                               <  2:2212295/100‑1 (MQ=255)
               aCGGAAAACGATCGTCGACGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTCTTTTGc                                                                              <  2:1744330/100‑1 (MQ=255)
                cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAATTTCTTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                              <  2:865356/99‑1 (MQ=255)
                cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                             >  2:1219178/1‑100 (MQ=255)
                cGGAAAACGATCGTCGCCGCGACAAACGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                >  2:500952/1‑97 (MQ=255)
                 ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                             <  2:881977/99‑1 (MQ=255)
                  gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTg                                                                           >  1:1862717/1‑100 (MQ=255)
                    aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGACTGTc                                                                         >  1:784484/1‑100 (MQ=255)
                     aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCg                                                                        <  1:2384514/100‑1 (MQ=255)
                      aCGATCTTCGTCGTGCTAACCGCGCACCGAAAGACGATTAGTTTCGTTTTTTAGTCATAAAAATCCTGGCAGCGATGGCAGGTTTTTTTTTGc                                                                              <  2:419466/93‑1 (MQ=255)
                       cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                      >  1:85177/1‑100 (MQ=255)
                       cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                      >  1:516463/1‑100 (MQ=255)
                       cGATCGTCCCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                      <  2:1981760/100‑1 (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)
                                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)
                                             gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCTATGTGAGTAtc                                                <  2:1252625/100‑1 (MQ=255)
                                                 cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                            <  2:1080359/100‑1 (MQ=255)
                                                  gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                            <  1:1058188/100‑1 (MQ=255)
                                                  gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                           <  1:1714969/100‑1 (MQ=255)
                                                  gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                           >  1:2118636/1‑100 (MQ=255)
                                                    aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                         >  1:196289/1‑100 (MQ=255)
                                                      gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                       >  1:2246089/1‑100 (MQ=255)
                                                       aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                      <  2:633825/100‑1 (MQ=255)
                                                       aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                      <  2:2233911/100‑1 (MQ=255)
                                                        cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                     <  2:775407/100‑1 (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)
                                                              aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                               <  1:2149283/100‑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)
                                                                 tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAACTCTGTGGtt                            >  2:1829524/1‑100 (MQ=255)
                                                                   gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                          <  1:566052/100‑1 (MQ=255)
                                                                      ttttAGTGATAGAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAAGGCGAGTATCTACCTTATAAGTATGTGGTTACGTc                       >  2:1255027/1‑100 (MQ=255)
                                                                      ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                       <  2:1580655/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)
                                                                             cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                >  2:337890/1‑100 (MQ=255)
                                                                              aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga               <  1:657961/100‑1 (MQ=255)
                                                                              aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga               <  1:2291129/100‑1 (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)
                                                                                  aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           >  1:1310266/1‑100 (MQ=255)
                                                                                   aaaCCTGCCATCGCTGGCAGGTTTTTTTTTTTCTGTCTTTAGAGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTAAGTAGAGAACttt          >  2:2075930/1‑100 (MQ=25)
                                                                                    atccTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  1:1288137/98‑1 (MQ=255)
                                                                                    aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  1:1109698/100‑1 (MQ=255)
                                                                                    aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTCATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  2:1008442/100‑1 (MQ=255)
                                                                                     aCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt        >  1:102551/1‑100 (MQ=255)
                                                                                      ccTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc       <  2:2094641/100‑1 (MQ=255)
                                                                                      ccTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc       <  1:1268155/100‑1 (MQ=255)
                                                                                        tGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac     >  1:63360/1‑100 (MQ=255)
                                                                                         gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca    <  1:406759/100‑1 (MQ=255)
                                                                                         gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca    >  2:712848/1‑100 (MQ=255)
                                                                                           cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt  <  2:2271524/100‑1 (MQ=255)
                                                                                           cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt  <  2:2338132/100‑1 (MQ=255)
                                                                                             gCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa   <  2:1497342/97‑1 (MQ=255)
                                                                                                |                                                                                               
CGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAT  >  NC_000913/4060168‑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: