Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,608,416 G→A A656T (GCC→ACC)  zntA → Zn(2(+))/Cd(2(+))/Pb(2(+)) exporting P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,4160GA91.2% 191.6 / ‑5.0 68A656T (GCC→ACC) zntAZn(2(+))/Cd(2(+))/Pb(2(+)) exporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/37);  minor base C (0/6);  total (25/43)
Fisher's exact test for biased strand distribution p-value = 7.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGTCAAAGGGGTGATCCTCACCGGCGATAATCCACGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG  >  NC_000913/3608175‑3608658
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
gggCGTCAAAGGGGTGATCCTCACCGGCGATAATCCCCGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACAGGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGAGGGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACgcag                                                                                                                                                                                                                                           <  2:165812/251‑1 (MQ=255)
    gTCAAAGGGGTTATCCTCACCGGCGATAATCCACGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACTCGCCAGCGATGAAATCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCAtt                                                                                                                                                                                                                                       <  1:215158/251‑1 (MQ=255)
      cagtggggTGTTCCTCACCGTCGATAATCCTCGCGCAGCGGCGGCAATTGCCGGGGGGCTGGGGCTGGAGTTTAAAGCGGGACTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTTAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCAGCGGTATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTaa                                                                                                                                                                                                                                     <  2:94073/247‑1 (MQ=255)
            ggTGATCCTCACCGGCGATAATCCACGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATa                                                                                                                                                                                                                               <  2:147090/251‑1 (MQ=255)
             gTGATCCTCACCGGCGATAATCCACGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATaa                                                                                                                                                                                                                              <  1:47807/251‑1 (MQ=255)
                     cACCGGCTATAATCCACGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTgc                                                                                                                                                                                                                      <  2:235825/251‑1 (MQ=255)
                       ccGGCGATAATCCACGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTgcgc                                                                                                                                                                                                                    <  1:106078/251‑1 (MQ=255)
                               aaTCCACGCGCAGCGGCGGCAATTGCCGGGGCGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGGTTGCAATGGGTAGCGGCCCAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTCCGCGGCCTggg                                                                                                                                                                                                            >  2:170717/1‑250 (MQ=255)
                                     ggcgcagcggtggtaggtggcGGGGGGCTGGGGCGGGAGTTTAAAGCGGGCGTGTTGCCGGAAGATAAAGGCAAAGCTGTGACCGAGCTGAATCAACATGCGCCGCGGGCGATGGTCTGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATACCCACCTGCGCGGCCTTGTGCAAAt                                                                                                                                                                                                      <  2:97507/231‑1 (MQ=255)
                                         cagggggggcAATTGCCGGGGAGCTGGGGCGGGAGTTTAGAGCGGGCCTGTGGGCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGAGTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGAtt                                                                                                                                                                                                  <  2:197792/244‑1 (MQ=255)
                                         aaGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGAATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGGGGCGCTGGAAACCACCGACGCAGCATGAACCCATAACCACCTGCGCGGCCTGGTGCAAATGAtt                                                                                                                                                                                                  >  2:142290/2‑251 (MQ=255)
                                                   aaTTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCAc                                                                                                                                                                                        <  1:90685/251‑1 (MQ=255)
                                                           gggAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACt                                                                                                                                                                                >  1:128470/1‑251 (MQ=255)
                                                                 gggggCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGGAGCATTATCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGcc                                                                                                                                                                          <  1:112648/250‑1 (MQ=255)
                                                                    ggCTGGAGTTTAAAGCGGGCCTGTTGCGGGAAGAGAAAGGCAAAGCGGGGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATAAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCCCCGGGGCGTCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAat                                                                                                                                                                       <  1:175144/251‑1 (MQ=255)
                                                                         gAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAACTGATTGAACTGGCACGCGCCACTCACGCCAATATccc                                                                                                                                                                  >  2:124511/1‑250 (MQ=255)
                                                                              taacgcgcgCCTGTTGCCCGAAGATAACGTCCAAGCGCTGACCGTGATGCATCAACATGCGCCGCTGGCGATGATAGGAGACGGTATTAACGAAGCGCCAGCGATGAAAGCTGCCTCCATCGGGATTGCAATGAGAAGCGGCACAGACGTGGCGCTGGCAACCACCGACGCAGCATTAACCCATAACCACCTGACAGGCCTGGTGCAACAGATTGAAACGGCACGCGCCACTCCCGCCAATATCCCCCAgc                                                                                                                                                             >  2:111421/9‑250 (MQ=255)
                                                                                   cGGGCCTGTTGCCGTAAGATAAAGTCAAAGCGGTGACCGAGCTAAATCAACATGCGCCGCTGTCGATGGTGGGTGACGGTATTAACGACGCGCCAGCTATGAAAGCTGCCGCCAGCGGGATGGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGGCGCTGCATTAACCCATAACCACCTGCGCGGCCTGGTTCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATc                                                                                                                                                        <  1:278728/251‑1 (MQ=255)
                                                                                   cGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGGGGCGCTGGAAACCACCGACGCAGCATTACCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATc                                                                                                                                                        <  2:40411/251‑1 (MQ=255)
                                                                                     ggCCTGTTGCCGGCAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCAACTGCGCGGCCTGGTGCAAATGATTGCACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCAc                                                                                                                                                      >  2:227934/1‑251 (MQ=255)
                                                                                          gTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTg                                                                                                                                                 <  1:111323/251‑1 (MQ=255)
                                                                                            tgcaggAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGCCGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTgcg                                                                                                                                               >  1:137366/5‑251 (MQ=255)
                                                                                            ggCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTgcg                                                                                                                                               <  1:82187/250‑1 (MQ=255)
                                                                                            gccctgaaaataaagcaaaatcGGTGACCGAGCGGATTCAACATGGGCCTCTGGGTATGGGCGGTGACGTTATTAACGACGCGCCAGCGATAAAAGCTGCCGCAACGGGGATTTCAATGGGTTGCGTCACAGCCGTGGCGCTGGAAACCACCGCCTAAGCATTAACCCAAAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACGCACGCCAATATCCGCCAGAACATCACTATTgcg                                                                                                                                               <  2:175926/230‑1 (MQ=255)
                                                                                                   aGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACCTGCGCCGCTGGCGATGGTCGGTGACGGACTTAAAGACGCGCCAGCGATGCAACCTCCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCCTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGc                                                                                                                                        >  2:134291/1‑251 (MQ=255)
                                                                                                       aaaGTCAAAGCGGTGACCGAGCTGAGTCAACAGGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCCCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGa                                                                                                                                     <  1:88117/250‑1 (MQ=255)
                                                                                                          ggcAAAGCGGTGAGCGAGCAGACTCACCCGGCGCCGCTGGGGAGGGTGGTGAGCGGTAGAAACGACGCGCAAGGAAGGAAAGTTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCAAGGCAGCATTACCCCATACCCACCGGCGCGGCCGGGTGCAAATGATGGACGTGGCTCGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAgg                                                                                                                                 <  2:3355/249‑1 (MQ=255)
                                                                                                            cAAAGCGGTGACCGAGCGGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGa                                                                                                                               <  1:208066/251‑1 (MQ=255)
                                                                                                               aGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCt                                                                                                                            >  2:77966/1‑251 (MQ=255)
                                                                                                                  ggTGCCCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTcc                                                                                                                         <  2:128470/251‑1 (MQ=255)
                                                                                                                  ggTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTcc                                                                                                                         >  1:197894/1‑251 (MQ=255)
                                                                                                                      cccGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGt                                                                                                                     <  1:253330/250‑1 (MQ=255)
                                                                                                                         gagctgcctCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTACAGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGCGCAAATGATTGAACTGGCACGCGCCACTAACGACAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTcac                                                                                                                  >  2:251412/9‑251 (MQ=255)
                                                                                                                         gAGCTGAATCACCATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGACAGCTGCCGCCGTCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGCCACGCGCCACTCACGCCAATATCCGCCAGAACCTCACTATTGCGCTGGGGCTGTAAGGGATCATCCTcatcgc                                                                                                                  >  2:53715/1‑246 (MQ=255)
                                                                                                                            cTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTc                                                                                                                       <  1:243828/243‑1 (MQ=255)
                                                                                                                            cTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTc                                                                                                                       >  2:243828/1‑243 (MQ=255)
                                                                                                                             tGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACg                                                                                                              >  1:107034/1‑251 (MQ=255)
                                                                                                                                   aaCATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGtt                                                                                                         >  1:150848/1‑250 (MQ=255)
                                                                                                                                   aaCAGGCGCCGTGGGCGAGGGTCGGTGACGGTATTAACGACGCGCCGGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCATCGACGCAGCATTAACCCATAACCCCCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTa                                                                                                        <  2:60852/251‑1 (MQ=255)
                                                                                                                                                      ggTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGtggc                                                                                     >  1:244860/1‑251 (MQ=255)
                                                                                                                                                      ggTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGtgg                                                                                      <  1:7436/250‑1 (MQ=255)
                                                                                                                                                             gACGGTATTAACGACGCGCCAGCGATGAAAGCGGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCCCCGACGCAGCATTAACCCATAACCCCCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGTGATGACCGGGTTGTGGCTGGCAGt                                                                              <  2:18246/251‑1 (MQ=255)
                                                                                                                                                                  gTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTg                                                                          <  1:179488/251‑1 (MQ=255)
                                                                                                                                                                      tAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGc                                                                        >  2:225018/1‑249 (MQ=255)
                                                                                                                                                                          gacgCGCCAGCGATGCGAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGGAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGacac                                                                  >  2:16192/1‑248 (MQ=255)
                                                                                                                                                                          gacgCGCCAGCGATGAAAGCGGCCGCCAGCGGGATGGCAAGGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGGTGACCGGGTTGTGGCTGGCAGTGCTGGCAGATAc                                                                  <  1:150032/251‑1 (MQ=255)
                                                                                                                                                                                ccAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCCCCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGc                                                            <  1:77966/251‑1 (MQ=255)
                                                                                                                                                                                        gAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGGGGCGCTGGAACCCACCGACGCAGCATTAACCCATAACCCCCTGCGCGCCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTg                                                     <  1:218306/250‑1 (MQ=255)
                                                                                                                                                                                        aaaaccTGCCGCCATCGGGATTGCAAGGGGTAGCGGCACAACCGTGGCGCTGGAATCCCCCGCCGCAGCATTACCCCATACCCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCTCGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCGCGTCCCCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGc                                                    <  2:137366/246‑1 (MQ=255)
                                                                                                                                                                                         aaaGCTGCCGCCATCGGGATTGCGATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCt                                                   >  1:210979/1‑251 (MQ=255)
                                                                                                                                                                                           aGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGGGGCGCTGGAAACCACCGACGCAGCATTACCCCATAACCACCTGCGCGGCCGGGTGCAAATGATTGACCTGGCACGCGCCACGCACGCCAATATCCGCCAGAACATCTCTAGTGCGCTGGGGCTGAAAGGGATCTGCCGCGTCACCACGCTGGTAGGGATGAGCGGGGTGTGGCTGGCAGTGCTGGCAGATACGTGGGCGACGGTGCTgg                                                 <  2:113264/251‑1 (MQ=255)
                                                                                                                                                                                            gCTGCCGCCACCGGGATTGCAAGGGGTACCGGCACCGACGTGGCGCTGGAACCCACCGTCGCTGCATTAACCCCTACCCCCCTGCGCGGCCGGGTGCAAATGATTGAACTGGCTCGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCATGCGGGTAGGGATGACCGGGTTGTGGCTGGCAGTCCTGGCCGATCCGGGGGAGACGGTGCTGGt                                                <  2:10406/251‑1 (MQ=255)
                                                                                                                                                                                                   ccATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAACCCACCGACGCAGCATTAACCCATAACCCCCTGCGCGGCCTGGTGCAAATGATTGAACTGGCCCGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCg                                         <  1:198154/251‑1 (MQ=255)
                                                                                                                                                                                                     aTCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGa                                        >  1:136186/1‑250 (MQ=255)
                                                                                                                                                                                                      tCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCCCCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGaa                                                                                                                                    <  2:216678/157‑1 (MQ=255)
                                                                                                                                                                                                      tCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGaa                                                                                                                                    >  1:216678/1‑157 (MQ=255)
                                                                                                                                                                                                       cGGGATTGCAATGTGTAGCGGCACAGTGGGGGCTCTGGAAACCCCCGCCGCAGCATTACCCCATAACCACCTGCGCGGCCGGGTGCAAATGATGGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTGCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAAt                                      <  2:181586/250‑1 (MQ=255)
                                                                                                                                                                                                           aTTGCAATGGGGAGCGGAACAGACGTGGCGCTGGACACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGCACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCCGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGtt                                 >  2:269246/1‑251 (MQ=255)
                                                                                                                                                                                                                    ggTAGCGGCACAGAGGGGGCGCGGGAAACCACCGACTCAGCATTAACCCATACCCACCTGCGCGGCCTGGTGCAAATGAGTTATCGGGCACGCGCCACTCACGCCAATAGCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGGTAGGGAGGACCGGGTTGTGGCTGTCAGTGCGGGCATATATGGGGGCGACGGTGCTGGTGACAGCCAAGGCGTTAAGAttgtt                        <  2:248788/251‑1 (MQ=255)
                                                                                                                                                                                                                      tagcgccaccgccccggCGCTGGCAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGACCTGCCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGACTGAAAGGGATCTTCCTCGCCCCCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTCGCAGATACGGGGGCGACGGGGCTGGTGACAGCGAACGGGTTAAGATTGTTgc                      >  2:252686/16‑251 (MQ=255)
                                                                                                                                                                                                                       aGCGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTgcg                     <  1:137874/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGGCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCCGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTgcgc                    >  1:150608/1‑250 (MQ=255)
                                                                                                                                                                                                                          ggCACAGCCGTGGCGCTGGAAACCCCCGCCGCAGCATTACCCAATACCCCCCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAg                  <  2:145489/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggCACAGACGTGGCGCTGGAAACCACCGACGCAGCATTAACCCATACCCACCGGCGCGGCCTGGGGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAg                  <  2:164587/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gTGGCGCTGGAAACCACCGACGTATCATTACCCCATAACCCCCTGCGCGGCCTGGTGCAAATGATTGGACTGGCCCGCGCCACTCCCGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCTGGGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAgt         <  2:160047/251‑2 (MQ=255)
                                                                                                                                                                                                                                   gTGGCGCTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCCGATACGGGGGCGACGGGGCTGGTGAACGCGAATGCGGTTAGATTGTTGCGCCGgagagcagg         >  2:4887/1‑246 (MQ=255)
                                                                                                                                                                                                                                      gcgcTGGAAACCACCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCaa      <  1:170559/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gggAAACCCCCGACGCAGCATTAACCCATAACCCCCGGCGCGGCCGGGTGAAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCg  <  2:170629/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
GGGCGTCAAAGGGGTGATCCTCACCGGCGATAATCCACGCGCAGCGGCGGCAATTGCCGGGGAGCTGGGGCTGGAGTTTAAAGCGGGCCTGTTGCCGGAAGATAAAGTCAAAGCGGTGACCGAGCTGAATCAACATGCGCCGCTGGCGATGGTCGGTGACGGTATTAACGACGCGCCAGCGATGAAAGCTGCCGCCATCGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG  >  NC_000913/3608175‑3608658

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