Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,154,675 T→A 5.1% I182F (ATC→TTC)  yqhC ← transcriptional activator of yqhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,154,6750TA 5.1% 242.7 / 2.7 78I182F (ATC→TTC) yqhCtranscriptional activator of yqhD
Reads supporting (aligned to +/- strand):  ref base T (29/44);  new base A (2/2);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

GCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCTCATCTTCGCCAATGTC  >  NC_000913/3154528‑3154822
                                                                                                                                                   |                                                                                                                                                   
gCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGgatgat                                                                                                                                                     <  7:322419/148‑1 (MQ=255)
 cTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCAGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCt                                                                                                                                                   >  8:363845/1‑149 (MQ=255)
  tCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTg                                                                                                                                                  >  2:3855/1‑149 (MQ=255)
        cTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGcc                                                                                                                                            <  2:162773/149‑1 (MQ=255)
           aGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCga                                                                                                                                          <  3:140679/148‑1 (MQ=255)
             gTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGaa                                                                                                                                       <  7:199534/149‑1 (MQ=255)
                ttCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGCAGTGAGGCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATa                                                                                                                                     <  3:209529/148‑1 (MQ=255)
                ttCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATa                                                                                                                                     <  5:137695/148‑1 (MQ=255)
                ttCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATa                                                                                                                                     <  7:188521/148‑1 (MQ=255)
                   ggTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATAcgc                                                                                                                                  <  8:186567/148‑1 (MQ=255)
                    gtgtATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATAcgcg                                                                                                                                 >  7:301083/1‑148 (MQ=255)
                      gtATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATAc                                                                                                                                    <  6:304832/143‑1 (MQ=255)
                      gtATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATAc                                                                                                                                    >  5:304832/1‑143 (MQ=255)
                       tatttattCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATc                                                                                                                             >  3:211856/1‑149 (MQ=255)
                        atttatcCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATc                                                                                                                             >  8:105939/1‑148 (MQ=255)
                         tttattCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATc                                                                                                                             >  4:212145/1‑147 (MQ=255)
                         tttattCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCa                                                                                                                           <  7:332538/149‑1 (MQ=255)
                          ttattCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGACAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCAGAGAATACGCGCATCCAg                                                                                                                          >  1:257574/1‑149 (MQ=255)
                           tattCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGt                                                                                                                         >  8:187908/1‑149 (MQ=255)
                            attCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGt                                                                                                                         <  5:244827/148‑1 (MQ=255)
                            attCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTg                                                                                                                        >  4:92818/1‑149 (MQ=255)
                               ctcAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTc                                                                                                                     >  2:136188/1‑149 (MQ=255)
                                tcAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCg                                                                                                                    >  8:197244/1‑149 (MQ=255)
                                 cAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGaa                                                                                                                                       <  4:179536/129‑1 (MQ=255)
                                 cAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGaa                                                                                                                                       >  3:179536/1‑129 (MQ=255)
                                  aaTCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGc                                                                                                                   <  1:216119/148‑1 (MQ=255)
                                  aaTCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGc                                                                                                                   <  7:122697/148‑1 (MQ=255)
                                  aaTCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCt                                                                                                                  >  4:331135/1‑149 (MQ=255)
                                  aaTCCGTTCCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTc                                                                                                                     <  2:80361/146‑1 (MQ=255)
                                    tCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTcc                                                                                                                <  3:170518/149‑1 (MQ=255)
                                    tCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTc                                                                                                                 >  1:295937/1‑148 (MQ=255)
                                       gTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCAt                                                                                                              >  3:42421/1‑148 (MQ=255)
                                       gTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATc                                                                                                             <  3:211821/149‑1 (MQ=255)
                                       gTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATc                                                                                                             <  2:24801/149‑1 (MQ=255)
                                           cAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACg                                                                                                          <  2:257574/148‑1 (MQ=255)
                                            aGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGCCCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTc                                                                                                        >  4:37427/1‑149 (MQ=255)
                                            aGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTc                                                                                                        >  1:181804/1‑149 (MQ=255)
                                            aGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTc                                                                                                        <  7:116647/149‑1 (MQ=255)
                                              cACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGa                                                                                                      <  4:42421/149‑1 (MQ=255)
                                                  cgGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTaa                                                                                                  <  2:25676/149‑1 (MQ=255)
                                                   ggCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAAc                                                                                                 <  4:211856/149‑1 (MQ=255)
                                                        aTCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGc                                                                                              <  2:181804/147‑1 (MQ=255)
                                                         tCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTcc                                                                                           >  7:152081/1‑149 (MQ=255)
                                                         tCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTcc                                                                                           <  1:3855/149‑1 (MQ=255)
                                                         tCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACATAGAACAGAATTTCGCGGATGATCTGTTTGCCGAGAATCCTCGCATACAGTGGTCGCTCCATCACGACGAGTAACCGCTcc                                                                                           >  8:336650/1‑149 (MQ=255)
                                                              cTGAAGGGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGGCCGGGCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCcgccgc                                                                                      <  5:333356/149‑1 (MQ=255)
                                                                   ctGAGTCTGGCGACTGACCAGCGCAAGTAACGCGCCGCCGAAAGGTACGGACAGCACGTAGTACAGAAATTCGCGGAAGAGCTGGTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCAt                                                                                  >  2:228196/2‑148 (MQ=255)
                                                                      aGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATaaaa                                                                              >  3:162982/1‑149 (MQ=255)
                                                                       gTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAAt                                                                             <  1:70851/149‑1 (MQ=255)
                                                                          tGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTc                                                                                                                 >  4:268576/1‑110 (MQ=255)
                                                                          tGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTc                                                                                                                 <  3:268576/110‑1 (MQ=255)
                                                                          tGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAAtctc                                                                          >  4:274568/1‑149 (MQ=255)
                                                                            gCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCtctt                                                                        <  7:146928/149‑1 (MQ=255)
                                                                            gCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCtctt                                                                        <  8:364590/149‑1 (MQ=255)
                                                                            gCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCtctt                                                                        <  1:136188/149‑1 (MQ=255)
                                                                             cGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCtcttc                                                                       <  2:261577/149‑1 (MQ=255)
                                                                              gACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCtcttct                                                                      >  2:22898/1‑149 (MQ=255)
                                                                                   aCCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGAt                                                                   <  1:22898/147‑1 (MQ=255)
                                                                                          ccAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCg                                                          <  7:105939/149‑1 (MQ=255)
                                                                                             gTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAg                                                       >  4:176377/1‑149 (MQ=255)
                                                                                                atgcgcCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGtt                                                     >  2:351979/3‑148 (MQ=255)
                                                                                                 cgcgCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTaa                                                   <  7:197244/149‑1 (MQ=255)
                                                                                                 cgcgCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTaa                                                   <  2:295937/149‑1 (MQ=255)
                                                                                                   cgccgccgcAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATc                                                 <  3:331135/149‑1 (MQ=255)
                                                                                                    gccgccgcAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATcc                                                <  4:162982/149‑1 (MQ=255)
                                                                                                    gccgccgcAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATcc                                                <  3:92818/149‑1 (MQ=255)
                                                                                                           cAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCCTCTGATAACGTGGCGGAGTTAATCCCgctggc                                         <  3:37427/149‑1 (MQ=255)
                                                                                                                    gTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGa                                >  1:171836/1‑149 (MQ=255)
                                                                                                                       aGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACgg                             <  5:152715/149‑1 (MQ=255)
                                                                                                                        gCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGc                            >  2:86118/1‑149 (MQ=255)
                                                                                                                         cACGTAGTACAGAATTTCGCGGATGAACTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCt                           <  1:164887/149‑1 (MQ=255)
                                                                                                                         cACGTAGTACAGAATTTCGCGGATGAACTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCt                           <  1:170313/149‑1 (MQ=255)
                                                                                                                         cACGTAGTACAGAATTTCGCGGATGAACTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCt                           >  2:164887/1‑149 (MQ=255)
                                                                                                                         cACGTAGTACAGAATTTCGCGGATGAACTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCt                           >  2:170313/1‑149 (MQ=255)
                                                                                                                          aCGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTg                          <  3:176377/149‑1 (MQ=255)
                                                                                                                                  cAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGc                   <  5:336871/148‑1 (MQ=255)
                                                                                                                                  cAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCt                  <  4:199774/149‑1 (MQ=255)
                                                                                                                                  cAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCt                  >  6:239079/1‑149 (MQ=255)
                                                                                                                                   aGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCTc                 <  5:283835/149‑1 (MQ=255)
                                                                                                                                               atgatCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCTCATCTTCGCCAAt     <  3:274568/149‑1 (MQ=255)
                                                                                                                                                  atCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAAc                                                                <  3:334914/87‑1 (MQ=255)
                                                                                                                                                  atCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAAc                                                                >  4:334914/1‑87 (MQ=255)
                                                                                                                                                  atCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCTCATCTTCGCCAATGTc  >  6:177865/1‑149 (MQ=255)
                                                                                                                                                   tCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAACATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCTCATCTTCGCCAATGTc  >  5:96522/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGGCGACTGACCAGCGCCAGTAACGCGCCGCCGCAAGGTCCGGTCAGCACGTAGTACAGAATTTCGCGGATGATCTGTTTGCCGAGAATACGCGCATCCAGTGGTCGCTCCATCACGTCGAGTAACCGCTCCGCCGCGCATAAAATCTCTTCTGATAACGTGGCGGAGTTAATCCCGCTGGCTGCCATCGACGGCTGGAAATGCTCATCTTCGCCAATGTC  >  NC_000913/3154528‑3154822

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