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→TTA) creDputative 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

CAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGC  >  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)
                                      gtttgtttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑At                                                                                                                                                                                                           >  1:547545/1‑251 (MQ=255)
                                      gtttgtttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑At                                                                                                                                                                                                           >  2:177907/1‑251 (MQ=255)
                                                   cGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGAtg                                                                                                                                                                                              <  1:161915/251‑1 (MQ=255)
                                                            gttGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTTAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATGGCTGTTGTTAATGCCTAGAATGCTg                                                                                                                                                                                                                                          <  2:46281/197‑1 (MQ=255)
                                                            gttGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTg                                                                                                                                                                                                                                          >  1:46281/3‑199 (MQ=255)
                                                            gttGCGCAACATGCCGGCAGGTGGCGCACTGGCCTCGATTCGCCTTCACCGTCCCTTCACATAGCTTCAACTTCTTCCCACATAGTCTTCACATCCTGCTGACATTGCACAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATg                                                                                                                                                                                     >  2:247505/3‑251 (MQ=255)
                                                                         cAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAg                                                                                                                                                                        <  1:569955/251‑1 (MQ=255)
                                                                           ggAAGGTGGCGTGCGGGCCTCGCTTCCATTTCCCCGTCACTTCACATAGATTCAAATTCTTCCCACATAGGTTTCTTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGTAGTATTGCTGTTGTTAATTCCGATAATGCTGACTCGTCAGGTGATTGTCGAACGTGCTG‑ATGACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTa                                                                                                                                                                      <  2:404009/251‑1 (MQ=255)
                                                                                  ggCGTGCTTGCCTCGCTTCGCCTTCACCGTCACTTCACATAGCTTCAAATTCTCCCCACATAGTACTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCg                                                                                                                                                                <  2:148973/250‑1 (MQ=255)
                                                                                    cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGc                                                                                                                                                             <  1:222604/251‑1 (MQ=255)
                                                                                     gTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACAGAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTTGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCGG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCggggc                                                                                                                                                            <  2:93120/251‑3 (MQ=255)
                                                                                          ggCCTCGCTTCGACTTCACCGTCACTTCACATAGCGTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCGG‑ATTACCGTAGCGATGTGGAAGAGGCGAGTCGCCAGAGGACCAGCGTGCCGCa                                                                                                                                                         <  2:272302/249‑1 (MQ=255)
                                                                                                     gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCACCTCCTGCTGCCATTGCCACGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAACATCCCTAGCCTGATAGGTGCAGTATTGCTGTTGTTCATTCCGTCAATGCTGACTAGGCAGGTGATTGTCGCACGTGCTGCATT‑CAGTAGCCCTGTAGATGATCCGATTCGCAAAAGTGCCCGCGGGAAGCAAAAACTCGTTg                                                                                                                                              >  1:542245/1‑249 (MQ=255)
                                                                                                             cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGTAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGGTGGGCCGCTCAt                                                                                                                                    >  2:4873/1‑251 (MQ=255)
                                                                                                              gTCACTGCCCATAGCTTCAAATTCTTCCCACATAGTCTTCCTATCCTGCTGCCATTGCGAAGGAGAAGTCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTGTGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATGCGGCAGGTGATTGTCCAACGTGCTG‑ACTACCGTAGCTATTTGGATGATGCGATTCTCCAAAGTAGCAGCGGGCCGCAAAAGCTCggttgtccgttcatc                                                                                                                                   <  2:95292/251‑15 (MQ=255)
                                                                                                                    tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTAtt                                                                                                                             >  1:462282/1‑251 (MQ=255)
                                                                                                                      acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTcc                                                                                                                           <  1:158408/251‑1 (MQ=255)
                                                                                                                          aGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCtgt                                                                                                                        <  2:191104/250‑1 (MQ=255)
                                                                                                                              tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGAc                                                                                                                     <  1:164003/249‑1 (MQ=255)
                                                                                                                                aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTTGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAg                                                                                                                 <  1:392004/251‑1 (MQ=255)
                                                                                                                                      ttcCCACATAGTCTTCGTATCATGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGGTTGGCCGCTCATCGATATTCCTGTGACCGAGCTTtat                                                                                                           >  2:185788/1‑251 (MQ=255)
                                                                                                                                       tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCttt                                                                                                             <  1:30783/248‑1 (MQ=255)
                                                                                                                                       gcccACATAGTCTTCTGCTCCTGCTGCCATTGCGAAGGAGAAGACTATGTTGATATCCTCCCTGTGCTGTAAAATGACTAGCATGTTTGGTGCAATATTGCTGTTGTTAATTCCGATAATGCTGATTCTGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGAGGCGATTCGCACAAGTACCAGCGTGCCGCAATAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTta                                                                                                            <  2:146086/248‑1 (MQ=255)
                                                                                                                                          cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGACAAGCCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCAGCTCATCGCTATTCCTGTGACCGAGCTTTATACgg                                                                                                       >  2:41358/1‑251 (MQ=255)
                                                                                                                                           acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAACACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGt                                                                                                      >  2:278454/1‑251 (MQ=255)
                                                                                                                                                      gTATCCTGCTGCCATTGCATAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGTTGTGGaa                                                                                                                                                                                         >  2:492709/1‑157 (MQ=255)
                                                                                                                                                      gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGaa                                                                                                                                                                                         <  1:492709/157‑1 (MQ=255)
                                                                                                                                                        agccTGCTGTCATTGCAAAGGAGAATACTATGTTGAAATCCCCCCTGTTATGGAAAATGACTAGCCTGTTTGGTTCAGTATTGCTGTTGTTAATTCCGATAATGCGGATTCGGCAGGTGATTGGCGAACGTGCTG‑ATTACCGTAGCGATGGGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCTTTGGGCCGCTCAGCGCTATTCCTGTGACCGAGCTTTATCCGGTGCAGTAAGAGGAt                                                                                         <  1:32981/249‑1 (MQ=255)
                                                                                                                                                           ctgccgacATTGCAAACGAGAACACTATGTTGCAATCCCCCCCGTTCTGGAAACTGACTAGCCTGTATGATGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGCTGTGGAAGATGCGATTTGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATACCTATTACTGTGACCGAGCTTTCCACGGTGCAGGCAGACGATaa                                                                                       >  2:102820/8‑250 (MQ=255)
                                                                                                                                                             gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGTGGATTGTCGAACGTGCTG‑ATTACCGGAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCTCTCAGCGCTATTCCTGTGACCGAGCTTTATTCGGTGCAGGAAGATGATaaaa                                                                                     <  2:524893/250‑1 (MQ=255)
                                                                                                                                                             gacgtcATTGCAAAGGAGAAGACTATGTTGAAATCGCCCCGGTTCTGGAAAATGTCTAGCCTGTTTGTGGCAGGATGGCTGTTGTTAATTCCGATAATGCTGATTCGTCAGGTGATTTTGGAACGTGCGG‑ATTACCGTAGTGATTTGGAAGGTGCTATTCGTCCAAGTGCCAGCGGGCCTGAAAAGCGTTTTGGGACTCCCATCTCTATTCCTGTGACCCTGCTTTATACGTTGCAGTATGAGGATaaaa                                                                                     <  2:278548/245‑1 (MQ=255)
                                                                                                                                                               tgcCATAGCAAAGGAGAAGACTATGTTGAAATACCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGTAGAGGATAAAACCg                                                                                  >  1:483077/1‑251 (MQ=255)
                                                                                                                                                                  cATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGCACGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                                               >  1:11802/1‑251 (MQ=255)
                                                                                                                                                                      gCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAATCTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGatcg                                                                           <  2:453562/251‑4 (MQ=255)
                                                                                                                                                                        aaaGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCTCCAAAGTACCAGCGGGACGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCgg                                                                          >  1:431358/1‑250 (MQ=255)
                                                                                                                                                                          aGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGaaa                                                                       <  1:330159/251‑1 (MQ=255)
                                                                                                                                                                             agaagcCTATGATGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGAATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGACGCTCATAGCTATTACTGTGACCGGGCTTTCTACGGTGCAGGAAGAGGATAAAACCGTGGCGCGGCAACg                                                                     >  2:300674/1‑250 (MQ=255)
                                                                                                                                                                                agaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAAGt                                                                 <  1:214760/251‑1 (MQ=255)
                                                                                                                                                                                           aaaTCCCCCCTGTTCTGGAACATGACTAGACTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAATCTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATCCGGTGCCGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCAttt                                                      >  1:432632/1‑251 (MQ=255)
                                                                                                                                                                                               ccccccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGATGTTGTTCATTCCGTTATTGCTGCTTGGACAGGTGCTTGGAGGGCGTGCTG‑CTTGCCGCAGCGATGAGGTAGATCCGCTTCGCAAAAGTACCAGCGGGCCGCTCAAGCTCGATGCGCCCCTCATCGCTCTTCCTGTGAACGAGATTTATACGGCGACGGAAGAGGATAAACCCGTGAACCGGCACCGAAGATTTATCCAttttttg                                                  >  1:585197/1‑249 (MQ=255)
                                                                                                                                                                                                  cccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGtt                                                >  1:340153/1‑250 (MQ=255)
                                                                                                                                                                                                         cTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAACAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTgcgt                                        >  1:197144/1‑248 (MQ=255)
                                                                                                                                                                                                         cTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGt                                        >  1:209481/1‑251 (MQ=255)
                                                                                                                                                                                                          tGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGCAGTTTTGTCCATTTTTGGTTACCTGAGTc                                       >  1:138840/1‑251 (MQ=255)
                                                                                                                                                                                                            gAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGGGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGGTTTCTCCATTTTTGGTTACCTGAGTCAt                                     >  2:449034/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                          <  1:19057/146‑1 (MQ=255)
                                                                                                                                                                                                                                gTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                          >  2:19057/1‑146 (MQ=255)
                                                                                                                                                                                                                                  gCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGaa               <  1:278454/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTg           >  1:40228/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGga         >  1:469336/1‑249 (MQ=255)
                                                                                                                                                                                                                                             gttgttAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTTCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGGGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTGACCTGAGTCATTGATGTGTGATGTCAATCAGAACGTGgaagaa     <  2:547545/250‑1 (MQ=255)
                                                                                                                                                                                                                                               tgttaATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGc  >  2:432976/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
CAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTG‑ATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGC  >  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: ‑