Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,237,434 A→G 6.8% L107P (CTG→CCG)  mglC ← methyl‑galactoside transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,237,4340AG 6.8% 170.7 / 3.1 61L107P (CTG→CCG) mglCmethyl‑galactoside transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (2/2);  total (28/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 = 3.50e-01

CCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAACTAAAAGAACGACGTAAATACCGCCCTCTTTCAGG  >  NC_000913/2237134‑2237718
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                            
ccGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCa                                                                                                                                                                                                                                                                                              <  2:588056/301‑1 (MQ=255)
 tttgctcttccgatctCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAg                                                                                                                                                                                                                                                                                             <  2:76670/286‑1 (MQ=255)
                 ccACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAAACAATAAGAATAACCaaagaacggaagagcgac                                                                                                                                                                                                                                                                             >  1:76580/1‑284 (MQ=255)
                  cacaACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAAACAGAGAAATCGGCAtcgtcg                                                                                                                                                                                                                                                                            >  2:512881/1‑301 (MQ=255)
                    caacaCCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACACAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTAACGTTGAGATAAGCGAAAATCAGACCGTTGATCAAACCGATCACCGCCCCAATGGCAAAAACAATCAGAATAACCAGAGAAATCGGCATCGTcccc                                                                                                                                                                                                                                                                          >  1:225171/1‑301 (MQ=255)
                         cccAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCACCGACAAAGCAGTAATAGAGCGAGTTGATGCAATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGAGAAACACGATCACCCCACCAAAGGCCCAGACAAACAGAAACACAAGCGCAAACGACAccgcacacatttc                                                                                                                                                                                                                                                                     >  1:384922/1‑291 (MQ=255)
                             gACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCggg                                                                                                                                                                                                                                                                 <  2:634549/301‑1 (MQ=255)
                                              tCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCACCATTTCCGGGAACACCTTGTTGGCCtt                                                                                                                                                                                                                                                >  1:210064/1‑301 (MQ=255)
                                                 aCGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTGGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATc                                                                                                                                                                                                                                             <  2:86196/301‑1 (MQ=255)
                                                      tAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATgg                                                                                                                                                                                                                                        <  2:272889/301‑1 (MQ=255)
                                                        gAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCAACGCACCCATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGCACACTTTATTGGCGTAACCCATCGAc                                                                                                                                                                                                                                      >  1:197083/1‑301 (MQ=255)
                                                         aaGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACt                                                                                                                                                                                                                                     <  2:250299/301‑1 (MQ=255)
                                                               aTGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCaata                                                                                                                                                                                                                               >  2:119851/1‑301 (MQ=255)
                                                                       gagaCGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTAGATGCAATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGAGCAATCGGCATAGACGACATTTCCGGGAACACATTGTAGGACTAAACCATGCACGACAATAACGTCGcc                                                                                                                                                                                                                       >  1:80314/1‑300 (MQ=255)
                                                                        agaCGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAg                                                                                                                                                                                                                      <  2:225414/301‑1 (MQ=255)
                                                                                aCTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACt                                                                                                                                                                                                              <  2:628751/301‑1 (MQ=255)
                                                                                 cTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCAATg                                                                                                                                                                                                             >  2:180807/1‑301 (MQ=255)
                                                                                   ccccAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGACACCAACACTgcc                                                                                                                                                                                                           >  2:250831/1‑301 (MQ=255)
                                                                                    cccAGGGGGACAAGGCCCGGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAAGTGGCGACGCCCAGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTgccg                                                                                                                                                                                                          <  1:431469/301‑1 (MQ=255)
                                                                                    cccAGCGCGACAAAGCCCTTATCAAAGGTAGAGATGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTgccg                                                                                                                                                                                                          <  2:196871/301‑1 (MQ=255)
                                                                                     caagtgggataaggCCCTGGGCAAAGGT‑GAGGAGGCAATGGTCAAAACAGGAAATGTGTGAGGCCCCGTGAAAGTGGTAATGGGGCGATGTGATGCCAGAGACGAGGATCATCGTGCCGAGCGTGTGAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTgccgc                                                                                                                                                                                                         <  1:350313/288‑1 (MQ=255)
                                                                                          gcgACAAAGCCTTGAGAAAAGTTAGATAAGCCACT‑GTCAAAAACAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGcc                                                                                                                                                                                                    <  2:385418/301‑1 (MQ=255)
                                                                                               aaaGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTAcc                                                                                                                                                                                               <  2:210307/301‑1 (MQ=255)
                                                                                                    ccTGGGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacg                                                                                                                                                                                          <  1:119727/301‑1 (MQ=255)
                                                                                                    ccTGATCAAAGGTAGAGAAGCCACT‑GTTAAAACCAGAAATTGGCCACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacg                                                                                                                                                                                          <  2:183449/301‑1 (MQ=255)
                                                                                                      tgtgctcttccgatctaGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacgac                                                                                                                                                                                        <  2:480586/285‑1 (MQ=255)
                                                                                                       ctgttgcccttcacgacggtcttccgatctccAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGAGGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCGGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacgac                                                                                                                                                                                        <  1:85213/271‑1 (MQ=255)
                                                                                                          aaaaGGTAGGTAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACcagc                                                                                                                                                                                    <  1:597037/300‑1 (MQ=255)
                                                                                                                tAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAg                                                                                                                                                                              <  2:516236/301‑1 (MQ=255)
                                                                                                                  gagaAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAAAt                                                                                                                                                                            >  1:279903/1‑301 (MQ=255)
                                                                                                                      aGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgccgacagatcggaagagcgtc                                                                                                                                                                        >  1:479985/1‑285 (MQ=255)
                                                                                                                           cT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTACCGCCAGCCCTACCTGACGACCAGCAGAAAAATCAGTcccac                                                                                                                                                                   >  2:184882/1‑299 (MQ=255)
                                                                                                                           actgggaaaaCCGAAAATTGGGGTCGCCCCGACAAGGTCGTAAGAGGGCTAGTTGATGCCATTGACGGTGATATGCGTGCCGAGCGTGGTAATTAACGGCGTCACTGTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATGTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTcccc                                                                                                                                                                    <  1:434477/295‑1 (MQ=255)
                                                                                                                              gTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGCGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGCACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGCCAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCCTGGACTGCACAAATGACGCAGCCCCCACTGCCGACAGCCATCCCTgacgacacgacgaacaaacgacacccccg                                                                                                                                                                 >  1:418669/1‑278 (MQ=255)
                                                                                                                               tCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTggg                                                                                                                                                                <  2:604589/301‑1 (MQ=255)
                                                                                                                                aaaaaCCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCAGCGTGCCGAGCGTGGTAATGAACGGCGTCACGTCGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGt                                                                                                                                                               <  2:119640/300‑1 (MQ=255)
                                                                                                                                     ccAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCGGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacgacagatcggaagagcaaaagtctgaactccaa                                                                                                                                                          >  2:85337/1‑271 (MQ=255)
                                                                                                                                           aTTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAAcc                                                                                                                                                    >  2:467746/1‑301 (MQ=255)
                                                                                                                                                                tCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGataa                                                                                                                               <  2:280206/301‑1 (MQ=255)
                                                                                                                                                                  ccgttctagCTAGTTGATGCCATTGACGATGATCCTCTTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGataata                                                                                                                             <  1:69721/294‑1 (MQ=255)
                                                                                                                                                                       ttttccgatctttGCCATAGACGATGATCATCGTGCCTAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACg                                                                                                                        <  2:15346/289‑1 (MQ=255)
                                                                                                                                                                         agCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGataataagatcgg                                                                                                                      >  2:69811/1‑295 (MQ=255)
                                                                                                                                                                                  aTGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGGCCCCTGGGTGAAAATTAACCCAGCCAccccacacacgaaaacaaaacgagctcggaacc                                                                                                             >  1:15345/1‑272 (MQ=255)
                                                                                                                                                                                      cATAGACGTTGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACTTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGt                                                                                                         <  1:260016/301‑1 (MQ=255)
                                                                                                                                                                                      cATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTgggg                                                                                                         >  2:260284/1‑300 (MQ=255)
                                                                                                                                                                                      cATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGt                                                                                                         >  2:151909/1‑301 (MQ=255)
                                                                                                                                                                                       aTAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTg                                                                                                        <  2:542164/301‑1 (MQ=255)
                                                                                                                                                                                           acggggatcatctgtccgtggggggTAATGAACGGTGTCGCGTTGAGCTAAGCGATAATCAGACCTTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCTGCGCACTCGGCATCGTCGCCATTTCCGGGAACACTTTGTTTGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCTGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGaa                                                                                                    <  1:244234/278‑1 (MQ=255)
                                                                                                                                                                                               tgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGAATGGGTGAGAatata                                                                                                >  2:61688/1‑300 (MQ=255)
                                                                                                                                                                                                            gAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGTCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAAc                                                                                   <  2:418170/301‑1 (MQ=255)
                                                                                                                                                                                                                     aaTGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTACGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGACGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATAGGTCCCATGGGTGACAATTAAACCTGCCACACCGAACGCGATCATACTACGAACCAATGACTGGGTCAGACAATTAATTAAAATCAAAAAACTTaaa                                                                          >  2:418461/1‑301 (MQ=255)
                                                                                                                                                                                                                       tGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGAACCGATGACTGGGTGAGAATATTACTTAAGTTAAACAGACTTaaaaa                                                                        >  1:139035/1‑301 (MQ=255)
                                                                                                                                                                                                                               gTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCGGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAAACTTAAAACTGTTGGGt                                                                >  1:619559/1‑301 (MQ=255)
                                                                                                                                                                                                                                      tGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTCACCCTGCCACACCGACCGCGATAATACTACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGCCCTGGaa                                                         >  1:461155/1‑301 (MQ=255)
                                                                                                                                                                                                                                            aaGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAAACTTAAAAATGTTGGGTCCCGGAaaaaaaa                                                   >  1:28240/1‑300 (MQ=255)
                                                                                                                                                                                                                                                 aTAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATca                                              >  2:427543/1‑300 (MQ=255)
                                                                                                                                                                                                                                                 aTAATCAGACCGTGGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCg                                              <  2:461719/301‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGc                                             <  2:139160/301‑1 (MQ=255)
                                                                                                                                                                                                                                                           cGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAAc                                    <  2:197326/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                aTCAAATCGATCAACGCACCAATGGTACAGACAATCAGAATAACCAGTGCAATCGTCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAACTaaaa                               <  1:467182/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAACTAAAGGAAAGACGTAAATa                 >  2:243314/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                               ccACCAATGGCACAGACAATCAGAATAACCGGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAACTAAAAGAACGACGTAAATAc                <  2:620417/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             gACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATATGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTAGGTCCTGGAAAATAATAATCGCCAGCAAAACTAAAAGAACGACGAAAAAACCGCCCTCTTTCaag  >  1:50496/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                            
CCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACT‑GTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAACTAAAAGAACGACGTAAATACCGCCCTCTTTCAGG  >  NC_000913/2237134‑2237718

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