Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 935,676 C→T 6.5% T818I (ACC→ATC)  ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913935,6760CT 6.5% 194.2 / 3.8 62T818I (ACC→ATC) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base T (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

TTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCTCGCCTGGTGGAAGCGCGTCTGGCTGATTTCCGTATTAAAGCCGATGTCGTCAATTACTCTCCGGGGC  >  NC_000913/935427‑935885
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                 
ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTCTCAGCAGCCCACACCGCAAGTTGCACCCACACCCCAGGATCCCCTCATTCCTCCGCAGTTGCTGCg                                                                                                                                                                                           >  4:359220/1‑274 (MQ=21)
      aaccgttACAGCAGCCGCAACAGCCGGTTGCACCGCAGCAGCAATAGCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGc                                                                                                                                                                                    <  1:136378/269‑1 (MQ=255)
                   gcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGCTACCCTGATTCATCCGCTGTTGATGCGTAATGGCGAAAGtcgtccgt                                                                                                                                                                       >  2:114282/1‑275 (MQ=255)
                   gcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccgt                                                                                                                                                                       >  2:333946/1‑275 (MQ=255)
                      gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGc                                                                                                                                                                    <  2:208964/275‑1 (MQ=255)
                                                                                      cagTATCAGCAGCCACAACCGCCGGTTGCGCCGCAGCCACAATAGCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                                                                                                     <  2:36179/274‑1 (MQ=255)
                                                                                       agTATCAGCAGCAACAACAGCCGGTTGCGCCGCAGCCACAAGATCAGCAGCCGCAACAACCTGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                                                                                                     <  2:437652/273‑1 (MQ=255)
                                                                                                 gcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCCTCCGCAGTTCATTCGTAATGGCTACCGCCGTCCCTTGCCTAAACCGGCGCCCCCGCTGCCTTCTCTGGATTTCCTTACGCCCCCGCCGGCCGACGTGGAGcccggcaatacct                                                                                         >  4:377337/1‑264 (MQ=25)
                                                                                                       acaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGc                                                                                                                                               <  1:89657/215‑1 (MQ=255)
                                                                                                               ttGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGGTATCCTGcggttg                                                                                                                                                                                                         <  7:182155/149‑6 (MQ=255)
                                                                                                                 gcgcCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACaa                                                                           <  2:377342/273‑1 (MQ=255)
                                                                                                                            acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                                                                                                                                                                                                           <  6:274144/134‑1 (MQ=255)
                                                                                                                            acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                                                                                                                                                                                                           >  5:274144/1‑134 (MQ=255)
                                                                                                                            acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGc                                                                                                                                                                                            >  7:150173/1‑149 (MQ=255)
                                                                                                                            acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGc                                                                                                                                                                                            >  5:226271/1‑149 (MQ=255)
                                                                                                                                     gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCg                                                                                                                                                                                   >  7:288918/1‑149 (MQ=255)
                                                                                                                                     gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCg                                                                                                                                                                                   >  7:214119/1‑149 (MQ=255)
                                                                                                                                        gcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCTCGCCTGGTGGAAGCGCGTc                                                  >  2:80572/1‑275 (MQ=255)
                                                                                                                                           gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGcc                                                                                                                                                                             >  7:204008/1‑149 (MQ=255)
                                                                                                                                              acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCTCGCCTGGTGGAAGCGCGTCTGGCTg                                            >  2:462831/1‑275 (MQ=255)
                                                                                                                                                aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg                                                                                                                                                                        >  7:400818/1‑149 (MQ=255)
                                                                                                                                                aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg                                                                                                                                                                        >  5:116505/1‑149 (MQ=255)
                                                                                                                                                 acCGGTTGCGCCGCAGCAGCAATATCAGCAGCCGCAACAACCCGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccgt                                                                                                                                                                       <  6:219792/149‑1 (MQ=255)
                                                                                                                                                  ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGtt                                                                                                                                                                      <  7:222062/149‑1 (MQ=255)
                                                                                                                                                     gTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa                                                                                                                                                                   <  7:293360/149‑1 (MQ=255)
                                                                                                                                                     gTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCTCGCCTGGTGGAAGCGCGTCTGGCTGATTTc                                       <  4:80571/273‑1 (MQ=255)
                                                                                                                                                          gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                                                                                              >  7:464220/1‑149 (MQ=255)
                                                                                                                                                          gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                                                                                              <  6:116505/149‑1 (MQ=255)
                                                                                                                                                                   gcaATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGAcgc                                                                                                                                                     <  6:429963/149‑1 (MQ=255)
                                                                                                                                                                     aaTATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCTCGCCTGGTGGAAGCGCGTCTGGCTGATTTCCGTATTAAAGCCGATgtc                     <  1:148175/275‑1 (MQ=255)
                                                                                                                                                                           agcagcCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCgacgac                                                                                                                                                       >  2:179529/1‑139 (MQ=255)
                                                                                                                                                                            gcagcCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                                                                                              >  7:303735/1‑131 (MQ=255)
                                                                                                                                                                            gcagcCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAAGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                                                                                              <  8:303735/131‑1 (MQ=255)
                                                                                                                                                                             cagcCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTg                                                                                                         >  2:516711/1‑183 (MQ=255)
                                                                                                                                                                               gcCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGAcgccg                                                                                                                                                   <  4:179527/139‑1 (MQ=255)
                                                                                                                                                                                ccGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTctct                                                                                                                                        >  5:267551/1‑149 (MQ=255)
                                                                                                                                                                                 cGCAACAACCGGTTGCGCCGCAGAANCAGTNCCAGCAGCCACAACAGCCAGTGGCCCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                                                                                                     <  4:516706/183‑1 (MQ=255)
                                                                                                                                                                                  gCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCCCAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTgg                                                                                                                                      <  8:400818/149‑1 (MQ=255)
                                                                                                                                                                                   caacaaCCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGa                                                                                                                                     <  6:137692/149‑1 (MQ=255)
                                                                                                                                                                                        acCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGc                                                                                                                                <  6:272345/149‑1 (MQ=255)
                                                                                                                                                                                        acCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGc                                                                                                                                <  6:203459/149‑1 (MQ=255)
                                                                                                                                                                                            gTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCTCGCCTGGTGGAAGCGCGTCTGGCTGATTTCCGTATTAAAGCCGATGTCGTCAATTACTCTCCGGGGc  <  2:359225/271‑1 (MQ=255)
                                                                                                                                                                                                      agcagcagTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACAccgccgccg                                                                                                                  >  6:88383/1‑149 (MQ=255)
                                                                                                                                                                                                       gcagcagTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGa                                                                                                                 <  8:150173/149‑1 (MQ=255)
                                                                                                                                                                                                          gcagTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCCCGCCTGGTGGAAGCGCGTCTGGCTGATTTCCGTATTaaa                              >  2:395762/1‑229 (MQ=255)
                                                                                                                                                                                                                  agcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAg                                                                                                      <  8:214119/149‑1 (MQ=255)
                                                                                                                                                                                                                   gcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                                                                                                     <  7:207767/149‑1 (MQ=255)
                                                                                                                                                                                                                   gcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                                                                                                     >  5:277260/1‑149 (MQ=255)
                                                                                                                                                                                                                   gcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                                                                                                     >  7:111029/1‑149 (MQ=255)
                                                                                                                                                                                                                      gcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCgg                                                                                                  <  7:1397/149‑1 (MQ=255)
                                                                                                                                                                                                                          caacaGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGa                                                                                              <  6:277260/149‑1 (MQ=255)
                                                                                                                                                                                                                             caGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATAc                                                                                           <  8:288918/149‑1 (MQ=255)
                                                                                                                                                                                                                               gCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGCGCCGGTAGATACCt                                                                                         <  6:267551/149‑1 (MQ=255)
                                                                                                                                                                                                                               gCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCt                                                                                         <  7:457641/149‑1 (MQ=255)
                                                                                                                                                                                                                                     tGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGcc                                                                                                                                                                             >  2:115154/1‑59 (MQ=255)
                                                                                                                                                                                                                                       cgcCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCtt                                                                                 <  8:464220/149‑1 (MQ=255)
                                                                                                                                                                                                                                        gcCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTg                                                                                >  8:472917/1‑149 (MQ=255)
                                                                                                                                                                                                                                         ccACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtcc                                                                                                                                                                         <  4:115153/59‑1 (MQ=255)
                                                                                                                                                                                                                                           acaACCGCAGGATATCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTctct                                                                                                                                        >  5:285819/1‑90 (MQ=255)
                                                                                                                                                                                                                                           acaACCGCAGGATATCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTctct                                                                                                                                        <  6:285819/90‑1 (MQ=255)
                                                                                                                                                                                                                                           acaACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAAc                                                                             >  5:305450/1‑149 (MQ=255)
                                                                                                                                                                                                                                             tgaagCAGGGTATCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACaa                                                                           >  5:477429/5‑149 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                 
TTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAAATGGCTCGCCTGGTGGAAGCGCGTCTGGCTGATTTCCGTATTAAAGCCGATGTCGTCAATTACTCTCCGGGGC  >  NC_000913/935427‑935885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑