Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,331 G→A 5.5% R302H (CGT→CAT)  ydgH → DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,3310GA 5.5% 289.5 / 4.4 91R302H (CGT→CAT) ydgHDUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (44/42);  new base A (2/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

AAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGC  >  NC_000913/1679184‑1679477
                                                                                                                                                   |                                                                                                                                                  
aaGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGt                                                                                                                                                   >  2:113138/1‑149 (MQ=255)
  gCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTgg                                                                                                                                                 >  4:107590/1‑149 (MQ=255)
    gACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGATACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGACAGTGGCAGGAACGTGGta                                                                                                                                               <  1:163798/149‑1 (MQ=255)
      ccGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGtaat                                                                                                                                             <  1:399584/149‑1 (MQ=255)
           gCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCt                                                                                                                                        <  4:95588/149‑1 (MQ=255)
             gatgatGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGa                                                                                                                                      <  3:180508/149‑1 (MQ=255)
               tgatgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGAcc                                                                                                                                    >  5:111838/1‑149 (MQ=255)
                 atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGt                                                                                                                                  <  1:130882/149‑1 (MQ=255)
                 atgGTCCCGTTCGACAGCATCAAATACTCTGGCAAATACGGCAACATGACCGATGTCTCCTATCAGGTTGCGAAACGTGCCGCTAAGAAAGGTGCAAAGTATTACCACATCACCCTACAGTGGCAGGAACGTGGTAATAACCTGaaagt                                                                                                                                  >  6:42585/1‑145 (MQ=255)
                    gTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCa                                                                                                                                <  8:331436/148‑1 (MQ=255)
                    gTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAg                                                                                                                               >  3:285517/1‑149 (MQ=255)
                    gTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGAAATAACCTGACCGTCa                                                                                                                                <  3:107590/148‑1 (MQ=255)
                       ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAAAGTGCCGCGAAGAAAGGTGCTAAGTATTACCAAATCACACGCCAGTGGCAGGAAAGTGGTAATAAACTGACCGTCAgct                                                                                                                             >  6:266071/1‑147 (MQ=255)
                        cGTTCGACAGCATCAAATTCTCTGGCAACTACTGCAACATGACAGAAGTCTCCTATCAGGTTCCGAAACGTGCCGCGAAGAAAGGTTCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGACAgcg                                                                                                                             >  5:81846/1‑147 (MQ=255)
                        cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcg                                                                                                                             >  5:81844/1‑147 (MQ=255)
                            cGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGAt                                                                                                                        <  7:342669/148‑1 (MQ=255)
                             gACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCt                                                                                                                      >  5:392032/1‑149 (MQ=255)
                             gACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCt                                                                                                                      <  6:111838/149‑1 (MQ=255)
                               cAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGt                                                                                                                    <  8:69622/149‑1 (MQ=255)
                                 gCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAAAGTGCCGCGAAGAAAGGTGCTAATTATTACACCATCACCCGCCAGTGGCAGGAACGTGGTAATAACATGACCGTCAGCGCAGATCTGtat                                                                                                                  <  1:121517/149‑1 (MQ=255)
                                    tCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGAcc                                                                                                                                    <  3:136199/128‑1 (MQ=255)
                                    tCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGAcc                                                                                                                                    >  4:136199/1‑128 (MQ=255)
                                       aaTTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaataa                                                                                                            <  5:367430/149‑1 (MQ=255)
                                        aTTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAAc                                                                                                           >  6:153753/1‑149 (MQ=255)
                                            tctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAtt                                                                                                       >  5:82609/1‑149 (MQ=255)
                                              tGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAt                                                                                                        <  6:392032/146‑1 (MQ=255)
                                                gCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCa                                                                                                     <  1:139742/147‑1 (MQ=255)
                                                 cAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt                                                                                                  <  8:151512/149‑1 (MQ=255)
                                                  aaCTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTa                                                                                                 >  7:380486/1‑149 (MQ=255)
                                                  aaCTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACATGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTa                                                                                                 >  6:214242/1‑149 (MQ=255)
                                                  aaCTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACATGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTa                                                                                                 >  6:219528/1‑149 (MQ=255)
                                                       cGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTc                                                                                              <  5:173950/147‑1 (MQ=255)
                                                        ggCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACCTCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCta                                                                                           >  2:207331/1‑149 (MQ=255)
                                                         gCAACATGTCCGTATGCTCGTATCCTGGTGCGATTCGTGCCGCGAAGAATGGTTCTTAGTATTACCACATCACCCGCCAGTGGCAGGAACGTTGTTATATCCTGTCATTCAGCTCAGATATGTATAAATAACTATTCAGTTACTCCtat                                                                                          <  6:274818/149‑1 (MQ=255)
                                                           aaCATGACCGAAGTCTCCTATCAGGTCGCAAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGGCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa                                                                                        <  2:279706/149‑1 (MQ=255)
                                                             cATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCtat                                                                                          <  1:5191/145‑1 (MQ=255)
                                                             cATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCtat                                                                                          >  2:5191/1‑145 (MQ=255)
                                                              aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCATATCTGTATAAATAACTATTCAGTTACTCCTATAAggc                                                                                     <  2:76982/149‑1 (MQ=255)
                                                              aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAggc                                                                                     >  6:37554/1‑149 (MQ=255)
                                                              aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAggc                                                                                     >  3:191712/1‑149 (MQ=255)
                                                              aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAggc                                                                                     <  7:15884/149‑1 (MQ=255)
                                                                gACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAggcgg                                                                                   >  4:297030/1‑149 (MQ=255)
                                                                  ccGATGTCTCCTATCAGGTTGCGAAACGTGCCGCGATGAAAGGTTCTAAGTATTACCACATAACACGCAAGTTGCAGGAACGTGGTTATAACCTGACCGGCATCGCAGATATGTATAAATATCTATTCAGTTACTCCTATAAGGCGGCt                                                                                 <  5:42585/149‑1 (MQ=255)
                                                                  ccGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCt                                                                                                                      >  1:327906/1‑112 (MQ=255)
                                                                  ccGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCt                                                                                                                      <  2:327906/112‑1 (MQ=255)
                                                                     aaGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACATGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGa                                                                              <  2:184680/149‑1 (MQ=255)
                                                                      aGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGAt                                                                             <  2:386716/149‑1 (MQ=255)
                                                                      aGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGAt                                                                             <  4:59195/149‑1 (MQ=255)
                                                                        tctcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTgg                                                                                                                                                 <  1:129711/79‑1 (MQ=255)
                                                                        tctcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTgg                                                                                                                                                 >  2:129711/1‑79 (MQ=255)
                                                                        tctcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgc                                                                           >  3:364964/1‑149 (MQ=255)
                                                                          tcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccg                                                                         >  5:159891/1‑149 (MQ=255)
                                                                          tcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccg                                                                         >  8:96750/1‑149 (MQ=255)
                                                                            cTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgcc                                                                       >  3:340789/1‑149 (MQ=255)
                                                                            cTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgc                                                                        <  5:177905/148‑1 (MQ=255)
                                                                            cTATCAGGTTGCGAAACGTGCAGCGAAGAAAGGTGCTAAGTATTACCACCTCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTCTAACTAACTATTAAGTTACTCCTATAAGGCGGCTTGATgccgcc                                                                       <  5:266071/149‑1 (MQ=255)
                                                                             tATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCCGGAACGTGGTAATAACCTGACCGTCAGCGCAGCTCTGGATAAATAACTATTCAGTTACTCCTATAAGGCGGGTTGATCCCGCCt                                                                      >  2:186131/1‑149 (MQ=255)
                                                                              aTCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGCCCGTCAGCGAAGATCTGTATAAATACCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCtt                                                                     >  6:358477/1‑149 (MQ=255)
                                                                              aTCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCtt                                                                     >  4:63859/1‑149 (MQ=255)
                                                                                 aGGTTGCGTTACGTGCCGCGTAGTAAGGTGCTAAGTATTTCCAGCTCAGCGGTCAGTGGCTGGAACGTTGTAATAACCTGACCGTCAGCGCAGATATGTATAAATAACAATTCAGTTACTCCTTTAAGGAGGCTTGATGCCGACTTTc                                                                   <  6:81846/148‑1 (MQ=255)
                                                                                 aGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                   <  6:81844/148‑1 (MQ=255)
                                                                                 aGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                   <  4:80451/148‑1 (MQ=255)
                                                                                  ggTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGCCCGTCCGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCg                                                                  >  5:12126/1‑148 (MQ=255)
                                                                                   gTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCa                                                                >  8:164909/1‑149 (MQ=255)
                                                                                   gTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCa                                                                >  6:136864/1‑149 (MQ=255)
                                                                                   gTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCa                                                                >  7:58417/1‑149 (MQ=255)
                                                                                       cGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAttt                                                             <  5:37554/148‑1 (MQ=255)
                                                                                        gAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATATCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAattt                                                            >  6:108203/1‑148 (MQ=255)
                                                                                         aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAtttt                                                            >  3:38332/1‑147 (MQ=255)
                                                                                            cGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTc                                                       <  2:289503/149‑1 (MQ=255)
                                                                                                   cgAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAtt                                                >  6:93992/1‑149 (MQ=255)
                                                                                                   cgAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAtt                                                >  2:169646/1‑149 (MQ=255)
                                                                                                          aGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCcgc                                         <  6:159891/149‑1 (MQ=255)
                                                                                                          aGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCcgc                                         <  6:82609/149‑1 (MQ=255)
                                                                                                             tGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGAc                                      <  1:113138/149‑1 (MQ=255)
                                                                                                                tAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCAt                                   >  2:77700/1‑149 (MQ=255)
                                                                                                                  aGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCAtt                                  <  3:297030/148‑1 (MQ=255)
                                                                                                                    tATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgca                               <  5:153753/149‑1 (MQ=255)
                                                                                                                        aCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgc                            <  6:233451/148‑1 (MQ=255)
                                                                                                                         ccACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAg                          >  4:203445/1‑149 (MQ=255)
                                                                                                                             aTCACCCGCCAGTGGCAGGAACATGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGacac                      <  5:219528/149‑1 (MQ=255)
                                                                                                                             aTCACCCGCCAGTGGCAGGAACATGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGacac                      <  5:214242/149‑1 (MQ=255)
                                                                                                                                aCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACAtc                   >  1:167971/1‑149 (MQ=255)
                                                                                                                                  ccGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACAtctc                 >  5:167752/1‑149 (MQ=255)
                                                                                                                                     ccAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTc              <  4:340789/149‑1 (MQ=255)
                                                                                                                                      cAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCg             <  8:380486/149‑1 (MQ=255)
                                                                                                                                            cAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAAt       >  3:145279/1‑149 (MQ=255)
                                                                                                                                             aGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATc      >  2:365720/1‑149 (MQ=255)
                                                                                                                                               gAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATccc    >  2:227917/1‑149 (MQ=255)
                                                                                                                                                 aCGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCcgc  >  3:244414/1‑149 (MQ=255)
                                                                                                                                                 aCGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCcgc  >  5:389058/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGC  >  NC_000913/1679184‑1679477

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑