Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,586,450 C→A 6.4% L240L (CTC→CTA narQ → sensory histidine kinase in two‑component regulatory system with NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,586,4500CA 6.4% 267.3 / 2.1 85L240L (CTC→CTAnarQsensory histidine kinase in two‑component regulatory system with NarP
Reads supporting (aligned to +/- strand):  ref base C (49/29);  new base A (2/4);  total (51/34)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGGCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCGGTGAA  >  NC_000913/2586305‑2586597
                                                                                                                                                 |                                                                                                                                                   
acGGGCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAACCCTTTAAGCAGATGTCGAGCGAGCTGCATAAATTGTAGCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCAc                                                                                                                                                  >  1:270583/1‑149 (MQ=255)
  gggCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCAc                                                                                                                                                  >  1:362166/1‑147 (MQ=255)
    gCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAgg                                                                                                                                              >  6:382909/1‑149 (MQ=255)
        ttCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCATGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCCTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGAACTCCACGAGGCCaa                                                                                                                                          >  2:204190/1‑149 (MQ=255)
           gACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGCGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGcg                                                                                                                                       >  1:350972/1‑149 (MQ=255)
           gACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGcg                                                                                                                                       >  4:295283/1‑149 (MQ=255)
                gccgccgcTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCgtcgtc                                                                                                                                  >  8:371494/1‑149 (MQ=255)
                    ccgcTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGa                                                                                                                              <  7:368897/149‑1 (MQ=255)
                          gATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGtgtt                                                                                                                        <  7:59298/149‑1 (MQ=255)
                                  ccTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTAAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCTATTTTATTTTAATTATAAGTATAGAACCGAGATCTACACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAgt                                                                                                                <  6:199664/149‑1 (MQ=255)
                                        gAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGTGCTGCATAAATTTTACCTTTCTTTTTTTTTTTAATTATAAGTATAGAACCGAGATCTTCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGc                                                                                                          <  6:296843/149‑1 (MQ=255)
                                        gAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCg                                                                                                           >  4:108959/1‑148 (MQ=255)
                                         aaTGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGCTGTTGTATCAGTGTTca                                                                                                           >  2:16270/1‑146 (MQ=255)
                                         aaTGAGCTTGGTCTGATTGCAAAAACCTTTAACAAGATGTCGAGCGAGCTCAATAAATTGTACCGTTCGCTGGAAGCTTCAGTAGAAGAAAAGACCCGCGATATCCACAAGGCAAATCGACGACTGGATGTGTTCTATCAGTCTTcca                                                                                                          >  2:9299/1‑146 (MQ=255)
                                          aTGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAg                                                                                                        >  3:224387/1‑149 (MQ=255)
                                             aGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGgcg                                                                                                     >  5:130780/1‑149 (MQ=255)
                                             aGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTACACGATGCCAAGCGTCGTCAGGAGGTGTTGTATCAGTGTTCGCAGgcg                                                                                                     >  5:62907/1‑149 (MQ=255)
                                              gCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCTAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTCTTGTATCAGAGTTCGCAGgcgc                                                                                                    >  2:147168/1‑149 (MQ=255)
                                              gCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGgcgc                                                                                                    >  8:82236/1‑149 (MQ=255)
                                               cTTGGTCTGCTTTCAAAAACATTTAACCAGATGTCGATCGAGCTTCATAAATTGTACCGTTCGTTTGTTTTTTAATTATAAGTATAGACCCGAGATCTACACGAGGCCAAGAGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCt                                                                                                   <  6:62907/149‑1 (MQ=255)
                                               cTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGgcg                                                                                                     <  3:261429/147‑1 (MQ=255)
                                               cTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGgcg                                                                                                     >  4:261429/1‑147 (MQ=255)
                                                ttGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGTGTCGTCTGGAGGtgttg                                                                                                                       <  7:90723/128‑1 (MQ=255)
                                                ttGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGTGTCGTCTGGAGGtgttg                                                                                                                       >  8:90723/1‑128 (MQ=255)
                                                  ggTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGGTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCACGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGaa                                                                                                >  6:231070/1‑149 (MQ=255)
                                                   gTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGaa                                                                                                >  3:285333/1‑148 (MQ=255)
                                                   gTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAACCGTCGTCTGGAGGTGTTGAATCAGTGTTCGCAGGCGCTGaa                                                                                                >  2:400852/1‑148 (MQ=255)
                                                      tGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGATCGCTGGAAGCGTCAGTATAAGAAAAGACCCCCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACt                                                                                            <  1:204190/149‑1 (MQ=255)
                                                       gCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGgtgt                                                                                                                         >  7:206428/1‑119 (MQ=255)
                                                       gCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGgtgt                                                                                                                         <  8:206428/119‑1 (MQ=255)
                                                        cTTTCAAAAACCTTTAACCAGATGTAGAGCGAGGTTCATAAATTTTACCGTTAGATTGAATAGTAATTAGAAGAATAGAACCGAGATCTACACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTTAACACTa                                                                                           <  2:270583/148‑1 (MQ=255)
                                                        cTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACt                                                                                            <  7:82236/147‑1 (MQ=255)
                                                         ttGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGc                                                                                         >  6:214369/1‑149 (MQ=255)
                                                           gCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGc                                                                                         >  4:135884/1‑147 (MQ=255)
                                                           gCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGc                                                                                         >  8:303474/1‑147 (MQ=255)
                                                             aaaaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAATAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGa                                                                                     <  4:285333/149‑1 (MQ=255)
                                                             aaaaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGa                                                                                     <  5:382909/149‑1 (MQ=255)
                                                             aaaaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACAGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCAACGAGGCCAAGCTTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCAGAACACTAGc                                                                                         >  1:103196/1‑145 (MQ=255)
                                                              aaaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGAt                                                                                    >  8:347996/1‑149 (MQ=255)
                                                              aaaaCCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTAGTTTTTTTTTTAATTATAAGTATAGCACCGAGATCTACACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGAt                                                                                    <  4:364660/149‑1 (MQ=255)
                                                                  ccTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAt                                                                                <  6:130780/149‑1 (MQ=255)
                                                                  ccTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGAGTTCGCAGGCGCTGAACACTAGCCAGATTGAt                                                                                <  1:400852/149‑1 (MQ=255)
                                                                   cTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAtg                                                                               >  7:216844/1‑149 (MQ=255)
                                                                   cTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAtg                                                                               >  2:24824/1‑149 (MQ=255)
                                                                    tttAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGAACCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGTCGCTGAACACTAGCCAGATTGAtgt                                                                              >  6:85242/1‑149 (MQ=255)
                                                                      tAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGAtgt                                                                              <  1:24824/147‑1 (MQ=255)
                                                                        aCCAGATGTCGAGCGAGCTGCAGAAATTGTACCGTTCGCTGGCAGCGTCAGTAGAAGAAAAGACCCGCGATCTAAACGAGTCCAAGCGGCTTCTGGAGTTGTTGTATCAATGTTCGCAGGCGCTGAACACTAGCCAGAATTATGTTCAt                                                                          >  8:293849/1‑149 (MQ=255)
                                                                         ccAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCAtt                                                                         >  1:114860/1‑149 (MQ=255)
                                                                                    gcgagcTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCat                                                              <  8:216844/149‑1 (MQ=255)
                                                                                      gagcTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCACGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGCACACTAGCCAGATTGATGTGCATTGTTTCCGCCatat                                                            >  4:345931/1‑149 (MQ=255)
                                                                                       agcTGCATAATTTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATAtt                                                           <  5:66954/149‑1 (MQ=255)
                                                                                       agcTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATAtt                                                           <  5:214369/149‑1 (MQ=255)
                                                                                       agcTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATAtt                                                           >  3:134484/1‑149 (MQ=255)
                                                                                       agcTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATAtt                                                           >  4:174997/1‑149 (MQ=255)
                                                                                            cATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCa                                                      >  7:131785/1‑149 (MQ=255)
                                                                                            cATAAATTGTACAGTTCGCTTGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGTAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCATATTGATGTGCATTGTTTCCTCCATATTTTGCa                                                      <  5:85242/149‑1 (MQ=255)
                                                                                                     tACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGCACACTAGCCAGCTTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTcg                                             >  1:396722/1‑149 (MQ=255)
                                                                                                     tACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTcg                                             >  5:284596/1‑149 (MQ=255)
                                                                                                     tACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTcg                                             >  8:196328/1‑149 (MQ=255)
                                                                                                       ccGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTcgc                                            >  3:127845/1‑148 (MQ=255)
                                                                                                          ttCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGAc                                         >  1:295630/1‑148 (MQ=255)
                                                                                                             gCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAAt                                      <  2:396722/148‑1 (MQ=255)
                                                                                                                 gtaGCTTCAGTAGAAGAAAAGACCCGCTATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAg                                  <  4:127845/146‑1 (MQ=255)
                                                                                                                 gAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGaa                                   <  1:147168/147‑1 (MQ=255)
                                                                                                                 gAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAg                                  <  3:135884/148‑1 (MQ=255)
                                                                                                                 gAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGc                                 >  4:177076/1‑149 (MQ=255)
                                                                                                                          gTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATAAGTGTTCGAAGGCGCTGAACACTAGCCAGAGTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGGTCGCGCCAATGAAGCGGCTGAata                        >  8:119218/1‑149 (MQ=255)
                                                                                                                              aagaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATc                      >  2:145298/1‑147 (MQ=255)
                                                                                                                              aagaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATc                      <  1:145298/147‑1 (MQ=255)
                                                                                                                              aagaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTg                    <  7:347996/149‑1 (MQ=255)
                                                                                                                                gaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATAtgtg                    >  6:2805/1‑147 (MQ=255)
                                                                                                                                gaaAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGa                  <  4:134484/149‑1 (MQ=255)
                                                                                                                                  aaaGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGtt                                               >  8:60414/1‑118 (MQ=255)
                                                                                                                                  aaaGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGtt                                               <  7:60414/118‑1 (MQ=255)
                                                                                                                                   aaGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGtt               <  5:231070/149‑1 (MQ=255)
                                                                                                                                       cccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAAt           <  3:295283/149‑1 (MQ=255)
                                                                                                                                       cccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAAt           >  8:330255/1‑149 (MQ=255)
                                                                                                                                       cccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGAAAGAGTGATGTGCATTGTTTCCTCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAAt           >  8:264847/1‑149 (MQ=255)
                                                                                                                                        ccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCt                             <  2:192820/130‑1 (MQ=255)
                                                                                                                                        ccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCt                             >  1:192820/1‑130 (MQ=255)
                                                                                                                                        ccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATg          >  1:403840/1‑149 (MQ=255)
                                                                                                                                        ccGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGAAAATGAAGCGGCTGAATATCTGCAGTTAAATg          >  5:29102/1‑149 (MQ=255)
                                                                                                                                           cgATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCg       <  3:85220/149‑1 (MQ=255)
                                                                                                                                           cgATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCg       <  3:108959/149‑1 (MQ=255)
                                                                                                                                             aTCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCGGt     >  1:361145/1‑149 (MQ=255)
                                                                                                                                                tcCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGa                           <  6:148189/124‑1 (MQ=255)
                                                                                                                                                tcCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGa                           >  5:148189/1‑124 (MQ=255)
                                                                                                                                                 ccACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCGGTGaa  >  3:147190/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACGGGCAGTTCGACTCGCCGCCGCTGGATACCAACCTGCCGAATGAGCTTGGTCTGCTTGCAAAAACCTTTAACCAGATGTCGAGCGAGCTGCATAAATTGTACCGTTCGCTGGAAGCGTCAGTAGAAGAAAAGACCCGCGATCTCCACGAGGCCAAGCGTCGTCTGGAGGTGTTGTATCAGTGTTCGCAGGCGCTGAACACTAGCCAGATTGATGTGCATTGTTTCCGCCATATTTTGCAGATTGTTCGCGACAATGAAGCGGCTGAATATCTGGAGTTAAATGTCGGTGAA  >  NC_000913/2586305‑2586597

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