Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,351,519 C→T T48T (ACG→ACA dcuS ← sensor histidine kinase DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,351,5190CT96.9% 210.6 / NA 67T48T (ACG→ACAdcuSsensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (25/38);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTG  >  NC_000913/4351278‑4351761
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
gcgcTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGTGCTTCAGGATTCGAGTAGCGAATTCTTTTCATAGCTGTAACGACAATAAACATCAGAACGGTGCGTTTGCGTACGGTTTCCGCGATGTCCGGGATGCCACTCTCCTGCGGTTTTTTCTTCAATCCCTGCCGTATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCa                                                                                                                                                                                                                                           <  1:283182/251‑1 (MQ=255)
       aaGGAGGTCATCACCTTTAAATGGCTGACCAATGCGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCGGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGAtt                                                                                                                                                                                                                                     <  2:197760/250‑1 (MQ=255)
        aGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCGGGATGCCACGCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGGGTCATATCACTGATTTg                                                                                                                                                                                                                                   <  2:82197/251‑1 (MQ=255)
           aTGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGCTGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGAAGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCgg                                                                                                                                                                                                                                >  2:86388/1‑250 (MQ=255)
            tGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGTATTGCCGGCGAGTCGGCTAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGctgg                                                                                                                                                                                                                               >  1:308737/1‑248 (MQ=255)
              tcatcaCCTTTAAATGACTGACCAATACGCTGGGCTTCAGGATGCGTGTAGCGAAGACTTTGCATATCGGTAACGACAATAAAAAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCATGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCGAGCCAATCTCGCGTAATATCACTGATTTGCgaga                                                                                                                                                                                                                              >  2:372086/1‑250 (MQ=255)
               catcaCCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGaa                                                                                                                                                                                                                             <  1:202823/250‑1 (MQ=255)
                 tcaCCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTCGCCAAGACTTTGCATATCGCTCACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCCAGCCCTGACGGATTTCCGGCGCGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGt                                                                                                                                                                                                                           >  1:545338/1‑250 (MQ=255)
                  caCCTTTAAATGGCTGACCAATACGCGGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGATTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTaa                                                                                                                                                                                                                         <  1:546640/251‑1 (MQ=255)
                   cccTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTaaa                                                                                                                                                                                                                        <  1:279833/250‑1 (MQ=255)
                    ccTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGCAGTAAAt                                                                                                                                                                                                                       >  1:330159/1‑251 (MQ=255)
                     cTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGACACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATc                                                                                                                                                                                                                      >  2:474487/1‑251 (MQ=255)
                     cTTTAAATCGCTCACCCATACGCTGGGCTCCAGGATGCGAGTACCGACGACTCTGCATATCGGCACCGACAATGAACAGCAGACCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTCGCCCATCTCGTGTCATATCACTGCTTTGCGAGAAGTAAATc                                                                                                                                                                                                                      >  1:481987/1‑251 (MQ=255)
                                 gACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCGGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGGCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCaccacc                                                                                                                                                                                                          <  2:305537/251‑1 (MQ=255)
                                                  tCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGTCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCGGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAAc                                                                                                                                                                                           <  2:573213/249‑1 (MQ=255)
                                                  acAGGATGCGAGGAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGAGTACGGCTTCCGCGATGGCCTGGATGCCACTCGCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCTCCACCAATAGCACCGAGAACAg                                                                                                                                                                                         <  1:380234/250‑1 (MQ=255)
                                                    aGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTa                                                                                                                                                                                       <  1:326074/251‑1 (MQ=255)
                                                       atgcaggTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCg                                                                                                                                                                                    <  2:272767/245‑1 (MQ=255)
                                                          cGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTTCGTAAGGCTTCCGCGCTGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACAGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGTAGTAAATCAGATGCAACACCAATAGCACCGAGAACAGTATCGCAc                                                                                                                                                                                 >  2:246704/1‑251 (MQ=255)
                                                          cGAGTAGCGAAGACTTTCCATATCGGAAACGACGATAAGCAGAAGATCGTTGAGTTTGCGTACGGCTTACGCGATGGCCTGGATGCCACTATCATGAGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCAATGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCAc                                                                                                                                                                                 >  2:260705/1‑251 (MQ=255)
                                                            aGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTg                                                                                                                                                                               <  1:179555/251‑1 (MQ=255)
                                                               aGCGAAGACTTTTCATCTCGGTAACGACAATAAACTGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGGGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGat                                                                                                                                                                             <  2:496639/250‑2 (MQ=255)
                                                                   aaGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTa                                                                                                                                                                        <  2:390973/251‑1 (MQ=255)
                                                                     gACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTa                                                                                                                                                                        >  2:435052/1‑249 (MQ=255)
                                                                     aaCTTTGCATATCGGTAACGACAATAAACAGCAGAGCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCTCCGAGAACAGTACCGCACTGATCATTAAg                                                                                                                                                                      <  1:375835/250‑1 (MQ=255)
                                                                      aCTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTTCCGCACTGATCATTAAGa                                                                                                                                                                     <  1:397968/251‑1 (MQ=255)
                                                                         ttGCATATCGGTAACGACAATAAACAGCAGAGCGTTGCGTTTGCGGACGGCTTCCGCGATGGCCGGGATGCCACGCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGAt                                                                                                                                                                    <  1:353632/249‑1 (MQ=255)
                                                                            cATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCCCTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtg                                                                                                                                                               >  1:479939/1‑251 (MQ=255)
                                                                              tatCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgtg                                                                                                                                                             >  1:492039/1‑251 (MQ=255)
                                                                               atCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTAATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGTTCACtgggg                                                                                                                                                            >  2:437834/1‑248 (MQ=255)
                                                                                 cGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGt                                                                                                                                                           >  1:154406/1‑250 (MQ=255)
                                                                                 cGGTAACGACAATAAACAGCAGAGCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCGGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGt                                                                                                                                                           <  2:459152/250‑1 (MQ=255)
                                                                                   gTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGCTGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACt                                                                                                                                                        >  1:266898/1‑251 (MQ=255)
                                                                                          cAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCAACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACCGATCATTAAGATCACTGTGGTACTCAAttcc                                                                                                                                                 >  2:219684/1‑249 (MQ=255)
                                                                                           aaTAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGCGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCACTAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTc                                                                                                                                                 >  1:375866/1‑250 (MQ=255)
                                                                                               aaCAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCgg                                                                                                                                            >  2:389438/1‑251 (MQ=255)
                                                                                               aaCAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCGGCGGGTTTTGCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTTGTACTCAATTTCATCgg                                                                                                                                            <  2:308719/251‑1 (MQ=255)
                                                                                                           tttcgtttgcgtACGGCTTCCGCGATGGCCGTGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACAGAGAACAGTACCTCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTa                                                                                                                                 <  2:568784/247‑1 (MQ=255)
                                                                                                                 ttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGc                                                                                                                           <  1:425011/250‑1 (MQ=255)
                                                                                                                   gcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTGCCGCACTGATCAGTAAGATCACTGTGGTACTCAATTTCATCGGCCGTTTGCGTAACATGCgg                                                                                                                         <  2:333311/250‑1 (MQ=255)
                                                                                                                    cgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGAGTTGCGAGAAGTAAATCAGATGCGCCACCAATAGCACCGAGAACAGGACCGCACGGATCATTAAGAGCACTGTGGTACTCAATTTCATCGGAGGTGTGCGGAACATGCGGTa                                                                                                                       <  2:349193/251‑1 (MQ=255)
                                                                                                                     gtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAg                                                                                                                      <  1:441700/251‑1 (MQ=255)
                                                                                                                       aCGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGACCTGTCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCAATGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCTCCGAGAACAGTTCCGCACTGATCATTAAGATCACTGTGGTACTCATTTTCATCGGACGGTTGCGTAACATGCGGTAgg                                                                                                                     <  1:431329/250‑1 (MQ=255)
                                                                                                                         ggCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCa                                                                                                                  <  1:484862/251‑1 (MQ=255)
                                                                                                                             tCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCCTTAAGATCACTGTGGTACTCAATTTCATCGGACGCTTGCGTAACATGCGGTAGGGCaatg                                                                                                               >  2:482736/1‑250 (MQ=255)
                                                                                                                                   aTGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtct                                                                                                        <  2:512945/251‑1 (MQ=255)
                                                                                                                                     ggCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCtca                                                                                                      <  1:317531/251‑1 (MQ=255)
                                                                                                                                          ggATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCtcatca                                                                                                   <  1:76541/249‑1 (MQ=255)
                                                                                                                                           gATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCGCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCCCCAAGAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCt                                                                                                <  1:334574/251‑1 (MQ=255)
                                                                                                                                                   cccccggCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgtg                                                                                        <  2:361284/245‑1 (MQ=255)
                                                                                                                                                       cTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGAc                                                                                    >  1:440168/1‑251 (MQ=255)
                                                                                                                                                         gCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACa                                                                                   >  1:280159/1‑250 (MQ=255)
                                                                                                                                                         gCGGTTTTTCCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCAGCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACa                                                                                   <  2:23963/250‑1 (MQ=255)
                                                                                                                                                          cGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACaaa                                                                                 <  1:416398/251‑1 (MQ=255)
                                                                                                                                                              tttttCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAAttt                                                                              <  1:259394/250‑1 (MQ=255)
                                                                                                                                                                  tCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCATCTTCCTTGTGTGACAAATTTCTTa                                                                          >  2:325358/1‑250 (MQ=255)
                                                                                                                                                                     gCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGc                                                                       >  1:251257/1‑250 (MQ=255)
                                                                                                                                                                           ccTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAtct                                                                >  1:521022/1‑251 (MQ=255)
                                                                                                                                                                             tGACGGATTTCCGGCGAGTCGGCGAGGGTCCGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAtct                                                                <  2:421133/249‑1 (MQ=255)
                                                                                                                                                                                 ggATTGCCGGCGAGTCGGCGAGGGTTCGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtgatg                                                          <  2:375866/251‑1 (MQ=255)
                                                                                                                                                                                           cGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTa                                                 >  1:342055/1‑250 (MQ=255)
                                                                                                                                                                                              gTCGGCGAGGGTACGCGCCACTGCCAGCGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTc                                             <  1:7693/251‑1 (MQ=255)
                                                                                                                                                                                                           cgcgCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGa                                <  2:330159/251‑1 (MQ=255)
                                                                                                                                                                                                                       aGTGCATTGTCGGCTAGACCATCTCGTGTCATATCAATGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCAACGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGAAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTTAAAGGGAGTAAGaaaggttggttat                    >  2:552582/1‑240 (MQ=255)
                                                                                                                                                                                                                              tgttgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAgg               <  1:9254/249‑1 (MQ=255)
                                                                                                                                                                                                                               gttgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGa            <  1:86467/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ccATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGAGGAAGAATGATTGGCTATATAGGGGGAGAGActatg  >  1:490997/1‑248 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
GCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTG  >  NC_000913/4351278‑4351761

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