Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,543 C→A 100% P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA100.0% 270.3 / NA 81P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (47/34);  total (47/34)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAA  >  NC_000913/4184317‑4184767
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
tgcGAGAATAAACTCGAAGCGAAACGACGAAAAATCACCTAGGGCGACTATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                                                                                                                 <  1:542907/227‑1 (MQ=255)
  cGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTCATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGAGGAAAACGGTACGCCGGAAGAAAACGTAATAAACCAGCt                                                                                                                                                                                                                                >  2:264500/1‑228 (MQ=255)
       agaaACGCGAAGCGGAACGGCGGAAAAGCACCCAGGGCGACGATCTGGCACCGGGCGTAGTGAATATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                          <  1:425139/229‑1 (MQ=255)
            cgcGAAGCGATGCGCCTCAAAATCACCCAGGGCGCCGATCGGGCACCGGGCGTGCTGAAGATTGTCAAGGTATATCTGGCGGTTAAACGACGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTGCTAAGATCCACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                                                                                                                     <  1:258795/227‑1 (MQ=255)
            cTCGAAGGTAAACGCCGGAAAATCCCCTAGTGCTGCGATCTGGCCCCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATACAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTGCTAAGATGAACCCGATCGAAGATCAGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                                                                                                                     <  1:498560/229‑1 (MQ=255)
                      aaCGGCGAAAAAACACCCAGGGCGACGATCCGGCCCCGGGCGTGGTGAAGATTGGTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                                                                                                           <  1:285597/229‑1 (MQ=255)
                       aCGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGGATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCACCTGGGGGTACCGCCCCGTATGac                                                                                                                                                                                                           >  1:206335/1‑227 (MQ=255)
                        cgccgcAAAATCCCCCAGGGCGACGAGCTGGCACCGGGCGTGCTTAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                                                                                                         <  1:388620/229‑1 (MQ=255)
                        agcggcAAAAGCCCCCAGGGGGCCGATCTGGCACTGGGCTTGCGGAAGTTTGTTAAGGGATATCTGGTGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                                                                                                         <  2:135201/225‑1 (MQ=255)
                                 aTCACCCGGGGCGACGATTTGGGACCGGGGGTGCTGAGGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTTTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                                                                                                                <  2:85163/229‑1 (MQ=255)
                                 aTCACCCAGGGCGACGATCGGGCACCGGGCGTGCTGAAGATTGGTAAGGTATAGCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGAGGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                                                                                                                <  1:189709/229‑1 (MQ=255)
                                     cccAGGGCGACGATTTGGCCCCGGGCGTGCCGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                                                                                                                            <  2:219953/229‑1 (MQ=255)
                                          ggCGACGAGCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTCCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                                                                                                       <  1:524981/229‑1 (MQ=255)
                                            cgagggtCTGGCCCCGGGCGTGATGAAGATTGGTAAGGTATATTGGGCGGTTAAACGCCGGATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                                                                                                                     <  2:213761/223‑1 (MQ=255)
                                               cgaTCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCTGATGAACATCGGTCAGATCCTCGaacacca                                                                                                                                                                                   >  1:247626/1‑223 (MQ=255)
                                               cgaTCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCGGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                  <  1:162181/229‑1 (MQ=255)
                                                            ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGCTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATACGCCTTACGATGACAACGGTACGCCGGTAGACATCGTACTGAACCAACTGGGCGTACCGCCTCGTATGACAACCGGCCAGATCCCCGAAACCCCCCTGGGGATGGccg                                                                                                                                                                     >  2:399387/1‑227 (MQ=255)
                                                            ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTTGGCGTACCGTCTCGTATGAACGTCGGTCAGATCCTCGAAACCCACCTGGGTTTGGCt                                                                                                                                                                      >  1:29845/1‑228 (MQ=255)
                                                            ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCCCCTGGGGATGGCt                                                                                                                                                                      >  1:386214/1‑228 (MQ=255)
                                                            ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGAGGAAAACGGTACGCCGGTAGACATCCTACTGAACCAGCAGGGCGTACCGTCCCTGATGAAAAACGGCCAAAACCTCGAAACCCACCTGGGGATGGcgg                                                                                                                                                                     >  2:659732/1‑227 (MQ=255)
                                                            ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGGACTGAACCAGCTGGGCGTACCGTCTCGCATGAACATCGGTCAGATCCCCGAAACCCACCTGGGTACGGcg                                                                                                                                                                      >  1:297354/1‑227 (MQ=255)
                                                                          tGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGAACCGTCCCGTATGAACATCGGTCAAATCCTCGAAACCCCCCTGGGCATGGCTGCGAAAGGGATccga                                                                                                                                                       >  2:683078/1‑226 (MQ=255)
                                                                           gTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACACCGGTAAGATCCTCGACACCCACCTGGGGATGGCTGCGAAAGGTACCGGc                                                                                                                                                       >  1:206962/1‑228 (MQ=255)
                                                                            ttAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                                                                                                                        <  1:472273/162‑1 (MQ=255)
                                                                            ttAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                                                                                                                        >  2:472273/1‑162 (MQ=255)
                                                                                   atatTTGGCGGTTAAACGCCGTATGCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                                                                                                               <  2:312908/228‑1 (MQ=255)
                                                                                    tatCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGGACATCGGTCAGATACTCCGAACCCACCTGGGAATGGCTGCGAAAGGTATTGGCGACCAGaaaa                                                                                                                                             >  2:749348/1‑226 (MQ=255)
                                                                                        gggCGGTTAAACGCGGTTTCAGGGCTGGTGACAAGGTGGCAGGGGGTACCGGGAACAAGGGTGTAATTTCCAAGATCAACCCGATGGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                                                                                                                         <  2:245629/228‑1 (MQ=255)
                                                                                          gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTTTGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGGGAAAGGAATTGGGAACaaaacacaaccca                                                                                                                                       >  2:478525/1‑218 (MQ=255)
                                                                                          gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACAGCCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAAGGATCGGCGCACACACCAACGCCa                                                                                                                                      >  2:478553/1‑229 (MQ=255)
                                                                                           cGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGCCGACAAGATCAACGCCAt                                                                                                                                     >  1:496468/1‑229 (MQ=255)
                                                                                                aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGTAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCCCCTGGGCATGGCTGCGCAAGGGATCGGCGACAAGGTCAACGCCATTCTGa                                                                                                                                >  1:394219/1‑229 (MQ=255)
                                                                                                  aCGGCGTAGCCAGGCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTGTAAGATCAACCCGATCGAAGATATGCCCTACGAGGAAAACGGTACGGCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGTTATCGGCGACAAGATCAACGCCATGCTGaaa                                                                                                                              <  1:2787/229‑1 (MQ=255)
                                                                                                  aCGCCGTATCCAGGCTGGTGACAAGATGGTAGGGAGTCAAGGTAACAAGGGTGTAATTTCTAAGGTTAACCCGGTCGAAGATAGGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                                                                                                              <  2:386214/229‑1 (MQ=255)
                                                                                                  aCGCCGTATCCAGCCTGGTGACAAGAAGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                                                                                                              <  1:170859/229‑1 (MQ=255)
                                                                                                             cgCCTGGTACAAAGATGGCGGGGCGGCACGGTAACTAGGGGTTAATTTCTAAGAACAACCCGATCGAAGATATGCCTTACGATGAAAACGGTGCGCCGGTAGGCATCGTACTGAACCATCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAACCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                                                                                                                   <  2:397430/228‑1 (MQ=255)
                                                                                                                cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGAGTGGCGGCGAAAGGAATCGGCGACAAGACAAAAGCCACGCTGAAACAGCAGCaagaagg                                                                                                                >  2:736071/1‑228 (MQ=255)
                                                                                                                cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGAACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGCACGGGCGACAAGATCAACGCCATCCTAAAACAACAGCaaaaaaa                                                                                                                >  2:688855/1‑227 (MQ=255)
                                                                                                                          gATGGCAGGTCGTCACGGTAACAAGGGTGTAAGTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa                                                                                                         <  1:632760/226‑1 (MQ=255)
                                                                                                                           aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTg                                                                                                      >  1:64695/1‑228 (MQ=255)
                                                                                                                            tGGCAGGTCGTCCAGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg                                                                                                    <  1:378234/229‑1 (MQ=255)
                                                                                                                            tGGCAGGTCGACGAGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg                                                                                                    <  1:556595/229‑1 (MQ=255)
                                                                                                                             ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCCCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGCTATGGCTGCGAAAGGTATCCGCGACAAGATCAACCCCAATGTTGAACCACAACAAAAAGGAGCGAAAccacgg                                                                                                   >  2:300565/1‑224 (MQ=255)
                                                                                                                                 ggTCGGTACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa                                                                                                <  1:153456/228‑1 (MQ=255)
                                                                                                                                  gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGCCAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTGGCGAAACTGCGCGAAtt                                                                                              >  1:271355/1‑229 (MQ=255)
                                                                                                                                  gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGGGAAtt                                                                                              >  1:34401/1‑229 (MQ=255)
                                                                                                                                       cACGGTAACAAGGGTGTAATTTCCAAGAGCAACCCGATCGAAGATATGCATTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATcc                                                                                         <  2:642462/229‑1 (MQ=255)
                                                                                                                                        aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGCCAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGAGCGAATTCATCCa                                                                                        >  1:161638/1‑229 (MQ=255)
                                                                                                                                             aaCAAGGGTTTAATTTCTAAGATCAAGCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcg                                                                                     <  1:419892/227‑1 (MQ=255)
                                                                                                                                               cAAGGGTTAAATTTCAAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTGCGCCGGTAGACATCGTACTGAACCAACTGGGCGTACCGTCTCGTATGAACATTGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcgtgcg                                                                                 <  2:554745/229‑1 (MQ=255)
                                                                                                                                               cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAAATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGACGCGGAACTGCGCGCATTCATCCAgcgtgcg                                                                                 >  1:399053/1‑229 (MQ=255)
                                                                                                                                                    gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACACCATGCTGAAACAGCAGCAAGGAGTCGCGAAAACGCACGAAATCATCCAGGGAGCGAAcaa                                                                            >  2:247470/1‑227 (MQ=255)
                                                                                                                                                      gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTGGCGAAAATGCGCGCATTCATCCAGCGTGGTTACGaa                                                                           >  1:628338/1‑227 (MQ=255)
                                                                                                                                                      gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCATCGTGGGTACGATc                                                                          >  1:521303/1‑229 (MQ=255)
                                                                                                                                                      gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAAGCTGAAACAGCAACAAGAAGTCGCGAAACTGCCGCAATTCATCCAacgtcgggagaaac                                                                          >  2:91959/1‑219 (MQ=255)
                                                                                                                                                       tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATc                                                                          <  1:399387/228‑1 (MQ=255)
                                                                                                                                                          tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGAAAGAAGTCCCGAAACTGCGCGAATTAATCCAGCGGGCGCACGATCTgg                                                                       >  1:412911/1‑228 (MQ=255)
                                                                                                                                                            tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTgggg                                                                     >  1:207919/1‑227 (MQ=255)
                                                                                                                                                               aaGATCAACCCGAACGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTCCCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa                                                                                                <  1:417858/198‑1 (MQ=255)
                                                                                                                                                                                      gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGTAAGTTGACCTg                                          >  1:618354/1‑229 (MQ=255)
                                                                                                                                                                                         ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCAAGCGCGCGAACGATCTGGGCGCTGACGTTCGTCAGAAAGGTGACCTgtgt                                       >  1:704371/1‑229 (MQ=255)
                                                                                                                                                                                         ttACGATGAAAACGGTACGACCGTAAAGTGCTGCATGACACAGCTGGGCGACCCGGTTCGTATGAACATCGATCAGTTCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTCACGTTCGTCCGAAAGTGGTCATGAGt                                       <  2:394219/229‑1 (MQ=255)
                                                                                                                                                                                          tACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTCGACCTGAgca                                      >  1:684644/1‑227 (MQ=255)
                                                                                                                                                                                            cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGCCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGCTCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGCAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTGCGTCCGAAAGTTCACCTggctaac                                    >  1:146496/1‑223 (MQ=255)
                                                                                                                                                                                            cGATGAAAACGGCACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACCAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTCCGCGAATTCATCCCGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTTACCTGAGTAcc                                    >  1:121972/1‑229 (MQ=255)
                                                                                                                                                                                                 aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                                                                                                       >  1:670828/1‑78 (MQ=255)
                                                                                                                                                                                                 aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                                                                                                       <  2:670827/78‑1 (MQ=255)
                                                                                                                                                                                                 aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATAAACGCCATGCTGCAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGGCCGTACGACCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTCAGTACCGTCAg                               >  1:266444/1‑229 (MQ=255)
                                                                                                                                                                                                    aCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGCAGACAAGATCAACGCCATGCTGAAACAGCAGCCAGAAGTCGCGACACTGCGCGAATCAATCCAGCGTGCGTACGATCGGGGCGCTGACGTTCGCAAGCAAGTTGACCTGGGGACCTTAAgggc                            >  2:278171/1‑226 (MQ=255)
                                                                                                                                                                                                      ggTACGCCGGTAGCCCTCGTACTTAACCAGCTTGGCGTACCGGCTCGTATGAACATCGGTCAGATCCCCGAAACCCACCTGGGTAGGGCTGCGAAAGGTATCGGCTACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATg                          <  2:161638/229‑1 (MQ=255)
                                                                                                                                                                                                       gTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTCGGTACCTCCACCGATgt                         >  1:527850/1‑228 (MQ=255)
                                                                                                                                                                                                                     aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAACGGCGACAAGATCAACGCCATGCTAAAACAGCAGCAAGAAGGACCGAAACTGCGCGAATTCCTCACGCGTGCGGACGCACTGGGCGCTGACGGTCGACGTGAAGGTGACCTGGGTACCTTCAGCGGTGAGAACGTTAtcctt            >  2:445244/1‑224 (MQ=255)
                                                                                                                                                                                                                       cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                                                                                                                    <  1:504804/91‑1 (MQ=255)
                                                                                                                                                                                                                       cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGAGCTGGGCGCTGCCGTCAGTCAGAAAGTTGACCTGAGGACCTTCAGCGTGGAAGAAGTTATCGGTCt          >  1:678810/1‑228 (MQ=255)
                                                                                                                                                                                                                       cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGAGAAAGGTATCGGCGAc                                                                                                                                                    >  2:504804/1‑91 (MQ=255)
                                                                                                                                                                                                                         tACTGCACCAGCTGTGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAAAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGCCGCTGACGTTCGCCAGAAAGTTGACCTGAGTACCTTCATCCATGAAGAAGTTATTCGACTgt        >  1:197730/1‑227 (MQ=255)
                                                                                                                                                                                                                           cTGAACCAGGCGGGCGGTCCGTTGTGTATAAACATCGGTCAGATCCCCGAAACCCGCCTGGGTATTGCTGACAAAGGTATCGGCGACAAGAACAACGCCATGCCGAAACAGCAGCAAGAAGGCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTg     <  2:115211/229‑1 (MQ=255)
                                                                                                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACACGATCAACGCCATCCTGAACCAGCAGCAAGAAGTCGCGAAACTGCGCGACCTCATCCAGCGTGCGTACGATCTGGGCGCTGCCGTTCGTCAGTAAGTTGCCCTGAGTACCTTCTGCGATGCAGAAGTTATGCGTCTGGCTGaaa  >  1:437034/1‑229 (MQ=255)
                                                                                                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCACAGCCATGCTGAAACAGCGGCAAGAAGTCGCGAAACTGCGCGAATTCCTCCAGCGTGAGTACGACCGGGGCGCGGACGTTCGTCATAAAGTTGGCCTGCGTACCTTCAGCGATGAAGACGTTATGGGTCTGGCTGaaa  >  2:156979/1‑229 (MQ=255)
                                                                                                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGGGCGTACGAGCTGGACGCTGACGTTCGTCAGAAAGTTGCCCTGAGTACCTTCAGCGATGAGGAAGTTATGAGTCTGGCTGaaa  >  1:690907/1‑229 (MQ=255)
                                                                                                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGAATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCCAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATATGGGGGCTGACGGTCGGCAGAAAGCTGACCTGAGTACCTTCAGCGTTGAAGCAGGTATGCGTCTGGCTGaaa  >  1:384898/1‑229 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
TTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAA  >  NC_000913/4184317‑4184767

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