Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,603 C→T 8.9% intergenic (+22/+16) ygiC → / ← ygiD ATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,6030CT 8.9% 159.3 / 5.3 56intergenic (+22/+16)ygiC/ygiDATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base T (2/3);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

CAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCGCCCACGGATACGGTGAACTATCACCGTGCCACTTCACTGTGCGCAGGT  >  NC_000913/3181303‑3181902
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           
aaaaTGGAAAATGTTTTGGTTGTACCTATTGTCTCCCGTGAATGCGCATACGTGTCGATCATTGAGAACGGCAAAACCAGTGAAGANGCGGAAGGTCTNNNGTGCGAAGAAGGGATNNNTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCTACAGCTATATGCTGATTGGGAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTc                                                                                                                                                                                                                                                                                                             <  1:451639/300‑1 (MQ=255)
        aaaaTATGTTGTTAAACCGATCTTTTCCCTGGAAGGGGCAAACTGGTCGATCATTAAGAATGGCAAATCCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTGCAGCAATTGCACCCGTTACCGAAATGCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATGCGTGAAGACCGTGCATTGATCACCCAGGATATGGCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCg                                                                                                                                                                                                                                                                                                     <  1:362287/301‑1 (MQ=255)
            tatGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGTat                                                                                                                                                                                                                                                                                                 <  2:267580/301‑2 (MQ=255)
                     aaaCCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGACGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACAGACAGGCACTCGCCATCCGGTAATTGTTAGcc                                                                                                                                                                                                                                                                                        >  1:431655/1‑301 (MQ=255)
                                           gcgcAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAACCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAACGACAACATACTc                                                                                                                                                                                                                                                                  >  1:198913/1‑301 (MQ=255)
                                                   gtgtCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGcc                                                                                                                                                                                                                                                          <  2:549660/301‑1 (MQ=255)
                                                                agaACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAAGCTGCCCATTTCAATAccc                                                                                                                                                                                                                                             >  1:67377/1‑301 (MQ=255)
                                                                 gaACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCAAACTCAGGCTGCCCATTTCTAAAccac                                                                                                                                                                                                                                            >  2:130958/1‑299 (MQ=255)
                                                                             cATTGATGCAGCTGAAGGTCGGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATgg                                                                                                                                                                                                                                <  1:71503/301‑1 (MQ=255)
                                                                             cATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCGGGCTGCCCATTTCTATACCCTCAACCGGAATgg                                                                                                                                                                                                                                <  2:432219/301‑1 (MQ=255)
                                                                              aTTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGCGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAACggg                                                                                                                                                                                                                               >  1:156178/1‑300 (MQ=255)
                                                                               ttGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTa                                                                                                                                                                                                                              <  2:359579/301‑1 (MQ=255)
                                                                                   ggcggcGGAAGGTCCGTATGGGGAAGAAGGGATGATTGTTCATCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATTCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                                                                                                                                                                          <  1:556741/297‑1 (MQ=255)
                                                                                     cagcGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                                                                                                                                                                                        <  2:67467/301‑1 (MQ=255)
                                                                                             ggTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAACCCCAACCGGAATGGTAATTTGACCCTGccc                                                                                                                                                                                                                >  1:271838/1‑301 (MQ=255)
                                                                                              gTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAAAGGTAATTGGCACCTGcccc                                                                                                                                                                                                               >  1:190145/1‑300 (MQ=255)
                                                                                                 cGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGATAATTGGCTCCTGCCCAAcc                                                                                                                                                                                                            >  2:74192/1‑301 (MQ=255)
                                                                                                     tGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACg                                                                                                                                                                                                        <  2:432338/301‑1 (MQ=255)
                                                                                                                       tGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATacaaca                                                                                                                                                                                      <  2:64486/301‑1 (MQ=255)
                                                                                                                           aagcgagttcccccGTTACCTATGTTGTGGGACAGCTATATTCTGTTTTGTAGCTTGCTGGTGAAAGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACgg                                                                                                                                                                                  <  1:130833/290‑1 (MQ=255)
                                                                                                                            aGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGc                                                                                                                                                                                 >  1:312410/1‑301 (MQ=255)
                                                                                                                                   ccACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCAACGCACCTAACACAACCAACAACGGCAaataaa                                                                                                                                                                          >  2:110034/1‑300 (MQ=255)
                                                                                                                                               gAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTg                                                                                                                                                              <  2:359584/301‑1 (MQ=255)
                                                                                                                                                  tttCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGt                                                                                                                                                           <  2:272141/300‑1 (MQ=255)
                                                                                                                                                   ttCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGtt                                                                                                                                                          <  2:199156/301‑1 (MQ=255)
                                                                                                                                                      gcttttccgattttGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCTTTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTATTTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGa                                                                                                                                                       <  1:119054/288‑1 (MQ=255)
                                                                                                                                                           cAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACg                                                                                                                                                  >  2:599680/1‑301 (MQ=255)
                                                                                                                                                             gttattTGGTGATTGGTAGTTGGTTGGTGAATGTTTAACTCGCATGTATTGGCATTCGTGAAGACCGTGCATTGATCACCCATGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTg                                                                                                                                                <  2:337745/299‑1 (MQ=255)
                                                                                                                                                                   tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCTTTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTATTTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGaggatcggaagagc                                                                                                                                          >  2:119178/1‑288 (MQ=255)
                                                                                                                                                                              aGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGATGGGTTCGATAACGTGCAAACTTCAAGGTCGAg                                                                                                                               >  1:174639/1‑301 (MQ=255)
                                                                                                                                                                                 tggcttgTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTa                                                                                                                            >  2:460949/1‑301 (MQ=255)
                                                                                                                                                                                         gAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCGGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATACTGCTCAGGCGCTTAGCACTATAACGTTCAACCTACATgccccagcacatacaccc                                                                                                                    >  1:598625/1‑284 (MQ=255)
                                                                                                                                                                                           aggATCAACCCGCCGGAATTGGCTTTCGTGAAGACTGTGCATTGATCACCCATGATATGTCTCGGTTTTATCCACAGATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCaga                                                                                                                  <  1:109910/299‑1 (MQ=255)
                                                                                                                                                                                                aaaCCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAGCGGCAAATAATGCGCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGAtg                                                                                                             <  1:460385/300‑1 (MQ=255)
                                                                                                                                                                                                              tGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTACTTGGCTCCTGCCCATCCCACGCACCTAACACATACACCAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTAGAGGTAATTCACCAGAGGATGTTGTTCCACTGGcc                                                                                               >  1:299663/1‑301 (MQ=255)
                                                                                                                                                                                                               ggCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGccc                                                                                              <  2:174836/301‑1 (MQ=255)
                                                                                                                                                                                                                 cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAATCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCtt                                                                                            >  1:63973/1‑301 (MQ=255)
                                                                                                                                                                                                                 cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCtt                                                                                            >  2:36564/1‑301 (MQ=255)
                                                                                                                                                                                                                          agtccGTGTATTGATCACCCAGTTTTTGTCTCGGGTTTTTTCACATATTTTTGTTGATTAAGCCACGATACCAGATTGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCa                                                                                   <  1:262052/298‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCACTTGCCACGTCAACTTCGCatta                                                                        >  2:65699/1‑300 (MQ=255)
                                                                                                                                                                                                                                      gATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTCTCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCAAGTCAAATTAGAtttaa                                                                       >  2:364612/1‑299 (MQ=255)
                                                                                                                                                                                                                                       ttCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTcac                                                                      <  2:190387/300‑1 (MQ=255)
                                                                                                                                                                                                                                        tCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAATCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTcaca                                                                     <  2:64041/301‑1 (MQ=255)
                                                                                                                                                                                                                                              gggATGTGTTTCGGTTTTATCCACATATTGTTGTTGAATAGGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCa                                                               <  1:106445/300‑1 (MQ=255)
                                                                                                                                                                                                                                                             tttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGATTTCACAAACTCATTAAacaacgccacc                                                >  1:177505/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                           tttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAAAGTACCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCACCCACGGATACGGGGa                                  >  1:456374/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                 gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCACCCAAGGAAACGGTAAACTaat                            >  1:262330/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                    tAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGCTTCGATAACGTGCAACCTTCATGGTCGAGTTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCCGATTCCATTCCACAAACTCATTAAACGCCGTCCCCCACGGATACGGTGAACCATCAcc                         >  1:176448/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     aaGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACCCATTAAACGACGTCGCCCACGGACACGGTGAACTCTCAccc                        >  1:388780/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                       gCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACCCATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTCCCACCTTCATGGTCGAGGTACTTCACCAGAGGATGTTGTTCCACTGCCCCTTCCCACGTCAGACTCCCCTCCACATACCCATTAAACCACCTCCCCCCCCGATCCCGTTCACTAtacccctt                      >  1:324423/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                                       gCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCGCCCACGGATACGGTGAACTATCACCGTg                      <  2:389294/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         cccgttatcGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCGCCCACGGATACGGTGAACTATCACCGTGcc                    <  2:156329/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         cACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCGCCCACGGATACGGTGAACTATCACCGTGcc                    <  2:300067/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               aCCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCGCCCACGGATACGGTGAACTAACACCCATCCACTTcc              >  1:171992/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                     tGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCTCGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCCCGTCAGATTCGCTTTCACATCCTCATTAAACGACGTCGCCCACGGAAACGGTGAACTAACACCGTTCCACACCACtgtcc        >  1:18153/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                          cTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCGCCCACGGATACGGTGAACTATCACCGTGCCACTTCACTGTGCGCAgg   >  2:458970/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                           tCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGCTTCGCTTTCACATACTCATTAAACGACTTCGCCCACGGACACGGTGAACTATCACCGACCAACTTCACTATGCGAAGGt  >  2:230240/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           
CAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTTGCCACGTCAGATTCGCTTTCACATACTCATTAAACGACGTCGCCCACGGATACGGTGAACTATCACCGTGCCACTTCACTGTGCGCAGGT  >  NC_000913/3181303‑3181902

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