Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,292,472 G→A 5.2% intergenic (+107/+37) galU → / ← hns glucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,292,4720GA 5.2% 568.8 / 9.5 170intergenic (+107/+37)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base G (85/75);  new base A (5/4);  total (90/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 173 total aligned reads displayed.

GATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGAAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGGCTATTCAGCAGTTCGTTCGGGTCAATACCGTCAGCGATCAG  >  NC_000913/1292200‑1292730
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                   
                                                                                                                     caagtaaTGTAAAGCCTGGGTTTAAGAAGAGATGGGCCTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGGGCAATCTACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGc                                                                                                                                              <  1:34523/270‑1 (MQ=255)
                                                                                                                               aaaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTg                                                                                                                                                                                                                                                                  <  5:529740/149‑1 (MQ=255)
                                                                                                                                     tGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGa                                                                                                                                                                                                                                                            >  8:873506/1‑149 (MQ=255)
                                                                                                                                      ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAt                                                                                                                                                                                                                                                           <  8:127609/149‑1 (MQ=255)
                                                                                                                                      ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                                                                                                                                           <  8:883040/149‑1 (MQ=255)
                                                                                                                                      ggCTTGAAGAAGAGATGGGCATAAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGCTTATTGCTTGATCAGGAAATCGTCGCGGGATTTACCTTGCTCATCCATTGCTTTTTTTATTCACGCTTGAGGACGGAGTTGGCCATTCCCGGTTTTCGtt                                                                                                                             >  1:49187/1‑275 (MQ=255)
                                                                                                                                         ttGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAtttt                                                                                                                                                                                                                                                        <  5:1191954/149‑1 (MQ=255)
                                                                                                                                           gaagagGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCGCGAAAAGGCGTTGAGCAATCGCCGCAGTTTTTTTATAGCTTATTCTTAAAAAATTGTCTTAACCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                                                                                                                                      <  6:514811/149‑1 (MQ=255)
                                                                                                                                            aacaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGATTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGCCGGGCTTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTACAGTGATTTACCTTGCCCATCTATTGCTTTTTTTATTACAGCTGGAGTACGGNCTTTGNNCGACCTAGTTTTTGttttcctg                                                                                                                       >  3:16471/4‑270 (MQ=255)
                                                                                                                                              gtagagATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                                                                                                                                           >  8:1000462/3‑141 (MQ=255)
                                                                                                                                              gaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                                                                                                                                           <  7:1000462/141‑1 (MQ=255)
                                                                                                                                              gaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                                                                                                                                           >  5:446038/1‑141 (MQ=255)
                                                                                                                                                  agaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                                                                                                                                      <  7:135343/142‑1 (MQ=255)
                                                                                                                                                   gaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGc                                                                                                                                                                                                                                              >  7:12864/1‑149 (MQ=255)
                                                                                                                                                      gggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAAt                                                                                                                                                                                                                                           >  8:570058/1‑149 (MQ=255)
                                                                                                                                                      gggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAAt                                                                                                                                                                                                                                           <  8:871693/149‑1 (MQ=255)
                                                                                                                                                          aTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTg                                                                                                                                                                                                                                               >  8:236921/1‑141 (MQ=255)
                                                                                                                                                          aTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTg                                                                                                                                                                                                                                               <  7:236921/141‑1 (MQ=255)
                                                                                                                                                          aTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTAc                                                                                                                                                                                                                                       >  6:820700/1‑149 (MQ=255)
                                                                                                                                                          aTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTAc                                                                                                                                                                                                                                       >  5:44022/1‑149 (MQ=255)
                                                                                                                                                              agaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                                                                                                                                                   <  6:1036618/149‑1 (MQ=255)
                                                                                                                                                                aagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                                                                                                                                                   >  5:1021516/1‑147 (MQ=255)
                                                                                                                                                                aagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGa                                                                                                                                                                                                                                 >  5:405819/1‑149 (MQ=255)
                                                                                                                                                                  gTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGatt                                                                                                                                                                                                                               <  7:873506/149‑1 (MQ=255)
                                                                                                                                                                        tCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                                                                                                                                                                                                         >  7:391434/1‑149 (MQ=255)
                                                                                                                                                                           gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGa                                                                                                                                                                                                                      <  7:796668/149‑1 (MQ=255)
                                                                                                                                                                            tATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                                                                                                                                                     >  8:1379095/1‑149 (MQ=255)
                                                                                                                                                                             aTCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                                                                                                                                                                                                                    <  5:175465/149‑1 (MQ=255)
                                                                                                                                                                                   gtTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaaa                                                                                                                                                                                                              <  7:1056421/149‑1 (MQ=255)
                                                                                                                                                                                    ttATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAttt                                                                                                                                                                                                                                                         >  4:331692/1‑105 (MQ=255)
                                                                                                                                                                                     tATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtc                                                                                                                                                                                                            <  8:391434/149‑1 (MQ=255)
                                                                                                                                                                                       tCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgt                                                                                                                                                                                                          <  5:820700/149‑1 (MQ=255)
                                                                                                                                                                                         cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                                                                                                                                                                                                        >  5:219454/1‑149 (MQ=255)
                                                                                                                                                                                          ggttaaCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                                                                                                                                       <  6:308217/145‑1 (MQ=255)
                                                                                                                                                                                          aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAAt                                                                                                                                                                                                                                           >  7:81778/1‑113 (MQ=255)
                                                                                                                                                                                          aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                                                                                                                                       >  7:1179704/1‑149 (MQ=255)
                                                                                                                                                                                          aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                                                                                                                                       <  8:1072692/149‑1 (MQ=255)
                                                                                                                                                                                          aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                                                                                                                                       >  7:1452150/1‑149 (MQ=255)
                                                                                                                                                                                            gAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAgg                                                                                                                                                                                                     <  5:1254286/149‑1 (MQ=255)
                                                                                                                                                                                                 ggCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAttt                                                                                                                                                                                                >  5:743278/1‑149 (MQ=255)
                                                                                                                                                                                                 ggCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAttt                                                                                                                                                                                                <  5:995810/149‑1 (MQ=255)
                                                                                                                                                                                                  gCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCTGACAATAAAAAATCCCGCCAGCGGCGGGa                                                                                                                                                                                                                                                            >  5:1193400/1‑88 (MQ=255)
                                                                                                                                                                                                    gTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAcc                                                                                                                                                                                             >  5:155325/1‑149 (MQ=255)
                                                                                                                                                                                                       gAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTg                                                                                                                                                                                          <  6:809463/149‑1 (MQ=255)
                                                                                                                                                                                                        aGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                                                                                                                                                         <  5:1256899/149‑1 (MQ=255)
                                                                                                                                                                                                             tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                                                                                                                                                    <  6:44022/149‑1 (MQ=255)
                                                                                                                                                                                                             tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                                                                                                                                                    >  7:534992/1‑149 (MQ=255)
                                                                                                                                                                                                             tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCTGTCCAGGTGTTATTTTCGGCGTTTTCGTCAACGTAGATATCTTTTGCCAGACGCTGAACACGTTTAGCTTTGGTGCCAGCTTTAAc                                                      >  2:240861/1‑275 (MQ=255)
                                                                                                                                                                                                             tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTggg                                                                                  >  4:239577/1‑245 (MQ=255)
                                                                                                                                                                                                                aCGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCAt                                                                                                                                                                                 >  5:1237997/1‑149 (MQ=255)
                                                                                                                                                                                                                    cGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCt                                                                                                                                                                             >  6:1015620/1‑149 (MQ=255)
                                                                                                                                                                                                                           ttATAGCTTATTCTTATTCAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCTGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTTCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCaa                                                                                                             <  2:329638/206‑1 (MQ=255)
                                                                                                                                                                                                                            tataGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAt                                                                                                                                                                     >  8:702379/1‑149 (MQ=255)
                                                                                                                                                                                                                              taGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAt                                                                                                                                                                     >  8:390974/1‑147 (MQ=255)
                                                                                                                                                                                                                              taGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAt                                                                                                                                                                     <  7:390974/147‑1 (MQ=255)
                                                                                                                                                                                                                              taGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTa                                                                                                                                                                   <  6:1466896/149‑1 (MQ=255)
                                                                                                                                                                                                                                    attcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTg                                                                                                                                                             <  7:1379095/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAg                                                                                                                                                          >  8:935129/1‑149 (MQ=255)
                                                                                                                                                                                                                                           ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                                                                                                                                      >  5:37723/1‑149 (MQ=255)
                                                                                                                                                                                                                                           ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGATCCCGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGGCTATTCAGCAGttcgttcg                        >  4:110902/1‑274 (MQ=255)
                                                                                                                                                                                                                                             aaaTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGc                                                                                                                                                    <  6:155325/149‑1 (MQ=255)
                                                                                                                                                                                                                                                ttGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCAtt                                                                                                                                                                                >  7:913148/1‑118 (MQ=255)
                                                                                                                                                                                                                                                ttGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCtt                                                                                                                                                 <  8:81733/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTg                                                                                                                                                <  5:1380382/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAg                                                                                                                                                                                                                                                 >  4:23001/1‑51 (MQ=255)
                                                                                                                                                                                                                                                   tCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                                                                                                                                                         <  5:659194/106‑1 (MQ=255)
                                                                                                                                                                                                                                                   tCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                                                                                                                                                         >  6:659194/1‑106 (MQ=255)
                                                                                                                                                                                                                                                     ttAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa                                                                                                                                                                                                   >  7:732165/1‑94 (MQ=255)
                                                                                                                                                                                                                                                     ttAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCa                                                                                                                                            >  5:841818/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      tAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCa                                                                                                                                                                                                                                             <  2:23001/51‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTc                                                                                                                                         <  6:1081869/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTc                                                                                                                                         <  5:1160124/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          ccGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCa                                                                                                                                       <  8:756785/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           cGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg                                                                                                                                      <  8:730512/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAgg                                                                                                                                     >  5:757380/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCCGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGGCTATTCAGCAGTTCGTTCGGGTCAATACCGTCAGCGAt     >  2:103816/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCt                                                                                                                                                  >  4:163915/1‑135 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGt                                                                                                                                    <  6:118008/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGt                                                                                                                              >  4:80466/1‑155 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCg                                                                                                                 >  2:340668/1‑168 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGTCTATTCAGCAGTTCGTTCGGGTAAATACCGTCAGCGAt     >  2:298015/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGGCTTTTTAGCAGTTCGTTCGTGGCAATCACGTCAGCGaa     >  4:302286/1‑274 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGCCCAGGTTTAATTTTCCCTGTTTTCGTCAAGGTAACACTATTTGGCGGCGCGCTGAGCCGGTTTAGCTTTGTTGACGGATTAACGGGCAGAACGGCATACGAGCCGTTCGTGAGGGCCCATACCTTAAGCGac     >  4:53989/1‑274 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATCACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGACGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGGCTATTAAGCAGTTCGTTCGGGTCAATACCGTCAGCGAt     >  2:144997/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGACCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCCCCGTTTTCGTCAACGTAGCTATCTTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGTTTTAACGGCAGCAAGGCTATTCAGCAGTTCGTTCGGGGCAAAACCGTTAACGAt     >  4:145144/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 gcttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTACTGCTTGATCAGGAAATCGACGAGGGATTTACCTTGCACATCCATTGCTTTTTTGATTACAGCTGGAGTATGGCCTTGGCCAGTCCAGGTTTTa                                                                                                                                >  6:1453809/5‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                                                                                                                                >  8:1252435/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                                                                                                                                <  8:1353476/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTCGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                                                                                                                                <  8:1373944/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGGCTATTCAGCAGTTCGTTCGGGTCAATACCGTCAGCGATCAg  >  2:126181/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                   ttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                                                                                                                                <  8:1185064/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                                                                                                                                >  7:1185064/3‑147 (MQ=255)
                                                                                                                                                                                                                                                                   ttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGt                                                                                                                              >  7:712643/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaaaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGa                                                                                                                                                                                                                                 >  6:1304638/1‑50 (MQ=255)
                                                                                                                                                                                                                                                                      aaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTc                                                                                                                           <  8:527162/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTc                                                                                                                           <  6:219454/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCg                                                                                                                          >  7:1397394/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        aTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGc                                                                                                                         >  7:100305/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            cgccgcTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGtttt                                                                                                                                 <  6:747348/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                            cgccgcTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGtttt                                                                                                                                 >  5:747348/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                   
GATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGAAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTCAACGTAGCTATATTTTGCCGGACGCTGAGCACGTTTAGCTTTGGTGCCAGATTTAACGGCAGCAAGGCTATTCAGCAGTTCGTTCGGGTCAATACCGTCAGCGATCAG  >  NC_000913/1292200‑1292730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: