Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,474,290 T→C K88R (AAA→AGA)  rpsL ← 30S ribosomal subunit protein S12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,474,2900TC97.3% 263.3 / ‑1.7 75K88R (AAA→AGA) rpsL30S ribosomal subunit protein S12
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (38/34);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

AGAATTTTACGCTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACG  >  NC_000913/3474044‑3474519
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
aGAAGTTTACGCTGACCAATGACGCGACGAGGTGGTATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGTACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGTTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAac                                                                                                                                                                                                                                   <  2:242025/251‑1 (MQ=255)
     tttACGCTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGAc                                                                                                                                                                                                                              >  1:13843/1‑251 (MQ=255)
           cTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCAc                                                                                                                                                                                                                        >  1:49038/1‑251 (MQ=255)
                cAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACggatc                                                                                                                                                                                                                   >  1:168513/1‑251 (MQ=255)
                cAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACggatc                                                                                                                                                                                                                   <  1:57471/251‑1 (MQ=255)
                  aTGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACggatc                                                                                                                                                                                                                   >  1:174334/1‑249 (MQ=255)
                     aCGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTCCGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACggatcaggat                                                                                                                                                                                                              >  1:3828/1‑251 (MQ=255)
                       gcgACGTCGGGGCATGGAAATACTCCTTTGTTAATTCAGGATTGTTCAAAACTCTACTAGTTTAGTTTGAAATTTATTTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCCCGCCATACTTGGACCGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCCCGTACGGTGTGGTAACGAACACCCGTGAGGTCTCTAACACGACCGCCACggatcaggatca                                                                                                                                                                                                            <  2:252804/251‑1 (MQ=255)
                         gacgacGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTGCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACg                                                                                                                                                                                                          <  2:49038/251‑1 (MQ=255)
                                 ggCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCATGGATCAGGATCACGGAGTGCTc                                                                                                                                                                                                  >  2:213436/1‑251 (MQ=255)
                                              tCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTCACACTACCGCCACGGATCAGGATCACGGAGTGCGCCTGCAGGTtgtg                                                                                                                                                                                      >  2:120892/1‑250 (MQ=255)
                                                 gttgttAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAACACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCCTACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCCGGTTGTGACCt                                                                                                                                                                                  >  2:25649/1‑251 (MQ=255)
                                                        aTTCAGGATTGTCCAAACCTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTcaccac                                                                                                                                                                           <  1:295870/251‑1 (MQ=255)
                                                          tCAGGATTGTCCAAAACTCTACTAGTTTAGTTTGACTTTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTCCGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCaccacc                                                                                                                                                                          <  1:132747/250‑1 (MQ=255)
                                                          tCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATGAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCaccacc                                                                                                                                                                          <  2:76592/250‑1 (MQ=255)
                                                               aTTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTATAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCATTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCAAGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGt                                                                                                                                                                     >  2:209787/1‑250 (MQ=255)
                                                                 tGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACAGCCACggaacag                                                                                                                                                                                                                 >  1:365258/1‑204 (MQ=255)
                                                                 gtcccAAAACCCTACGAGTTTAGTTTGACATTTAAGTTAAACCGTTTGCCCTTCCTTACCGGGAAACCATTAGGCCTTAGGACGCTTCCCGCCATACTGGGCACGACCCTGCTTACGGTCTTTAACGCCGGAGCATTCAAGCGCACCCCGTACGGTGTGGTAACGAACCCCCGGGAGGTCTCTCTCACGCCCGCCACggatcag                                                                                                                                                                                                                 <  2:365258/201‑1 (MQ=255)
                                                                   tCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGa                                                                                                                                                                 >  2:347285/1‑250 (MQ=255)
                                                                              aCGAGTTTAGTTTTGCATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGTTGTGGGAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGaa                                                                                                                                                     <  1:21970/251‑1 (MQ=255)
                                                                              aCGAGTTTAGTTGTACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGTTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTc                                                                                                                                                        <  2:398663/248‑1 (MQ=255)
                                                                                 agtttagtttGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAAcc                                                                                                                                                  >  1:115841/1‑251 (MQ=255)
                                                                                   tttagtttGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATCTAGGAAGTCACTTCGAAACCGt                                                                                                                                                >  1:313585/1‑251 (MQ=255)
                                                                                   tttagtttGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGt                                                                                                                                                >  1:43255/1‑251 (MQ=255)
                                                                                      agtttGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGtt                                                                                                                                               >  1:411013/1‑249 (MQ=255)
                                                                                                ttAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGaa                                                                                                                                    <  1:158330/250‑1 (MQ=255)
                                                                                                  aaGTTAAAACTGTTGGCTTTATTTAACGGAAAACCATTAAGCCGTAGGACCCTTCACGCCATACTTGGAACGAGCCTGATTCCGGTCTTTAACGCCGGAGCAGTCAAGCGCCCCACGTACGGTGTGGGAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAaca                                                                                                                                  <  2:160624/250‑1 (MQ=255)
                                                                                                   aGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACg                                                                                                                                >  1:156809/1‑251 (MQ=255)
                                                                                                    gTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACgg                                                                                                                               <  1:256751/251‑1 (MQ=255)
                                                                                                      tAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGCCGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTTCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCCGACGAACACGGc                                                                                                                              >  1:351296/1‑250 (MQ=255)
                                                                                                      tAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAAAACGGc                                                                                                                              >  2:332870/1‑250 (MQ=255)
                                                                                                      tAAAACGTTTGGCCTTACTTAACGGAGAACCAGTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGc                                                                                                                              <  1:148514/250‑1 (MQ=255)
                                                                                                       aaaaCCTTTGGCCTTAATTACCGGAGAACCATTAGGCTTTCGAAAGCTTCACGCCATACTTGGAACGAGCATGCTTACCGTCTTTAACGCAGGACCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCCCGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGACGTCACTTCGACACCGTTAGTCAGACGAACACGGCAt                                                                                                                            >  2:393579/1‑251 (MQ=255)
                                                                                                        aaaCGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGACCACGGCAt                                                                                                                            >  1:336154/1‑250 (MQ=255)
                                                                                                            gTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACttt                                                                                                                       >  1:13786/1‑251 (MQ=255)
                                                                                                               tGGCCTTTCTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTGGGAACGATCCTGTTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATCTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTAc                                                                                                                     <  2:313585/250‑1 (MQ=255)
                                                                                                                    ttacttaACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAgcg                                                                                                               >  1:281795/1‑251 (MQ=255)
                                                                                                                      acttaACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCCCCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGTCCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGcgcg                                                                                                             <  1:390670/251‑1 (MQ=255)
                                                                                                                          cacggagaaaccttacgacatAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGAACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAAACGAACACGGCATACTTTACGCAGCGCGGAGt                                                                                                         >  2:55209/21‑251 (MQ=255)
                                                                                                                          aaCGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGt                                                                                                         >  1:348742/1‑251 (MQ=255)
                                                                                                                          aaCGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGt                                                                                                         >  1:4840/1‑251 (MQ=255)
                                                                                                                          aaCGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTCGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGc                                                                                                                                                                                              <  1:260245/166‑1 (MQ=255)
                                                                                                                          aaCGGAGAACCATTAAGCCTAAGGACGCATCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCAGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGc                                                                                                                                                                                              >  2:260245/1‑166 (MQ=255)
                                                                                                                            cgagaaaCCATTAACCCTTAGGACGCTTCCCGCCATACTTGGACCGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTc                                                                                                       <  2:13786/246‑1 (MQ=255)
                                                                                                                             ggAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCCCCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCCTACTTTACGCAGCGCGGAGTTCg                                                                                                      >  2:69292/1‑251 (MQ=255)
                                                                                                                                        aaGCCTTAGGAAGCTACCCGCCATACTTGGAACGCGCTTACTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAggg                                                                                           >  2:66196/1‑250 (MQ=255)
                                                                                                                                         aGCCTAAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAg                                                                                          <  2:174334/251‑1 (MQ=255)
                                                                                                                                          gCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGt                                                                                         >  2:141582/1‑251 (MQ=255)
                                                                                                                                           ccTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTg                                                                                        <  1:21676/251‑1 (MQ=255)
                                                                                                                                              tAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGgt                                                                                      >  1:332186/1‑250 (MQ=255)
                                                                                                                                              tAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGgt                                                                                      >  1:207772/1‑250 (MQ=255)
                                                                                                                                                   cGCTTCACGCCATACTTGGAACGGGCCGGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACCCCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGtata                                                                                <  2:211365/251‑1 (MQ=255)
                                                                                                                                                    gCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGtatat                                                                               >  2:206274/1‑251 (MQ=255)
                                                                                                                                                       tCACGCCATACTGGGAACGAGCCTGCTTACGGGCTTTACCGCCGGAGCAGTCAGGCGCACCACGTCCGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGtatata                                                                              <  2:4840/249‑1 (MQ=255)
                                                                                                                                                            ccgtcctgggAACGACCCGGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTAGGGTTGGGTACCGAACGCCCGGGAGGTCCCGAACTCGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGTTTGTTACCTTCACCACCGATGTAGGATGTCACTTCGAAACCGTTAGTCAGACGAACACGGCAGTCTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGt                                                                       <  1:246172/243‑1 (MQ=255)
                                                                                                                                                                   tGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCCCCACGTACGTTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATAc                                                                 <  1:47536/250‑1 (MQ=255)
                                                                                                                                                                              tGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTCCGCCGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATAc                                                                 >  2:344616/1‑239 (MQ=255)
                                                                                                                                                                              tGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATAc                                                                 <  1:344616/239‑1 (MQ=255)
                                                                                                                                                                                     ggTCGTTAACGCCGGAGCAGTAAAGCGCACCACGTACGGGGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCAt                                              <  2:207772/251‑1 (MQ=255)
                                                                                                                                                                                              cGCCGGAGCAGTCAAGCGCACCACGGACGGTGGGGTAACGAACCCCGGGGAGGTCTCTAGCACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAgc                                     <  2:281795/251‑1 (MQ=255)
                                                                                                                                                                                               gCCGGCGCAGTCACGCGCACCTCGTACGGTGTGGCAACGAACACCCGGGAGGTCTCTAACACGAAAGCCAAGGATAAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTAACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGAATGCTTCCAgcg                                    >  2:124759/1‑251 (MQ=255)
                                                                                                                                                                                               gCCGGAGCAGTCAAGCGCACCACGGACGGTGTGGTAACGAACCCCCGGGAGGCCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACTCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAgcg                                    <  2:292586/251‑1 (MQ=255)
                                                                                                                                                                                                ccGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAgcgc                                   <  1:69292/251‑1 (MQ=255)
                                                                                                                                                                                                       aGTCAAGCGCCCCACGTACGGGGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACg                            <  2:380589/251‑1 (MQ=255)
                                                                                                                                                                                                         tCAAGCGCAACACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGtt                          >  2:286739/1‑251 (MQ=255)
                                                                                                                                                                                                          cAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGTAGGTCTCTAAAACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGGTGTAGGAAGTCACTTCGACACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGAGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGACATGCTTCCAGCGCAGGCACGTTg                         >  1:371239/1‑251 (MQ=255)
                                                                                                                                                                                                          aaagagcACCACGTACGGTGTGGAAAAGAACACCGGGGAGGGATATAACACGACCGCCACGGATCAGGATAACGGAGGGCTCATGAAGGTTGTGACCTGCACAAACGAGGTAGGAAGTCACTGCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTGCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTGTTGCGGGCATGCTTCCAGCGCAGGCACGTTg                         <  2:96482/246‑1 (MQ=255)
                                                                                                                                                                                                               gcACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCtttt                    <  1:214121/251‑1 (MQ=255)
                                                                                                                                                                                                                caccacGTACGGTGTGGTAACGAACGCCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCTCGTTGCTTTTc                   <  1:259024/251‑1 (MQ=255)
                                                                                                                                                                                                                 cccccgtccgtgggggTAACGAACACCCGGGCGGTCTCTAACACGCCCGCCACGGATCAGGAGCACGGAGTGCTCCTGCAGGTTGTGACCTGCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTCACGCAGCGCTGAGTTCGGTTGTTTAGGAGTGGTAGTATATACACGAGTACAGACGCCACGTTTTTGCGGGCATCCTTCCAGCGCAGGCACGTTGCTTTTCg                  <  2:404931/237‑1 (MQ=255)
                                                                                                                                                                                                                  ccacGTACGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGCAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGc                 >  1:317623/1‑251 (MQ=255)
                                                                                                                                                                                                                        aCGGTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACttt           >  2:157427/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggTGTGGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGCTGGGACCTAAACCACCGATGTCGGACGTCACTTCGAAACGGATAGTCCGACGAACACGGCATACTTTACGCAGCGCGGAGCGCGGTTTTTTAGGAGTGGAAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGc         >  1:368691/1‑251 (MQ=255)
                                                                                                                                                                                                                             gtgGTAACGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAg      <  2:156809/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCGAACACCCGGGAGGTCTCTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTCACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACg  <  1:347285/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     
AGAATTTTACGCTGACCAATGACGCGACGACGTGGCATGGAAATACTCCGTTGTTAATTCAGGATTGTCCAAAACTCTACGAGTTTAGTTTGACATTTAAGTTAAAACGTTTGGCCTTACTTAACGGAGAACCATTAAGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACG  >  NC_000913/3474044‑3474519

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