Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,507,249 C→T D30N (GAC→AAC)  yhfT ← uncharacterized protein YhfT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,507,2490CT98.7% 238.9 / NA 76D30N (GAC→AAC) yhfTuncharacterized protein YhfT
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (23/51);  minor base G (0/1);  total (23/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCAATGGTGCAAC  >  NC_000913/3507003‑3507487
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ggCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGAGCGGGCGGATGCCGTTATGa                                                                                                                                                                                                                                            <  1:160094/251‑1 (MQ=255)
    attAAGGAGCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATAACGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGCCCAAAGTCGATACTCCCCGCCTCGCGACGTTTCATATAGCCGTCAAGCAGTTGCGGCAGGAACGGGCGGATGCCGTTATGaaaac                                                                                                                                                                                                                                        <  2:471125/249‑2 (MQ=255)
    agcaagaaTCAACACGCCCCAGATAGCGCAAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCTGTAGAAAGAAAGAGTAACCATGCTTTGAGCAGCCCGGTTTTCAGGGTAAAAAAGATCCCCACCGAGGCCTCAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGaaaaa                                                                                                                                                                                                                                        <  2:488160/244‑1 (MQ=255)
            acaccacGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCCCCAGGACGACGAGGATATCGGTAGGAAGAAAGAGAAACCAGGCGTTGAGCAGCCAGGTTTTCAGGGTAAAAGAGACCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAgccgcc                                                                                                                                                                                                                                <  2:606844/247‑1 (MQ=255)
             cAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGTGCGGATGCCGTTATGAAAAAAAGccgccg                                                                                                                                                                                                                               >  1:313321/1‑251 (MQ=255)
              aaCACCCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGAAGGCAGAACGAGAAACCATGCGTTGAGCAGCCCGGATTTCCGAGTACAAGAGATCCCCACAGCGGCCACAAAACCAATACTCCGACCAAAAGCGATACACCCAGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAgcagccc                                                                                                                                                                                                                               >  1:101751/1‑249 (MQ=255)
                cacGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAAAAGCCGCCGAgg                                                                                                                                                                                                                            >  2:558163/1‑250 (MQ=255)
                   gCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGAt                                                                                                                                                                                                                         <  1:548533/251‑1 (MQ=255)
                     cccAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGCCGCCGCGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTACAAGAGATCCCCACCGAGGACACCAAACCAATGCTCAGACCAAAAGCGCTACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGCTGCAGTTATGAAAAACAGCCGCCGAGCGATgc                                                                                                                                                                                                                       >  2:641398/1‑251 (MQ=255)
                            aGCGCCAGGACCAAACGCAATCAGGCTGTTTATCGCCGGGAGGCCGAGGATATCGGAGGGAAGAAAGAGTAGCCAGGGGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAAGGCGCGGACCAAAAGCGATACGCCCCGCCGCGCGACGGTGCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCTGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGaa                                                                                                                                                                                                                <  2:366425/251‑1 (MQ=255)
                            aGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGCTATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGa                                                                                                                                                                                                                 >  1:637943/1‑250 (MQ=255)
                                    gACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGTAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTc                                                                                                                                                                                                        >  2:135061/1‑251 (MQ=255)
                                      ccAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCAt                                                                                                                                                                                                      >  1:636950/1‑251 (MQ=255)
                                           cGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCg                                                                                                                                                                                                 <  1:18345/251‑1 (MQ=255)
                                              catcaagatagttATAGCCAGGACGCTTAGGATATAGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGCCGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGGGGATGCAGATATGACAAACAGACCCCGAGCGAGGAGCCAGAAGCGCTGTCATACCCGTc                                                                                                                                                                                               >  2:462712/12‑250 (MQ=255)
                                               aTCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAg                                                                                                                                                                                             <  1:466761/251‑1 (MQ=255)
                                                              gCCAGAACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCACCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCTTACCCGTCAGGCACGccaccccgct                                                                                                                                                                              >  2:652157/1‑246 (MQ=255)
                                                               ccAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATa                                                                                                                                                                             <  1:77044/251‑1 (MQ=255)
                                                                  gtaCGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGTTTTGCAGGGTAAAAGAGATCCCCACCGAGGCCACAAACCCAATGCTCAGCCCAAAAGCGATCCTCCCCGCCTCGGGACTGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATc                                                                                                                                                                          <  1:288953/249‑1 (MQ=255)
                                                                  gggcGCCGAGGATATCGGGAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTAGGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATc                                                                                                                                                                          <  1:358136/248‑1 (MQ=255)
                                                                   gaccgcGAGGGTATCTGTAGGAAGACGGAGTAGCCATGCGTTGACCAGCCCGGTTTCCAGGGTAAAAGATACCCCAGGCGAGGCATCAAAACCAATGCCCACACCAAAAGCGAGACTCGCCGCCTCGCGCCGGTTCATATAGCCTTGAACCAGTTGCGGCAGGTTCGGGCGGATGCCGTTATGAAAAACAGCCGCCGATCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCt                                                                                                                                                                         <  1:698941/246‑1 (MQ=255)
                                                                      gCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCCGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGa                                                                                                                                                                       >  2:228737/1‑250 (MQ=255)
                                                                        cGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGATCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAATCGATCCTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGAGCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAata                                                                                                                                                                    <  1:74985/251‑1 (MQ=255)
                                                                           ggATATCGGTAGGAAGAAGGAGTACCCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCCCCGAGGGCCAAAAGCCAATTCTCTGACCAAAAGCGTTACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGAGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACTGCACCACGATTATCTGAATATAc                                                                                                                                                                 <  2:65766/251‑1 (MQ=255)
                                                                            gATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCGATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCCGGCACGCCACCACGATAATCTGAATATAc                                                                                                                                                                 >  1:191619/1‑250 (MQ=255)
                                                                               atCGGAAGGAAGTAAGAGTAACCAGGCGTTGAGCAGCCCGGTTTTCAGTGAAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCATCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGa                                                                                                                                                              <  2:215437/250‑1 (MQ=255)
                                                                                 cGGTAGGAAGTAAGTGTAACCATGCGTTGAGCAGCCCAGTTTGCATGGTAAAAGAGATCCCCGCGTAGGCCACAAATTCATTGCACTGACCAAAAGCTATACTCCCCGCCTCTGTACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCTGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCTATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATc                                                                                                                                                            <  1:408172/250‑1 (MQ=255)
                                                                                   gTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCTCCACGATAATCTGAATATACAGATCCAt                                                                                                                                                         <  2:255746/251‑1 (MQ=255)
                                                                                     aGGAAGAAAGAGTCACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAg                                                                                                                                                       <  1:597693/251‑1 (MQ=255)
                                                                                       gaagaaAGAGTAACCATGGGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGACCCCCACCGAGGCCACAAAACCAATGCGCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAggc                                                                                                                                                     <  1:186806/251‑3 (MQ=255)
                                                                                            aagaGTAACCTTGCGTTGAGCAGCCCGGTTTTCAGGGGAAAAGAGAACCCCACCGAGGCCACAAAACCAATGCTAAGACCAAAAGCGATACTCCCCGGCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCTCGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGcc                                                                                                                                                <  1:620420/251‑1 (MQ=255)
                                                                                            aagaGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCCGGATCGGGCGGACGCCGTTATGAAAAACAGCCGCCGAGCGATGCCCCAGAAGCGATGGCATACCCGACAGGCACGCCACCACGATACTCTGAATATACAGATCCATAGTCGTGcc                                                                                                                                                >  2:190640/1‑251 (MQ=255)
                                                                                                  aaCCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCCCCGAGGCAACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGAGCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATTCAGATCCATAGTCGTGCCTTTTa                                                                                                                                           <  1:190640/250‑1 (MQ=255)
                                                                                                  aaCCAGGCGTGGAGCAGCCCGTTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAATGCGATCCGCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGAGCCATAGTCGTGCCTTTTa                                                                                                                                           <  2:11546/250‑1 (MQ=255)
                                                                                                            tGAGCAGCCCGGGTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCGTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGAGGTCATACCCGTCAGGCACGCCACCACTATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATgg                                                                                                                                 <  1:462712/250‑1 (MQ=255)
                                                                                                                cagcCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACACAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGACTTTTATTTCAGATGGTTAAt                                                                                                                            >  2:432914/1‑251 (MQ=255)
                                                                                                                cagcCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGa                                                                                                                                                                                                                                            >  2:277505/1‑139 (MQ=255)
                                                                                                                cagcCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGa                                                                                                                                                                                                                                            <  1:277505/139‑1 (MQ=255)
                                                                                                                   cccGGTTTTCAGGGTAAAAGAGATCCCCCCCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAg                                                                                                                         <  2:367362/251‑1 (MQ=255)
                                                                                                                       gTTTTAAGGGTAAAAGAAATCCCCACCGGTGCCACAAATCCAATGCTCAGTCCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGTCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCAta                                                                                                                     <  1:490900/251‑2 (MQ=255)
                                                                                                                          ttCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGAGGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCgg                                                                                                                   <  1:323399/250‑1 (MQ=255)
                                                                                                                             aGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATc                                                                                                               <  1:46406/251‑1 (MQ=255)
                                                                                                                             aGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTGATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATc                                                                                                               <  1:658312/251‑1 (MQ=255)
                                                                                                                               ggTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATcgc                                                                                                             <  1:21521/251‑1 (MQ=255)
                                                                                                                                  aaaaGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCgtg                                                                                                          <  1:176291/251‑1 (MQ=255)
                                                                                                                                  aaaaGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCgtg                                                                                                          >  1:441751/1‑251 (MQ=255)
                                                                                                                                            cccACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCCCGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGc                                                                                                >  2:660571/1‑251 (MQ=255)
                                                                                                                                               cccGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCa                                                                                             <  1:185620/250‑1 (MQ=255)
                                                                                                                                               aCCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTc                                                                                              >  1:52610/1‑250 (MQ=255)
                                                                                                                                                         cacaAAACCAATGCTCTGACCAAATGGTGTACTCCCCGCCTCGCGACGGTTCATCTAGCCTTCAATCAGTTGCGGCAGGAGCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAAGCAGCCTCGGAAACGCGTGTTCGACGTTCTCAACGGAAAg                                                                                     <  2:74085/251‑1 (MQ=255)
                                                                                                                                                             aaaCCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCa                                                                                 <  1:432914/251‑1 (MQ=255)
                                                                                                                                                                   aTGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGGGCTCAACGGAAAGGCCAAACGcc                                                                           >  2:492919/1‑251 (MQ=255)
                                                                                                                                                                    ggCTCAGACCAAAAGCTATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATATGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAAGCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCa                                                                          <  2:637943/250‑1 (MQ=255)
                                                                                                                                                                       ccAGACCAAAAGCGATACCCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTGGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACtt                                                                       <  2:441751/250‑1 (MQ=255)
                                                                                                                                                                         agcccAAAAGCGATCCCCCCCGCCTCGCGACTGTTCAGATAGCCTTCAATCAGTTGCGGCATGATCGGGCGGATGCCGGTATGAAAAACAGCCTCCGAGCGATGCGCCAGAAGCGATGTCATTCCCGTCAGGCACGCCTCCACGATAATCGGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGGGCTCAATGGAAAGGCCAAACGCCTCtttt                                                                     <  2:317032/248‑1 (MQ=255)
                                                                                                                                                                            ccAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCGCGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCCCTTTTCCt                                                                  <  1:26838/251‑1 (MQ=255)
                                                                                                                                                                             cAAAAGCGATTCTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCtt                                                                 <  2:333105/251‑1 (MQ=255)
                                                                                                                                                                                  gCGATACCCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCATTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCg                                                            <  1:541686/251‑1 (MQ=255)
                                                                                                                                                                                    gatcccccccGCCGCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCTTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAGCGCCACTTTGCCTTCAGCGAt                                                          <  2:636950/245‑1 (MQ=255)
                                                                                                                                                                                     aTACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGTCAGGAACGGGCGGATGCCGTTATGAACAACAGCCGCCGAGCGATGCGCCAGAAACGATGTCCTACCCGTCAGGCACGCCACCACGCTAATCTGCATATACAGATCCATCGTCGTGCCTTTTATTTCAGATGGTTAATCCGCATCGGAATCCCGTTTTAGACGTGCTCAACGGAAAGGCCAACCGCCACTTTTCCTTCCGCGATc                                                         >  2:405253/1‑251 (MQ=255)
                                                                                                                                                                                         cccccGCCGCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGAGGGTTAATCAGCAGCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCAttt                                                     <  1:641398/250‑1 (MQ=255)
                                                                                                                                                                                           cccGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTg                                                   <  1:659681/251‑1 (MQ=255)
                                                                                                                                                                                                       cGGTTCATATAGCCTTCAAACAGTTGCGTCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCCAGCGATGCGCCAGAAGCGATGTCATTCCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTcc                                       <  2:27410/251‑2 (MQ=255)
                                                                                                                                                                                                        ggTTCATATAGCCTTCAATCAGTTGCGGCAGGATGGGGCGGACGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATGCCCGTCAGGCACGCCGCCACGATAAACTGAATAGACAGATCCATAGTCGTGCCTTTTATTTCATATGGTTAATCAGCATCGGAATCTCGTTGTCGACTTGCTCAACGGAGAGGCCAAACTCCTCTGTTCCGTCGGCGGTCATTGTGGCGATATGCTCGt                                      <  2:363191/251‑1 (MQ=255)
                                                                                                                                                                                                              tataGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCg                                >  1:126784/1‑251 (MQ=255)
                                                                                                                                                                                                               ataGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCTCTTTTCCGGCAGCGATCATTTTGGCGATATGCTCGTCTTTCg                                <  1:208133/250‑1 (MQ=255)
                                                                                                                                                                                                               ataGCCGTCAATCAGTTTCGGCAGGATCGGTCGGATGCCGTTATGAAAAACAGCCTCCGAGCGATGCGCCAGAAGCGATGTCAGACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGGGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGc                               <  1:33041/251‑1 (MQ=255)
                                                                                                                                                                                                                taGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCt                              >  1:593092/1‑251 (MQ=255)
                                                                                                                                                                                                                   ccTTCAATCAGTGGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTa                           <  1:558163/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGCAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGc             >  2:150494/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ggCAGGAGCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCGGGTTTGGCa            <  1:150494/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gccgtaTCGGGCGGATTCCGTTATGAAAAACAGCCGCCGATCGATGCGCCAGAAGCGATGTCATTCCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTTCCTTTTATTTCAGATGTTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCaa           <  1:225695/246‑1 (MQ=255)
                                                                                                                                                                                                                                   gCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCaa           <  2:75445/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGGTGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCAATGGTGCAAc  <  1:228737/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCCCGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTc                                                                                              >  2:649063/1‑154 (MQ=255)
                                                                                                                                                                                                                                                 aTGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTc                                                                                              <  1:649063/154‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
GGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCAATGGTGCAAC  >  NC_000913/3507003‑3507487

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