Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,894,906 T→A 22.6% T236S (ACC→TCC)  yiaN ← 2,3‑diketo‑L‑gulonate transporter large permease YiaN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,894,9060TA22.6% 140.6 / 10.6 65T236S (ACC→TCC) yiaN2,3‑diketo‑L‑gulonate transporter large permease YiaN
Reads supporting (aligned to +/- strand):  ref base T (44/4);  new base A (3/12);  total (48/17)
Fisher's exact test for biased strand distribution p-value = 2.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACGCGGTGAATCGACCAACGGTTGCAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGTTGCGCCCATCATCATGCCAGGGGCAATGCCCGCCATAAACAGCTTACTGATGGATAATCCGCTGGAAACGCCGA  >  NC_007779/3894670‑3895136
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
cagACGCGGTGAATCGACCAACGGTTGCAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTaa                                                                                                                                                                                                                                      >  1:1498893/1‑239 (MQ=255)
     gcgGTGAATCGACCAACGGTTGCAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGACATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGAAAACAg                                                                                                                                                                                                                                  >  1:2152957/1‑238 (MQ=255)
       ggTGAATCGGCCAACGGTTGCAGTAAATCCGAAACCCTCATCGGCAATTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCCCTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGcc                                                                                                                                                                                                                               <  2:890862/239‑1 (MQ=255)
       ggTGAATCGACCAACGGTTGCAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGc                                                                                                                                                                                                                                <  1:675778/238‑1 (MQ=255)
                     cGGTTGCAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAAAAACTGGGGCGACAAACAGCGCATAGAAAGAAGCAACGGCCCCGGCTTCGGAGGGGGAAAAAAGGCCGGGGCGAACCCCg                                                                                                                                                                                                                  >  2:217227/1‑238 (MQ=255)
                       gTTGCAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCa                                                                                                                                                                                                               >  1:1680126/1‑239 (MQ=255)
                           cAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCaataa                                                                                                                                                                                                           >  1:522319/1‑239 (MQ=255)
                           cAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCaata                                                                                                                                                                                                            >  1:2040490/1‑238 (MQ=255)
                           cAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGGGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGGAAACAGGCCGGAGCGAAACCCGCCaataa                                                                                                                                                                                                           >  2:186217/1‑239 (MQ=255)
                                aaTCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAatgatg                                                                                                                                                                                                      >  1:1260171/1‑239 (MQ=255)
                                      aaaCCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGAGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGCAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGga                                                                                                                                                                                                >  1:1658066/1‑239 (MQ=255)
                                        aCCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGgaag                                                                                                                                                                                              >  1:489599/1‑239 (MQ=255)
                                            tcatcGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAAc                                                                                                                                                                                          >  1:2205974/1‑239 (MQ=255)
                                             catcGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACa                                                                                                                                                                                         >  1:819159/1‑239 (MQ=255)
                                                    aGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTATTCGCCGCGCCAATGAGAACATGCCAAAGAGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACa                                                                                                                                                                                         >  1:1884594/1‑232 (MQ=255)
                                                     gTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGgaaga                                                                                                                                                                                             >  1:1663661/1‑227 (MQ=255)
                                                      ttCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGa                                                                                                                                                                                >  1:296346/1‑239 (MQ=255)
                                                        caacagaggTAAACGGCGAGGAGAAAACGTGCGCCGAGCCACCAAAAAACATAACCACGGGGGTGGGTTCCGCCGCGCAAAGGAGAACATGCAAAGGGGGGGCGAAGTTCATTCCCGGGAAAATAACGGGGGCGACAAACAGCGCAAAGAAGGAAGCACCGCCCCGGGTTTGGGGAGGGGAAAACAGGCCGGAGCGAACCCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATa                                                                                                                                                                              <  2:1944701/232‑1 (MQ=255)
                                                                gTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGAc                                                                                                                                                                      >  1:1433455/1‑239 (MQ=255)
                                                                        ccAGGCGGAAACTTGAGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGt                                                                                                                                                              >  1:305932/1‑239 (MQ=255)
                                                                         cAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGt                                                                                                                                                              >  1:727254/1‑238 (MQ=255)
                                                                               gAAAGTTGGGCCGAGGCACCCAGAAACATACCCACGGAGGTGTTTTTCCCCGCGCAAATGAGTACATGCAAAAGGGGGGCGAAGTTCATTTCACGGTAAATAACTGGGGCGACAAACAGCGCATAAAAAGCAGCACCGGCCCCGGCTTGGGGGGGGGTAAACAGGCCGGAGAGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGa                                                                                                                                                        <  2:2205974/238‑1 (MQ=255)
                                                                                       tgCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCACCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAAACCGCCAATAATGCTGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCa                                                                                                                                               >  1:1071197/1‑239 (MQ=255)
                                                                                           ggggCCCCCAAAAACACAACCCGGGGGGGGGTTTCGGCCGCGACAGAGAGTACAGCCAAAGTGGGGGCAAAGGTATTTACAAGGAAAAAAGGGGGGGGAAAAAACAGCGAAAAGAAGGCAGCAGCGGCCCGGGCTGGGGGAGGGGAAAAAGGGCCGGAGGAAACCCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCg                                                                                                                                            <  2:2026171/236‑1 (MQ=255)
                                                                                            aGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGa                                                                                                                                                        >  1:427609/1‑225 (MQ=255)
                                                                                                         aCATAACCGGTGGGGTGTTTTTCCCCGCACAAAGGAGTACAGCCCAAGGGTGGGCAAAGTTCATTCCCGGGTAAATAACGGGGGCAACAAACGGCACAAAAAAAGCCCTCGCGGCCCGGGTTTGGGGGGGAGAAAACGGGGCGGGGAGACCCCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAAAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTgg                                                                                                                              <  2:690224/239‑1 (MQ=255)
                                                                                                             aaCCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCggg                                                                                                                          >  1:1471659/1‑239 (MQ=255)
                                                                                                                  cTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACaaaa                                                                                                                                                                  >  1:827310/1‑194 (MQ=255)
                                                                                                                  cTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGAAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATa                                                                                                                                                                              >  1:2088306/1‑182 (MQ=255)
                                                                                                                   tGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGaaaa                                                                                                                    >  1:1368114/1‑238 (MQ=255)
                                                                                                                    gAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAAtt                                                                                                                   >  1:536307/1‑239 (MQ=255)
                                                                                                                     aGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAAATa                                                                                                                  >  1:880198/1‑239 (MQ=255)
                                                                                                                     aGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAAATa                                                                                                                  >  1:487901/1‑239 (MQ=255)
                                                                                                                          gtTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAATTAAggag                                                                                                             >  1:520938/1‑237 (MQ=255)
                                                                                                                          gtTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGGAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGGGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAAtttaatca                                                                                                             >  1:293727/1‑233 (MQ=255)
                                                                                                                                 ccgcGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCAGGCGGGGGAAAATAAAGAGACTgg                                                                                                        >  1:1011797/1‑237 (MQ=255)
                                                                                                                                    cgcCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGACACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAAATAAGGCGACTGGCCTg                                                                                                    >  1:212182/1‑238 (MQ=255)
                                                                                                                                     cgcAATGAAGTACAGGCAAAATGGGGGAGAAGGTTATCTCCGGGTAAATAACTGCGGAGACAAACAGCGAAAAGAAAGCGGCACCGCCCCCGGTTGGGGGGGGGGAAAACGGGCCGGAGAGAACCCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGc                                                                                                   <  2:1433455/237‑1 (MQ=255)
                                                                                                                                        aaTGAGTACACGAAAGAGTGGGGAGAAGGTACTTGCACGGTAAATAGTGGGGGCAAAAAACAGCGCAGAGAAGGCAGCACCGGCCCCGGCTGGGGGGGGGGAAAACAGCCCGGAGAGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGccacc                                                                                               <  2:489599/239‑1 (MQ=255)
                                                                                                                                                 aggagacgggtgaGGAGAAGGTCATTTCACGGTAAATAGCGGCGGCGAAAAACAGCGCAAAAAAGGCCGCACCGCCCCGGGTTGCGGTGGGGAAAAACAGGCCGGAGAAAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTa                                                                                       <  2:456340/230‑1 (MQ=255)
                                                                                                                                                     caaaaGTGGGGGAAAGGTAATTTCGCGATAAATAACGGGGGGAACAAAGAGCACATAGAAACGCGCACCGCCCCCGGTTTCGGTGGGGGAAAACGGGGCGGAGCAAACCCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTa                                                                                       <  2:1174292/233‑1 (MQ=255)
                                                                                                                                                       aaaGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCCCCAGGTAAGCATc                                                                                 >  1:1633727/1‑239 (MQ=255)
                                                                                                                                                         aGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCCCCAACAGGGAAgcagca                                                                                >  1:667944/1‑238 (MQ=255)
                                                                                                                                                         aGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAAATAAAGCGGCTGGCCTCCCCCCCCCCGGGAAACAtac                                                                                >  2:1438382/1‑236 (MQ=255)
                                                                                                                                                              tgGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAAATAAGGCGACTGGCCTGCCCCCACCAGGTAAGCCTCAGCGtt                                                                           >  1:660183/1‑239 (MQ=255)
                                                                                                                                                                ggaaaaGGTCTTTTCAGGGAAAAAACGTGTGGCGAAAAACAGCGAAGAGAAAGCAGACACGCCCCGGGTTGGGGTGGGGGAAAACGGGCCGGAGCAAACCCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGtt                                                                           <  2:536307/234‑1 (MQ=255)
                                                                                                                                                                ggCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACATGCCTGCCCACACCAGGGAAGCATCAGCGTTg                                                                          >  1:1392378/1‑239 (MQ=255)
                                                                                                                                                                ggCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCCGCGtt                                                                           >  1:221046/1‑238 (MQ=255)
                                                                                                                                                                ggCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGCGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAATTAAGGAGACCGGCCTTCCACCCCCAGGTAAGCCTCAGCGtgg                                                                          >  1:2016176/1‑237 (MQ=255)
                                                                                                                                                                ggCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCACATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAATTAAGGCGACTGGCCTCCCCCCCCCCGGTAAGCATCAGCggt                                                                           >  1:773002/1‑236 (MQ=255)
                                                                                                                                                                       ttCATTTCAAGAAAATAAACGGTGGCGAAAAACGGCGAAAAAAACGCAGCACCGGCCCCGGTTGCGGTGGGGAAAAACAGGCCGGAGCAAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGACCTGCCACCACCAGGTAAGCATCAGCGTTGCGccc                                                                     <  2:660183/236‑1 (MQ=255)
                                                                                                                                                                       gTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACaaaa                                                                                                                                                                  >  1:1238349/1‑144 (MQ=255)
                                                                                                                                                                       gTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGtctggctctt                                                                   >  1:456340/1‑230 (MQ=255)
                                                                                                                                                                         aattaaggggaaaaaaaCTGGGGCAAAAAAGAGCGCAAAAAAGGAAGCCCCGGCCGGGGTGGCGGAGGGGGAAAACGGCCGGGACCAAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGTTGCGCCcatca                                                                 <  2:2152957/228‑1 (MQ=255)
                                                                                                                                                                          aTTTCAGGGTAAAAACTGGGGGCAAAAAACAGCAAAAAAAAAGCAGCAAGGGCCCGGGCTGGGGAGGGGGTAAACAGGCGGGAGCAAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGTTGCGCCcatcat                                                                <  1:217227/239‑1 (MQ=255)
                                                                                                                                                                            ttCACAGTAAAAAGGTGTGGGCAACAAACAGCGACAAGAGAGCAGCACCGCCCCGGGCGGGGGAGGGGGAAAAGGGGCGGGAGCAAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGTTGCGCCcatcatc                                                               <  2:1368114/239‑1 (MQ=255)
                                                                                                                                                                                      aTAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGCCTGGCCTGCCACCACCAGGGTAGCATCAGCGGTGCGCCCATCATCATGCCAaggg                                                     >  1:1869763/1‑238 (MQ=255)
                                                                                                                                                                                         aagcgggggCGACAAACGGCGCAAAGAAAGCAGCACCGGCCCGGGTTGGGGGGGGGAAAAACGGGCGGGAGAGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGTTGCGCCCATCATCATGCCAGGGGCa                                                   <  2:2162464/232‑1 (MQ=255)
                                                                                                                                                                                              tgGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCCCCACCAGGGAAGCCTTAGAGTTGCGCCCCTCATCATGCCAGGGGCAAAGccc                                             >  1:1274183/1‑239 (MQ=255)
                                                                                                                                                                                                 gCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGTTGCGCCCATCATCATGCCAGGGGCAAAGCCCg                                            >  1:684958/1‑238 (MQ=255)
                                                                                                                                                                                                                 aaaaaGCAGCACCGGCCCCGGTGGGGGGGGGGAGAAACGGGCGGGAGAGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGc                                                                                                       <  2:197305/161‑1 (MQ=255)
                                                                                                                                                                                                                      gcagcaACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAATTAAGGCGACTGGCCTGCCACCACCAGGGAAGCATCAGCGTTGCGCCCCTCATCAGGCCAGGGGCAATGCCCCCCACAAACAGCTTACTGATGGa                      >  1:101462/1‑239 (MQ=255)
                                                                                                                                                                                                                           aaCGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAATTAAAGCGACTGGCCTGccacccc                                                                                             >  1:197305/1‑161 (MQ=255)
                                                                                                                                                                                                                             cGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGGAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCg                                                                             >  1:856995/1‑177 (MQ=255)
                                                                                                                                                                                                                                          ggTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGAGACTGGCCTGCCACCACCAGGGAAGCATCAGCGTTGCGCCCATCATCATGCCAGGGGCAAAGCCCGCCATAAACAGCTTACTGATGGAGAATCCGCTGGAAACGCCGa  >  1:1506783/1‑239 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
CAGACGCGGTGAATCGACCAACGGTTGCAGTAAATCAGAAACCATCATCGGCAGTTCAGCAATGGTAATCAGCCAGGCGGAAACTTGTGCCGAGGCCACCAGAAACATAACCACTGAGGTGGTTTTCGCCGCGCCAATGAGTACATGCCAAAGTGTGGCGAAGGTCATTTCACGGTAAATAACTGTGGCGACAAACAGCGCATAGAAAGCAGCAACGGCCCCGGCTTCGGTAGGGGTAAACAGGCCGGAGCGAAACCCGCCAATAATGATGACTGGAAGAAACAGCGCCCAGATACCGGAGACAAAAGAGTGCCAGATTTCCTGCATCGTTGCTTTTTGCTGGCGGGGTAAATTAAGGCGACTGGCCTGCCACCACCAGGTAAGCATCAGCGTTGCGCCCATCATCATGCCAGGGGCAATGCCCGCCATAAACAGCTTACTGATGGATAATCCGCTGGAAACGCCGA  >  NC_007779/3894670‑3895136

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: