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.

CGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAT  >  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)
     aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                        <  1:1104868/100‑1 (MQ=255)
     aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                        <  1:2103501/100‑1 (MQ=255)
       gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttt                                                                                      <  1:123857/100‑1 (MQ=255)
        tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttt                                                                                     <  1:327596/100‑1 (MQ=255)
        tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttt                                                                                     <  2:1727355/100‑1 (MQ=255)
        tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttt                                                                                     <  1:1409471/100‑1 (MQ=255)
         cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttt                                                                                    <  1:2341504/100‑1 (MQ=255)
          acttGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                   <  1:742081/97‑1 (MQ=255)
          aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                   >  2:1268142/1‑100 (MQ=255)
           ccTGACGTAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttttt                                                                                   <  1:159576/99‑1 (MQ=255)
            cTGCCGGAAACCGATCGTCGCCGCGCCAACCGCGCACCGACAGACGATTAATTTCGTTTTTTAGTCCTAACAGCTCTGCCATCGCTGGCAGG‑tctttt                                                                                   <  2:1645086/98‑1 (MQ=255)
            cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttat                                                                                 <  1:34171/100‑3 (MQ=255)
            cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttttt                                                                                   <  1:2145851/98‑1 (MQ=255)
             tGACGGAAAACGATCGTCGCCGCGCCAACCTCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttttt                                                                                 <  2:560832/99‑1 (MQ=255)
             tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatg                                                                                <  1:1727543/100‑4 (MQ=38)
             tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCATG‑ttttttttt                                                                                >  2:474311/1‑100 (MQ=255)
             tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                <  2:628595/100‑1 (MQ=255)
             tGACGAAAAACGATCATCGAAGCGACAACCACGCACCGAAAGACAATTAATTTCGTGTTTTAGACATAAAAAACCTGCCAGCACTGGCAGG‑ttttttt                                                                                  >  2:229833/1‑98 (MQ=255)
              gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCTAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                               <  2:2212295/100‑1 (MQ=255)
               aCGGAAAACGATCGTCGACGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTCTTTTGc                                                                              <  2:1744330/100‑1 (MQ=255)
                cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAATTTCTTGCCAGCGATGGCAGG‑TTTTTTTTTGc                                                                              <  2:865356/99‑1 (MQ=255)
                cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                             >  2:1219178/1‑100 (MQ=255)
                cGGAAAACGATCGTCGCCGCGACAAACGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                >  2:500952/1‑97 (MQ=255)
                 ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                             <  2:881977/99‑1 (MQ=255)
                  gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTg                                                                           >  1:1862717/1‑100 (MQ=255)
                    aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGACTGTc                                                                         >  1:784484/1‑100 (MQ=255)
                     aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCg                                                                        <  1:2384514/100‑1 (MQ=255)
                      aCGATCTTCGTCGTGCTAACCGCGCACCGAAAGACGATTAGTTTCGTTTTTTAGTCATAAAAATCCTGGCAGCGATGGCAGG‑TTTTTTTTTGc                                                                              <  2:419466/93‑1 (MQ=255)
                       cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTg                                                                      >  1:85177/1‑100 (MQ=255)
                       cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTg                                                                      >  1:516463/1‑100 (MQ=255)
                       cGATCGTCCCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTg                                                                      <  2:1981760/100‑1 (MQ=255)
                        gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGa                                                                     >  1:893494/1‑100 (MQ=255)
                         aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTTTCGTGAg                                                                    >  1:775416/1‑100 (MQ=255)
                         aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAg                                                                    <  1:449930/100‑1 (MQ=255)
                               gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGca                                                                >  2:778904/1‑98 (MQ=255)
                                ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCacac                                                             <  1:878992/100‑1 (MQ=255)
                                        aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGATCGCACACCAATGCGa                                                     >  1:2365409/1‑100 (MQ=255)
                                         ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                     <  1:1847038/99‑1 (MQ=255)
                                         ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                    <  1:2013139/100‑1 (MQ=255)
                                          cgcgGACCGAAAGACGATTAATTTCGTTTTTTAGTCATGAAAAACCTGCCGGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                    <  2:1761080/99‑1 (MQ=255)
                                           gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                    <  1:194818/98‑1 (MQ=255)
                                             gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCTATGTGAGTAtc                                                <  2:1252625/100‑1 (MQ=255)
                                                 cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                            <  2:1080359/100‑1 (MQ=255)
                                                  gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                            <  1:1058188/100‑1 (MQ=255)
                                                  gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                           <  1:1714969/100‑1 (MQ=255)
                                                  gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                           >  1:2118636/1‑100 (MQ=255)
                                                    aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                         >  1:196289/1‑100 (MQ=255)
                                                      gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                       >  1:2246089/1‑100 (MQ=255)
                                                       aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                      <  2:633825/100‑1 (MQ=255)
                                                       aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                      <  2:2233911/100‑1 (MQ=255)
                                                        cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                     <  2:775407/100‑1 (MQ=255)
                                                          aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                   >  1:2338160/1‑100 (MQ=255)
                                                          aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                   <  1:985860/100‑1 (MQ=255)
                                                           ttATTTTTGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                         <  1:13011/93‑1 (MQ=255)
                                                             aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                >  2:1951111/1‑100 (MQ=255)
                                                              aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                               <  1:2149283/100‑1 (MQ=255)
                                                               tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                              <  2:12109/100‑1 (MQ=255)
                                                               tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                              <  2:2212198/100‑1 (MQ=255)
                                                               tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                              <  2:164760/100‑1 (MQ=255)
                                                                ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTTTTTTGAGCGCACACCAATGCGTGTATCTCTTTTATAAGTCtgt                                >  2:1167701/1‑97 (MQ=255)
                                                                 tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAACTCTGTGGtt                            >  2:1829524/1‑100 (MQ=255)
                                                                   gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                          <  1:566052/100‑1 (MQ=255)
                                                                      ttttAGTGATAGAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAAGGCGAGTATCTACCTTATAAGTATGTGGTTACGTc                       >  2:1255027/1‑100 (MQ=255)
                                                                      ttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                       <  2:1580655/100‑1 (MQ=255)
                                                                      ttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                       <  2:1502376/100‑1 (MQ=255)
                                                                       tttAGTCATAAAAAACCTGCCAGCGAGGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                      <  2:39886/100‑1 (MQ=255)
                                                                         tAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                    <  2:1168565/100‑1 (MQ=255)
                                                                           gTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGATTATCCCCCTTATAAAGCTGTGGTTACGTCAATa                   >  1:350533/1‑99 (MQ=255)
                                                                           gTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                   >  1:2267683/1‑99 (MQ=255)
                                                                             cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                >  2:337890/1‑100 (MQ=255)
                                                                              aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga               <  1:657961/100‑1 (MQ=255)
                                                                              aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga               <  1:2291129/100‑1 (MQ=255)
                                                                                aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc             >  1:1668102/1‑100 (MQ=255)
                                                                                 aaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt            >  1:1054003/1‑100 (MQ=255)
                                                                                  aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTTTTTTTTTTGAACGCACACAAATCCGAGTATCTCAAATATAAGTTTGTGTTTAAGGCAATaga                 >  2:1738674/1‑94 (MQ=9)
                                                                                  aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTTTCGTGAGCGCACACCAATGGGAGTATCTCCCTTTTAAGTCTGTGGTTACGTCAATag                  >  2:1597351/1‑93 (MQ=37)
                                                                                  aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           >  2:260495/1‑100 (MQ=255)
                                                                                  aaaaCCTGCCAGTGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           >  1:901430/1‑100 (MQ=255)
                                                                                  aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCCCACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGGCAATAGAGAGCtt           >  1:1724692/1‑100 (MQ=255)
                                                                                  aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           <  2:1887675/100‑1 (MQ=255)
                                                                                  aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt           >  1:1310266/1‑100 (MQ=255)
                                                                                   aaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTTTCTGTCTTTAGAGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTAAGTAGAGAACttt          >  2:2075930/1‑100 (MQ=25)
                                                                                    atccTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  1:1288137/98‑1 (MQ=255)
                                                                                    aaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  1:1109698/100‑1 (MQ=255)
                                                                                    aaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTCATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg         <  2:1008442/100‑1 (MQ=255)
                                                                                     aCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt        >  1:102551/1‑100 (MQ=255)
                                                                                      ccTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc       <  2:2094641/100‑1 (MQ=255)
                                                                                      ccTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc       <  1:1268155/100‑1 (MQ=255)
                                                                                        tGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac     >  1:63360/1‑100 (MQ=255)
                                                                                         gCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca    <  1:406759/100‑1 (MQ=255)
                                                                                         gCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca    >  2:712848/1‑100 (MQ=255)
                                                                                           cAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt  <  2:2271524/100‑1 (MQ=255)
                                                                                           cAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt  <  2:2338132/100‑1 (MQ=255)
                                                                                             gCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa   <  2:1497342/97‑1 (MQ=255)
                                                                                                |                                                                                               
CGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAT  >  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: ‑