Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,284,646 G→A A86V (GCC→GTC)  actP ← acetate/glycolate:cation symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,284,6460GA98.6% 217.1 / NA 70A86V (GCC→GTC) actPacetate/glycolate:cation symporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (31/37);  minor base C (1/0);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTGGGCGATAAGGTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCT  >  NC_000913/4284405‑4284896
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
caTCTGGGCGATAAGGTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGAAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACatat                                                                                                                                                                                                                                                    >  1:166495/1‑251 (MQ=255)
     gggCGATAAGGTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGt                                                                                                                                                                                                                                                 <  1:327234/249‑1 (MQ=255)
         gATAAGGTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCgg                                                                                                                                                                                                                                           <  1:236778/251‑1 (MQ=255)
          aTAAGGTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATAAAGTCCCCGGc                                                                                                                                                                                                                                          >  1:255806/1‑251 (MQ=255)
             aGGTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAAt                                                                                                                                                                                                                                       <  1:470587/251‑1 (MQ=255)
              ggTAAAGCGCCACCACCACAAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATc                                                                                                                                                                                                                                      >  1:174454/1‑251 (MQ=255)
               gTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCg                                                                                                                                                                                                                                     <  1:300549/251‑1 (MQ=255)
                       ccaccaccaccaGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCg                                                                                                                                                                                                                             >  1:590680/1‑251 (MQ=255)
                        caccaccaccaGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGt                                                                                                                                                                                                                            >  1:455264/1‑251 (MQ=255)
                                    agaACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAg                                                                                                                                                                                                               <  1:530710/251‑1 (MQ=255)
                                       aCCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTg                                                                                                                                                                                                             >  1:381305/1‑250 (MQ=255)
                                         cacaGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGCAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATAAGGCGACGGAGATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTgtt                                                                                                                                                                                                           >  2:333317/1‑248 (MQ=255)
                                                               cGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCACTCAGGAACAAAATGATCGGCCAGCCCACCAGGCAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAACCACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGCTATTGCCGCCTGCGGTGTagtagt                                                                                                                                                                                    >  2:568450/1‑251 (MQ=255)
                                                                    cTTGTTTCAGACGGTAAGTGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATAGGACAACCCACCCGGAAGCACAGCGAGTAACTTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGAACCCCAAGAATGAAGCGGCGGACATATAATAACCGGAAATCGCCAGCCCGTCCTGGAACCACGTGCTCTTGCCGCCCGCggacaagtggccgct                                                                                                                                                                                >  2:5411/1‑237 (MQ=255)
                                                                      tGTTTCAGACGGGAAGAGGCCACATCGGCAAAGGTGTAGCGACCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCACAGCGAGTAAATTAAGCCGGCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTa                                                                                                                                                                               <  1:103341/249‑1 (MQ=255)
                                                                         ttCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACgag                                                                                                                                                                           <  2:493053/250‑1 (MQ=255)
                                                                             gACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACgg                                                                                                                                                                      >  2:232208/1‑250 (MQ=255)
                                                                              aCGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTACCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGGGTAGTAGTCGCTCCGAGAACGt                                                                                                                                                                      >  2:30769/1‑250 (MQ=255)
                                                                              aCGGTAAGAAGCCACATAGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGAAATCAGGAACAAAATGATCGGCCAGCCCACCAGGACGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGGGTAGTAGTCGCTACGAGACCGt                                                                                                                                                                      >  2:600156/1‑250 (MQ=255)
                                                                               cGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGACCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTAc                                                                                                                                                                    <  1:134176/251‑1 (MQ=255)
                                                                                ggTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTGACGCAGCCGTTCGGCAATCAGGAACAAAATGAGCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACGCCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGAGATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTAcg                                                                                                                                                                   <  2:361023/251‑1 (MQ=255)
                                                                                 gTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTAcgc                                                                                                                                                                  <  2:597065/251‑1 (MQ=255)
                                                                                       ggCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTg                                                                                                                                                            <  1:428935/251‑1 (MQ=255)
                                                                                         ccACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGACAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTCCGCGTTTTGAt                                                                                                                                                          >  1:26250/1‑251 (MQ=255)
                                                                                          cacaTCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATg                                                                                                                                                         <  2:107592/251‑1 (MQ=255)
                                                                                           acaTCGGAAAAGGAGTAGCGCACCAAGTTAAGCAGACGTTAGGCAATCAGGAACAAAATGATCGGCAAGCCAACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCACCAAGAATGAGGCGACGGACATATAGACCCCGGCAATCGCCAGCCAGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTCCGCGtttggaggc                                                                                                                                                        >  2:84689/1‑245 (MQ=255)
                                                                                                    aaGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAAGGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGt                                                                                                                                               <  2:542315/251‑1 (MQ=255)
                                                                                                            gcgcCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGa                                                                                                                                       <  1:196329/251‑1 (MQ=255)
                                                                                                             cgcCCCAGGTTACGCAGACGTCCGGCAATCAGAAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATGAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAg                                                                                                                                      <  2:585517/251‑1 (MQ=255)
                                                                                                              gcCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGc                                                                                                                                     <  1:284271/251‑1 (MQ=255)
                                                                                                              gcCCCAGGCTACGCAGGCGTTCGGCAACCAGGACCAAAAAGATCCGCCAGCCCACCAAGAAGCCCGGCGCGAACATTCAGACGTCATAGCCGGCGCCAACCACCCGCGCGGAGATCCCCCCGAATGAGGCGCCCGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATCTTGCCGCCTGCGGTGTAGTAGTCGCTCCGCGAACGTACGCGTTTTGAAGCCCAGTAGGTACTGCCGAGc                                                                                                                                     >  2:322932/1‑251 (MQ=255)
                                                                                                                       tACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACacg                                                                                                                            >  1:519090/1‑251 (MQ=255)
                                                                                                                       tACGCAGACGTTAGGCAATCAGGAACAAAATGCTCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGCCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCTCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTaacaacg                                                                                                                            >  1:188188/1‑246 (MQ=255)
                                                                                                                         cGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACacgac                                                                                                                          <  1:414979/251‑1 (MQ=255)
                                                                                                                          gCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACacgacg                                                                                                                         >  1:405427/1‑251 (MQ=255)
                                                                                                                                    ggCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGa                                                                                                               <  2:183482/251‑1 (MQ=255)
                                                                                                                                     gCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGaa                                                                                                              <  2:381305/251‑1 (MQ=255)
                                                                                                                                     gCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGaa                                                                                                              <  1:248659/251‑1 (MQ=255)
                                                                                                                                      aaaTCAGGAACAAAATGAGCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACAGCAGCGCGGATATCCCCAAGAAGGAGGCGAGGGATATATAGTCCCGGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGACATTGCCGCCTGCGGTGTAGTAGTCGGTTCGAGAACGTACGCTTTTAGATGACCAGTAGGCAAGGCCGTGCGTAAACACGACGaaatgaaggaac                                                                                                             <  2:398439/250‑10 (MQ=255)
                                                                                                                                          cAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCTGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACata                                                                                                          <  1:428018/250‑1 (MQ=255)
                                                                                                                                          cAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATGTAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAACGCCGAGCGTAAACACTACGAAAATCAGGAACata                                                                                                          <  1:431668/250‑1 (MQ=255)
                                                                                                                                            ggAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACataat                                                                                                        >  1:481691/1‑250 (MQ=255)
                                                                                                                                                     tGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGCATGAGGCGACGGACATATAGTCACCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTCCGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGc                                                                                              >  1:537892/1‑251 (MQ=255)
                                                                                                                                                       aTCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCa                                                                                            <  1:232208/251‑1 (MQ=255)
                                                                                                                                                               gCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCGGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCTAGCGGAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTtg                                                                                    <  1:613521/251‑1 (MQ=255)
                                                                                                                                                                cccACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTtgg                                                                                   >  2:36566/1‑251 (MQ=255)
                                                                                                                                                                cccACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTGCGTtgg                                                                                   <  2:26250/251‑1 (MQ=255)
                                                                                                                                                                   accaGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCt                                                                                                >  1:70198/1‑235 (MQ=255)
                                                                                                                                                                   accaGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGGTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCt                                                                                                <  2:70198/235‑1 (MQ=255)
                                                                                                                                                                      cggAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGgcg                                                                             <  1:551354/250‑1 (MQ=255)
                                                                                                                                                                          aGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGGGTAGTAGTCGCTACGAGAACGTACGCGTGTTGATGCCCAGTAGTTAATGCCGGGCGTAAACACGACGAAAATCAGGAACATAATGAGCGCCTGCCAGTTCGTTGGCGGGCGctcg                                                                         <  2:61907/251‑2 (MQ=255)
                                                                                                                                                                            cccAGCGAGTAAATTACGCCGTCACAGCCGGAGGTAAACACCAGCGCGGAGCTCCCCAAGAATGAGGCGACGGACAGACAGTCCCCGGCAATCGCCAGCCCGGTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGGAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTGCGTTGGCTGGCGCTCTAc                                                                       <  2:430361/251‑1 (MQ=255)
                                                                                                                                                                                          gaaGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATa                                                          <  1:240171/249‑1 (MQ=255)
                                                                                                                                                                                              ccGTCATAGCCGGAGGTAAACCCCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATc                                                     <  2:455264/251‑1 (MQ=255)
                                                                                                                                                                                                 tCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCcgc                                                  >  1:480726/1‑251 (MQ=255)
                                                                                                                                                                                                   aTAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCgcgg                                                >  2:493319/1‑251 (MQ=255)
                                                                                                                                                                                                      gCCGGCGGTAAAAACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTCGTCGCTACGAGAACGTACGAGTTTTGATGCCCAGTAGGTAATGCCGAACGTAAACACGACGAAAATCAGGACCATAATAATCGCCTGCCAGTTCGTTGGATGGCGATCTACGGCCCCGCTAATAGCATCCGCGGCGt                                             >  2:624905/1‑251 (MQ=255)
                                                                                                                                                                                                        cGGAGGGAAACACCAGCGCGGAGAGCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTa                                           <  1:514014/251‑1 (MQ=255)
                                                                                                                                                                                                         ggAGGTAAACACCAGCGCGGAAATCCCCAAGAATGAGGCGACGGACATATAGTACCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATCATCGCCTGCCAGTTCGTTGGCTGGAGCTCTACGGCCCCGCTAATAGCATCCGCGGCGtt                                            >  2:578788/1‑249 (MQ=255)
                                                                                                                                                                                                           aGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCt                                        >  2:176492/1‑251 (MQ=255)
                                                                                                                                                                                                           aGGTAAACACCAGCGCGGAGATCCCCAAGAAGGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCt                                        <  2:94049/251‑1 (MQ=255)
                                                                                                                                                                                                            ggTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTg                                       >  2:293184/1‑251 (MQ=255)
                                                                                                                                                                                                                         gcgGAGATCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGGAgtgg                          >  1:332306/1‑250 (MQ=255)
                                                                                                                                                                                                                             agaTCCCCAAGAATGAGGCTACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGTAGTGTggcg                      <  2:520574/251‑1 (MQ=255)
                                                                                                                                                                                                                                tCCCCAAGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACAAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAACGGGAGTGTggcggc                    >  2:258102/1‑250 (MQ=255)
                                                                                                                                                                                                                                      aGAATGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGcc             >  1:191283/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaTGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACataa                                                                                                         >  1:174749/1‑157 (MQ=255)
                                                                                                                                                                                                                                        aaTGAGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACataa                                                                                                         <  2:174749/157‑1 (MQ=255)
                                                                                                                                                                                                                                            aGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGACGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTCCGCGTTTTGACGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTCAGa       <  1:333317/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aGGCGACGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTc          >  2:610988/1‑248 (MQ=255)
                                                                                                                                                                                                                                                 gcggACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATGGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGGAGTGTGGCGGCAATCGCCGTCAGAActct  <  1:492550/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
CATCTGGGCGATAAGGTAAAGCGCCACCACCACCAGAGAACCACAGGCCGAAAGAATACGAATCGGCCCTTGTTTCAGACGGTAAGAGGCCACATCGGCAAAGGTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCTGCCAGTTCGTTGGCTGGCGCTCTACGGCCCCGCTAATAGCATCCGCGGCGTTAGCTGCGAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCT  >  NC_000913/4284405‑4284896

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