Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,240 G→A L21L (TTG→TTA creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,2400GA96.4% 186.9 / NA 56L21L (TTG→TTAcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/27);  minor base C (2/0);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGC  >  NC_000913/4637997‑4638486
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
cagcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAACTGACTAGCCTGATTGGTGCAGTATTGCTGTTGTTAATTCCg                                                                                                                                                                                                                                                   >  1:381386/1‑250 (MQ=255)
              tGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCg                                                                                                                                                                                                                                     <  1:419055/250‑1 (MQ=255)
                 cGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAg                                                                                                                                                                                                                                 <  2:170721/251‑1 (MQ=255)
                        cAGTGAAGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGagtg                                                                                                                                                                                                                          >  2:462097/1‑248 (MQ=255)
                             aGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa                                                                                                                                                                                                                      >  1:196316/1‑250 (MQ=255)
                                      gtttgtttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGAt                                                                                                                                                                                                           >  1:547545/1‑251 (MQ=255)
                                      gtttgtttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGAt                                                                                                                                                                                                           >  2:177907/1‑251 (MQ=255)
                                                   cGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGAtg                                                                                                                                                                                              <  1:161915/251‑1 (MQ=255)
                                                            gttGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTTAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATGGCTGTTGTTAATGCCTAGAATGCTg                                                                                                                                                                                                                                          <  2:46281/197‑1 (MQ=255)
                                                            gttGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTg                                                                                                                                                                                                                                          >  1:46281/3‑199 (MQ=255)
                                                            gttGCGCAACATGCCGGCAGGTGGCGCACTGGCCTCGATTCGCCTTCACCGTCCCTTCACATAGCTTCAACTTCTTCCCACATAGTCTTCACATCCTGCTGACATTGCACAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATg                                                                                                                                                                                     >  2:247505/3‑251 (MQ=255)
                                                                         cAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAg                                                                                                                                                                        <  1:569955/251‑1 (MQ=255)
                                                                           ggAAGGTGGCGTGCGGGCCTCGCTTCCATTTCCCCGTCACTTCACATAGATTCAAATTCTTCCCACATAGGTTTCTTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGTAGTATTGCTGTTGTTAATTCCGATAATGCTGACTCGTCAGGTGATTGTCGAACGTGCTGATGACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTa                                                                                                                                                                      <  2:404009/251‑1 (MQ=255)
                                                                                  ggCGTGCTTGCCTCGCTTCGCCTTCACCGTCACTTCACATAGCTTCAAATTCTCCCCACATAGTACTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCg                                                                                                                                                                <  2:148973/250‑1 (MQ=255)
                                                                                    cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGc                                                                                                                                                             <  1:222604/251‑1 (MQ=255)
                                                                                     gTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACAGAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTTGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCGGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCggggc                                                                                                                                                            <  2:93120/251‑3 (MQ=255)
                                                                                          ggCCTCGCTTCGACTTCACCGTCACTTCACATAGCGTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCGGATTACCGTAGCGATGTGGAAGAGGCGAGTCGCCAGAGGACCAGCGTGCCGCa                                                                                                                                                         <  2:272302/249‑1 (MQ=255)
                                                                                                     gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCACCTCCTGCTGCCATTGCCACGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAACATCCCTAGCCTGATAGGTGCAGTATTGCTGTTGTTCATTCCGTCAATGCTGACTAGGCAGGTGATTGTCGCACGTGCTGCATTCAGTAGCCCTGTAGATGATCCGATTCGCAAAAGTGCCCGCGGGAAGCAAAAACTCGTTg                                                                                                                                              >  1:542245/1‑249 (MQ=255)
                                                                                                             cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGTAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGGTGGGCCGCTCAt                                                                                                                                    >  2:4873/1‑251 (MQ=255)
                                                                                                              gTCACTGCCCATAGCTTCAAATTCTTCCCACATAGTCTTCCTATCCTGCTGCCATTGCGAAGGAGAAGTCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTGTGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATGCGGCAGGTGATTGTCCAACGTGCTGACTACCGTAGCTATTTGGATGATGCGATTCTCCAAAGTAGCAGCGGGCCGCAAAAGCTCggttgtccgttcatc                                                                                                                                   <  2:95292/251‑15 (MQ=255)
                                                                                                                    tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTAtt                                                                                                                             >  1:462282/1‑251 (MQ=255)
                                                                                                                      acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTcc                                                                                                                           <  1:158408/251‑1 (MQ=255)
                                                                                                                          aGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCtgt                                                                                                                        <  2:191104/250‑1 (MQ=255)
                                                                                                                              tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGAc                                                                                                                     <  1:164003/249‑1 (MQ=255)
                                                                                                                                aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTTGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAg                                                                                                                 <  1:392004/251‑1 (MQ=255)
                                                                                                                                      ttcCCACATAGTCTTCGTATCATGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGGTTGGCCGCTCATCGATATTCCTGTGACCGAGCTTtat                                                                                                           >  2:185788/1‑251 (MQ=255)
                                                                                                                                       tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCttt                                                                                                             <  1:30783/248‑1 (MQ=255)
                                                                                                                                       gcccACATAGTCTTCTGCTCCTGCTGCCATTGCGAAGGAGAAGACTATGTTGATATCCTCCCTGTGCTGTAAAATGACTAGCATGTTTGGTGCAATATTGCTGTTGTTAATTCCGATAATGCTGATTCTGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGAGGCGATTCGCACAAGTACCAGCGTGCCGCAATAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTta                                                                                                            <  2:146086/248‑1 (MQ=255)
                                                                                                                                          cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGACAAGCCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCAGCTCATCGCTATTCCTGTGACCGAGCTTTATACgg                                                                                                       >  2:41358/1‑251 (MQ=255)
                                                                                                                                           acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAACACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGt                                                                                                      >  2:278454/1‑251 (MQ=255)
                                                                                                                                                      gTATCCTGCTGCCATTGCATAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGTTGTGGaa                                                                                                                                                                                         >  2:492709/1‑157 (MQ=255)
                                                                                                                                                      gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGaa                                                                                                                                                                                         <  1:492709/157‑1 (MQ=255)
                                                                                                                                                        agccTGCTGTCATTGCAAAGGAGAATACTATGTTGAAATCCCCCCTGTTATGGAAAATGACTAGCCTGTTTGGTTCAGTATTGCTGTTGTTAATTCCGATAATGCGGATTCGGCAGGTGATTGGCGAACGTGCTGATTACCGTAGCGATGGGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCTTTGGGCCGCTCAGCGCTATTCCTGTGACCGAGCTTTATCCGGTGCAGTAAGAGGAt                                                                                         <  1:32981/249‑1 (MQ=255)
                                                                                                                                                           ctgccgacATTGCAAACGAGAACACTATGTTGCAATCCCCCCCGTTCTGGAAACTGACTAGCCTGTATGATGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGCTGTGGAAGATGCGATTTGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATACCTATTACTGTGACCGAGCTTTCCACGGTGCAGGCAGACGATaa                                                                                       >  2:102820/8‑250 (MQ=255)
                                                                                                                                                             gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGTGGATTGTCGAACGTGCTGATTACCGGAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCTCTCAGCGCTATTCCTGTGACCGAGCTTTATTCGGTGCAGGAAGATGATaaaa                                                                                     <  2:524893/250‑1 (MQ=255)
                                                                                                                                                             gacgtcATTGCAAAGGAGAAGACTATGTTGAAATCGCCCCGGTTCTGGAAAATGTCTAGCCTGTTTGTGGCAGGATGGCTGTTGTTAATTCCGATAATGCTGATTCGTCAGGTGATTTTGGAACGTGCGGATTACCGTAGTGATTTGGAAGGTGCTATTCGTCCAAGTGCCAGCGGGCCTGAAAAGCGTTTTGGGACTCCCATCTCTATTCCTGTGACCCTGCTTTATACGTTGCAGTATGAGGATaaaa                                                                                     <  2:278548/245‑1 (MQ=255)
                                                                                                                                                               tgcCATAGCAAAGGAGAAGACTATGTTGAAATACCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGTAGAGGATAAAACCg                                                                                  >  1:483077/1‑251 (MQ=255)
                                                                                                                                                                  cATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGCACGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                                               >  1:11802/1‑251 (MQ=255)
                                                                                                                                                                      gCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAATCTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGatcg                                                                           <  2:453562/251‑4 (MQ=255)
                                                                                                                                                                        aaaGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCTCCAAAGTACCAGCGGGACGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCgg                                                                          >  1:431358/1‑250 (MQ=255)
                                                                                                                                                                          aGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGaaa                                                                       <  1:330159/251‑1 (MQ=255)
                                                                                                                                                                             agaagcCTATGATGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGAATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGACGCTCATAGCTATTACTGTGACCGGGCTTTCTACGGTGCAGGAAGAGGATAAAACCGTGGCGCGGCAACg                                                                     >  2:300674/1‑250 (MQ=255)
                                                                                                                                                                                agaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAAGt                                                                 <  1:214760/251‑1 (MQ=255)
                                                                                                                                                                                           aaaTCCCCCCTGTTCTGGAACATGACTAGACTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAATCTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATCCGGTGCCGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCAttt                                                      >  1:432632/1‑251 (MQ=255)
                                                                                                                                                                                               ccccccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGATGTTGTTCATTCCGTTATTGCTGCTTGGACAGGTGCTTGGAGGGCGTGCTGCTTGCCGCAGCGATGAGGTAGATCCGCTTCGCAAAAGTACCAGCGGGCCGCTCAAGCTCGATGCGCCCCTCATCGCTCTTCCTGTGAACGAGATTTATACGGCGACGGAAGAGGATAAACCCGTGAACCGGCACCGAAGATTTATCCAttttttg                                                  >  1:585197/1‑249 (MQ=255)
                                                                                                                                                                                                  cccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGtt                                                >  1:340153/1‑250 (MQ=255)
                                                                                                                                                                                                         cTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAACAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTgcgt                                        >  1:197144/1‑248 (MQ=255)
                                                                                                                                                                                                         cTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGt                                        >  1:209481/1‑251 (MQ=255)
                                                                                                                                                                                                          tGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGCAGTTTTGTCCATTTTTGGTTACCTGAGTc                                       >  1:138840/1‑251 (MQ=255)
                                                                                                                                                                                                            gAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGGGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGGTTTCTCCATTTTTGGTTACCTGAGTCAt                                     >  2:449034/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                          <  1:19057/146‑1 (MQ=255)
                                                                                                                                                                                                                                gTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                          >  2:19057/1‑146 (MQ=255)
                                                                                                                                                                                                                                  gCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGaa               <  1:278454/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTg           >  1:40228/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGga         >  1:469336/1‑249 (MQ=255)
                                                                                                                                                                                                                                             gttgttAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTTCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGGGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTGACCTGAGTCATTGATGTGTGATGTCAATCAGAACGTGgaagaa     <  2:547545/250‑1 (MQ=255)
                                                                                                                                                                                                                                               tgttaATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGc  >  2:432976/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
CAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGC  >  NC_000913/4637997‑4638486

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