Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 784,462 C→T G121D (GGC→GAC)  zitB ← Zn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913784,4620CT98.5% 204.2 / NA 66G121D (GGC→GAC) zitBZn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/32);  minor base G (0/1);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCGCTACGCAGAACCAGAAGCGACACCAGTATCGAGAGAATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCAGAAAGAAA  >  NC_000913/784219‑784706
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cACGCGCTACGCAGAACCAGAAGCGACACCAGTATCGAGCGAATGGGGTCAGAAGGTGTCCAGCAAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCAAGCAGGTCCCCGAGCACATGCAGTGCAGCGGCTCGACCGTTGAGGTCTTTCTCTTCGCTGCCGTGATGAAGTACCCAGAAAGAATGTATATTTGCCAGCAACCCAGCCACTGCAATTACCATAATCATACCGTCCTCgcc                                                                                                                                                                                                                                               >  2:328838/1‑249 (MQ=255)
         cGCAGAACCAGAAGCGACACCAGTATCGAGAGAATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACgcgg                                                                                                                                                                                                                                      >  2:370611/1‑251 (MQ=255)
             gAACCAGAAGCGACACCAGTATCGAGAGAATGGGGTCCGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGAACCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTACCCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACAGGACGCGGCGTg                                                                                                                                                                                                                                  >  2:276830/1‑251 (MQ=255)
                  aGAAGCGACACCAGTATCGAGAGAATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGa                                                                                                                                                                                                                              >  1:321488/1‑250 (MQ=255)
                     aGCGACACCAGTATCGAGAGAATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCg                                                                                                                                                                                                                          <  2:232159/251‑1 (MQ=255)
                             cAGTAGCGAGAGAATGGGGTCAGCAGGCGTCCAGCCAGTCCAGATAATAATCAGCGCGGCAATAATGGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCACGTACGTTTAGGTTTTTCCCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCACCCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGGGCGGAAGCGTTCTATCg                                                                                                                                                                                                                  <  2:156215/251‑1 (MQ=255)
                                 aTCGAGAGAATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGCTGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGACATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGACGCGTTCTATCGCCTc                                                                                                                                                                                                              >  2:63691/1‑251 (MQ=255)
                                    gagagaATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGAATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTccc                                                                                                                                                                                                            >  2:112967/1‑250 (MQ=255)
                                     ggagaATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAg                                                                                                                                                                                                          <  2:335734/250‑1 (MQ=255)
                                      gagaaagggggaagcagGTGTCCAGCCAGTCCAGATAATAACCAACGCGGCCATAATAGCACCAACCGCACCCAGCAGGTCCCCGAGCACATGAAGGGCCGCTGCCCGTACGTTGAGGTTTTTCTCTTAGCTGCCGAGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCACCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGCCGCGGCGTGCGGAAGCGTTCTACCGCCTccc                                                                                                                                                                                                            >  2:105664/13‑248 (MQ=255)
                                        gaATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACa                                                                                                                                                                                                       <  1:477050/251‑1 (MQ=255)
                                             ggtgccgcacgcggccacccggTCCACATAATAATGCACGCGGCATAAATACCCCCATCCGAGCCCAGCAGTCCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTGCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCGCGGCAATTGCCATCATCATGCCGTCCTCGACCTGACGCTGCGTGCGGACGCGTTCTATCGCCGCCCAGACAATCaa                                                                                                                                                                                                  <  2:469988/230‑1 (MQ=255)
                                                    cagGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTa                                                                                                                                                                                           <  2:182309/251‑1 (MQ=255)
                                                     agGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTa                                                                                                                                                                                           >  1:555389/1‑250 (MQ=255)
                                                     agGTGTCCAGCCAGGCCCGATAATAATCCACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTaa                                                                                                                                                                                          <  2:484471/251‑1 (MQ=255)
                                                       gtgtCCAGCCAGTTCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTTCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATTGTAAtt                                                                                                                                                                                        <  2:102473/251‑1 (MQ=255)
                                                               ccagTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAg                                                                                                                                                                                <  1:277319/251‑1 (MQ=255)
                                                                cagTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCTAACCCAGCAGGTCCCTGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCATCGCTGTCGTGATGAAGTAAACAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCGATCGCCTCCCAGACAATCAAAAGGGTAATTACCCCGAgg                                                                                                                                                                               <  1:524378/251‑1 (MQ=255)
                                                                   tCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCaaa                                                                                                                                                                                                 <  1:448287/230‑1 (MQ=255)
                                                                   tCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCaaa                                                                                                                                                                                                 >  2:448287/1‑230 (MQ=255)
                                                                              atCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAATCGTTCTATCGCCTCCCAGACAATCAAAAt                                                                                                                                                                                               <  1:171138/221‑1 (MQ=255)
                                                                              atCAACGCGGCAATAATAGCCCCAACCGAACCCAGAAGGTCCCCGAGCACATGCAGTGCCGCTGATCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAAt                                                                                                                                                                                               >  2:171138/1‑221 (MQ=255)
                                                                                cAANGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGGTGCTCGTACGTTGAGGTTGTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTTCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCTCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACaa                                                                                                                                                               <  2:586100/251‑1 (MQ=255)
                                                                                  aCGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAAc                                                                                                                                                             >  1:159670/1‑251 (MQ=255)
                                                                                  aCGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAAc                                                                                                                                                             <  1:302681/251‑1 (MQ=255)
                                                                                         aataatAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCcgcg                                                                                                                                                       >  2:32580/1‑250 (MQ=255)
                                                                                                             ccAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCa                                                                                                                                   >  2:160397/1‑250 (MQ=255)
                                                                                                                         cGAGCACATGCAGGGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGCCAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTgg                                                                                                                      >  2:142512/1‑251 (MQ=255)
                                                                                                                                      tGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTggg                                                                                                           >  2:385906/1‑249 (MQ=255)
                                                                                                                                       gCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGggagg                                                                                                        <  2:247370/251‑1 (MQ=255)
                                                                                                                                             gctCGTTCGTTGAGGTTTTTCTCTTCGCGTCCTTGATGAAGTAACCAGAAAGAAATTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGgacgac                                                                                                   <  2:34269/250‑1 (MQ=255)
                                                                                                                                             gctCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGacgacg                                                                                                  >  1:647349/1‑251 (MQ=255)
                                                                                                                                                  tACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAAtg                                                                                                                >  1:292366/1‑232 (MQ=255)
                                                                                                                                                  tACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAAtg                                                                                                                <  2:292366/232‑1 (MQ=255)
                                                                                                                                                   aCGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGaaa                                                                                             <  2:647349/250‑1 (MQ=255)
                                                                                                                                                   aCGTTGAGGGTTTTCTCTTCGCTGCAGTGATTAAGTAACCAGAAAGAAAGTATATTTGCCAGCAGCCCAGCCACGGAAATATCCAGCTGCATGCCGTCCTCGACCGGCCGCGTCGTGACGAAGTCTTCTATCGCCTCCCAGACAAGCAAAATGGTATTGACAACGAGGGCGATGTCATTCACAACCGCCGCGAGAGTGGTCAGTCTTAGACAACCGAAAGTGTGTCGAAGGGTAggacgaaggctttata                                                                                             <  2:23680/250‑14 (MQ=255)
                                                                                                                                                             ttttCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGAGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCg                                                                                  >  2:47080/1‑251 (MQ=255)
                                                                                                                                                                   ctTCGCTGCCGTGAAGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGAGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACAACGAGGGCTCTGGCATTCACAAACGCAGCGAGCGTGGTCAGTCTGAGCCAACCGAACGTGTGGCGAATGGTGGGAGGACGACGCGAAAATTGCACGGCGAANAgg                                                                            >  2:602891/1‑251 (MQ=255)
                                                                                                                                                                     tCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTAGATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAggg                                                                           <  2:648130/250‑1 (MQ=255)
                                                                                                                                                                      cGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGc                                                                          <  1:509740/250‑1 (MQ=255)
                                                                                                                                                                       gCTGCCGTGATGAAGTAACCAGAAAGAAAGGATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTACTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGAGGGAGGACGACGTGAAAATTGCACGGCGAGCAGGGCg                                                                         >  2:619391/1‑250 (MQ=255)
                                                                                                                                                                                  gaAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGggc                                                              <  2:243547/250‑1 (MQ=255)
                                                                                                                                                                                     gTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGggcggc                                                           >  1:634893/1‑250 (MQ=255)
                                                                                                                                                                                       aaCCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGt                                                        >  2:516021/1‑251 (MQ=255)
                                                                                                                                                                                        aCCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGt                                                        >  2:560858/1‑250 (MQ=255)
                                                                                                                                                                                                gaaagTATATTTGCCAGCAACCCATCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCaa                                               <  2:264357/251‑1 (MQ=255)
                                                                                                                                                                                                  aagTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGGGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAAc                                              >  2:418444/1‑250 (MQ=255)
                                                                                                                                                                                                   tgTATATTTGCCAGCATCCCAGCCACGGCAATTGCCATCATCATGCCGTCCTAGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGCAATTACCACGAGGGCTATGGCAGTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAAAGGTGGGAGGACGACGTGAAACTTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACat                                            >  2:102678/2‑251 (MQ=255)
                                                                                                                                                                                                   agTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACat                                            >  2:623195/1‑251 (MQ=255)
                                                                                                                                                                                                         tttGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTAACACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGACTGGGGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATgccc                                      >  2:35599/1‑248 (MQ=255)
                                                                                                                                                                                                           tGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTc                                                                                                                                                                                                                       >  2:98464/1‑72 (MQ=255)
                                                                                                                                                                                                           tGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTc                                                                                                                                                                                                                       <  1:98464/72‑1 (MQ=255)
                                                                                                                                                                                                            gCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTAACGGTCAACATATGACCCGc                                   >  2:217748/1‑251 (MQ=255)
                                                                                                                                                                                                                 cAACCCAGCCCCGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCgg                              <  1:145918/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaCCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGc                             <  2:28787/251‑1 (MQ=255)
                                                                                                                                                                                                                   aCCCAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGcc                            >  1:152724/1‑251 (MQ=255)
                                                                                                                                                                                                                    cccAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGGGCGTGGTCAGTCTTAGCCAGCCGAAAGTGTGGCGAATGGGGGGATGACGACGTGAAAATTGCACGTCGAGTGGGGCGAAAAGCAGGGCGGCAGTATCGGCGAACTTATTACCCGCATCGGCCa                           <  2:261028/251‑1 (MQ=255)
                                                                                                                                                                                                                      cAGCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGc                         >  1:482100/1‑251 (MQ=255)
                                                                                                                                                                                                                        gCCACGGCAATTGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCaa                       <  2:321488/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccACGGCAATTGCCATCATCATGCCGGCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCTCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAAt                      <  1:47080/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggCAATTGCCATCATCATGCCGTCCTCGACCGGACGAGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAAGCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCa                  <  2:131320/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ttGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTAACACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGAAGTGAAAATTGCAAGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGATAACATATGACCCGAATCGGCCAGCAATGCCAgaga              >  2:631507/1‑250 (MQ=255)
                                                                                                                                                                                                                                   tGCCATCATCATGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCTGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGaa             <  1:619391/250‑1 (MQ=255)
                                                                                                                                                                                                                                         catcatGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCagaa       <  1:516021/250‑1 (MQ=255)
                                                                                                                                                                                                                                           tcatGCCGTCCTAGCCTGGACGCGGAGAGAGGAAGAGTTATATCGCCTCCCAGACAATCAACATGGAAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTAAGTCTTAGCCAACCGAAAGCGTGGAGAATGGTGGGGGGACGACGAGAAAATTGCACGGGGCGTAGGGCGAAATGCAGGGCGGAAGAATCGGTCAACATATGAACCGCATCGGCCAGAAAtagccgcgaaacagaaga     >  2:362292/1‑233 (MQ=255)
                                                                                                                                                                                                                                             atGCCGTCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCAgaaagaaa  <  1:259049/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CACGCGCTACGCAGAACCAGAAGCGACACCAGTATCGAGAGAATGGGGTCAGCAGGTGTCCAGCCAGTCCAGATAATAATCAACGCGGCAATAATAGCCCCAACCGAACCCAGCAGGTCCCCGAGCACATGCAGTGCCGCTGCTCGTACGTTGAGGTTTTTCTCTTCGCTGCCGTGATGAAGTAACCAGAAAGAAAGTATATTTGCCAGCAACCCAGCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCAGAAAGAAA  >  NC_000913/784219‑784706

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