Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 927,612 Δ1 bp 6.5% coding (1584/1722 nt) cydC ← glutathione/cysteine ABC transporter export permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913927,6110C. 6.5% 220.0 / 7.7 62coding (1585/1722 nt)cydCglutathione/cysteine ABC transporter export permease/ATPase
Reads supporting (aligned to +/- strand):  ref base C (33/25);  new base . (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGAT  >  NC_000913/927333‑927898
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               
tAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAacac                                                                                                                                                                                                                                                                                        <  2:195625/288‑1 (MQ=255)
       gCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTACCACCGTTTtc                                                                                                                                                                                                                                                                                 >  1:427139/1‑288 (MQ=255)
       gCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTtct                                                                                                                                                                                                                                                                                <  2:315410/289‑1 (MQ=255)
       gCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTtc                                                                                                                                                                                                                                                                                 >  1:332324/1‑288 (MQ=255)
       gCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTtc                                                                                                                                                                                                                                                                                 >  1:400612/1‑288 (MQ=255)
       gCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAAGAATTTGCCCGTTGTCCATCAATATTATTTGTTGGAAACGAGAGAGTCCGCAAAGTCGTTGGGTGACCATTAAATCTGttttg                                                                                                                                                                                                                                                                                 >  1:299015/1‑287 (MQ=255)
         ccccAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTtctc                                                                                                                                                                                                                                                                               >  2:187937/1‑288 (MQ=255)
           ccAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACg                                                                                                                                                                                                                                                                            <  1:116827/289‑1 (MQ=255)
                         gCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTgcaa                                                                                                                                                                                                                                                              >  2:180756/1‑289 (MQ=255)
                         gCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTgca                                                                                                                                                                                                                                                               <  1:187937/288‑1 (MQ=255)
                           cGTTAGGTTAACGGAATGCGCCTGCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTgcaag                                                                                                                                                                                                                                                             <  1:258655/288‑1 (MQ=255)
                                  cTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTc                                                                                                                                                                                                                                                       <  2:482562/287‑1 (MQ=255)
                                     aCGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGa                                                                                                                                                                                                                                                  >  2:244896/1‑289 (MQ=255)
                                       gTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATc                                                                                                                                                                                                                                                <  2:387367/289‑1 (MQ=255)
                                              gcCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCttt                                                                                                                                                                                                                                         <  2:32541/289‑1 (MQ=255)
                                                       ggaagggaagGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATGTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTtgtg                                                                                                                                                                                                                                 <  2:299015/288‑1 (MQ=255)
                                                                            gCGAGAAAGCTGACACGATCGCAAGACCGCAACACTCCACTCAAATCGGATCGTTCAATAATAGCTAACAAACCCTGCTTGAACTGGTAATCACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGt                                                                                                                                                                                                            <  1:83825/288‑1 (MQ=255)
                                                                                         aCCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGGTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAacac                                                                                                                                                                                              <  2:331856/289‑1 (MQ=255)
                                                                                             ggCGCATGACCGCAGAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGt                                                                                                                                                                                          >  2:341695/1‑289 (MQ=255)
                                                                                               cgcATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTtgtg                                                                                                                                                                                                                                 >  2:215443/1‑248 (MQ=255)
                                                                                               cgcATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTtgtg                                                                                                                                                                                                                                 <  1:215443/248‑1 (MQ=255)
                                                                                               cgcATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTg                                                                                                                                                                                         <  1:244896/288‑1 (MQ=255)
                                                                                                  aTGACCGCAGAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGgcg                                                                                                                                                                                      <  1:341695/288‑1 (MQ=255)
                                                                                                  aTGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGgcg                                                                                                                                                                                      <  2:400612/288‑1 (MQ=255)
                                                                                                                                      taataGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGgcg                                                                                                                                                 <  2:332324/289‑1 (MQ=255)
                                                                                                                                            cTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTc                                                                                                                                           >  2:386251/1‑289 (MQ=255)
                                                                                                                                              tACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGgcg                                                                                                                                                 >  2:240164/1‑281 (MQ=255)
                                                                                                                                              tACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGgcg                                                                                                                                                 <  1:240164/281‑1 (MQ=255)
                                                                                                                                                cAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTcacc                                                                                                                                        >  2:119169/1‑288 (MQ=255)
                                                                                                                                                      cTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGgag                                                                                                                                 >  1:153338/1‑289 (MQ=255)
                                                                                                                                                       tGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGGGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGgaga                                                                                                                                >  1:259144/1‑289 (MQ=255)
                                                                                                                                                       tGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGgaga                                                                                                                                >  1:350398/1‑289 (MQ=255)
                                                                                                                                                       tGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGgaga                                                                                                                                >  2:168031/1‑289 (MQ=255)
                                                                                                                                                                  gTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTcc                                                                                                                      >  2:194980/1‑288 (MQ=255)
                                                                                                                                                                          gCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCAcc                                                                                                             <  2:153338/289‑1 (MQ=255)
                                                                                                                                                                            cccTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTa                                                                                                           >  2:368813/1‑289 (MQ=255)
                                                                                                                                                                            cccTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCAACGGACAGCTGGCGTCCGCCTTCAACTa                                                                                                           >  1:126466/1‑289 (MQ=255)
                                                                                                                                                                                gTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCa                                                                                                       >  1:40467/1‑289 (MQ=255)
                                                                                                                                                                                                       gTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTACCAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCCACTGTTGAGACCTGCCGCCTCGa                                                                                >  1:22459/1‑289 (MQ=255)
                                                                                                                                                                                                       gTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGa                                                                                >  2:400468/1‑288 (MQ=255)
                                                                                                                                                                                                       gTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGa                                                                                >  2:230019/1‑288 (MQ=255)
                                                                                                                                                                                                            ttGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTgcaag                                                                                                                                                                                                                                                             >  2:142813/1‑111 (MQ=255)
                                                                                                                                                                                                            ttGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTgcaag                                                                                                                                                                                                                                                             <  1:142813/111‑1 (MQ=255)
                                                                                                                                                                                                               cTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCttt                                                                        >  2:144480/1‑289 (MQ=255)
                                                                                                                                                                                                               cTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCttt                                                                        >  1:171673/1‑289 (MQ=255)
                                                                                                                                                                                                               cTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCCCGCGCAGCttt                                                                        >  1:56151/1‑289 (MQ=255)
                                                                                                                                                                                                                         tGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCaa                                                                                                      <  2:148127/249‑1 (MQ=255)
                                                                                                                                                                                                                         tGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCaa                                                                                                      >  1:148127/1‑249 (MQ=255)
                                                                                                                                                                                                                         tGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCa                                                               >  2:88493/1‑288 (MQ=255)
                                                                                                                                                                                                                                   ccATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCa                                                     <  1:180756/288‑1 (MQ=255)
                                                                                                                                                                                                                                   ccATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAg                                                   <  1:400468/289‑1 (MQ=255)
                                                                                                                                                                                                                                   ccATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAg                                                   <  1:230019/289‑1 (MQ=255)
                                                                                                                                                                                                                                             tattTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGAc                                          <  1:88493/289‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgttgGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACaga                                       >  2:427729/1‑288 (MQ=255)
                                                                                                                                                                                                                                                    tgGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACCGAGcc                                    >  2:430751/1‑288 (MQ=255)
                                                                                                                                                                                                                                                    tgGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGcc                                    >  1:342384/1‑288 (MQ=255)
                                                                                                                                                                                                                                                    tgGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGcc                                    >  2:17409/1‑288 (MQ=255)
                                                                                                                                                                                                                                                     ggAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCtc                                  <  2:350398/289‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaCGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCtcat                                <  1:386251/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCaggcgaggc             <  2:259144/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     aCCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGAt  >  1:202569/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                     aCCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGAt  >  2:14724/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               
TAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCGCATCATGTAACAGCGCACGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGAT  >  NC_000913/927333‑927898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: