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

CCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAACTAAAAGAACGACGTAAATACCGCCCTCTTTCAGG  >  NC_000913/2237134‑2237718
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                            
ccGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCa                                                                                                                                                                                                                                                                                              <  2:588056/301‑1 (MQ=255)
 tttgctcttccgatctCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAg                                                                                                                                                                                                                                                                                             <  2:76670/286‑1 (MQ=255)
                 ccACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAAACAATAAGAATAACCaaagaacggaagagcgac                                                                                                                                                                                                                                                                             >  1:76580/1‑284 (MQ=255)
                  cacaACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAAACAGAGAAATCGGCAtcgtcg                                                                                                                                                                                                                                                                            >  2:512881/1‑301 (MQ=255)
                    caacaCCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACACAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTAACGTTGAGATAAGCGAAAATCAGACCGTTGATCAAACCGATCACCGCCCCAATGGCAAAAACAATCAGAATAACCAGAGAAATCGGCATCGTcccc                                                                                                                                                                                                                                                                          >  1:225171/1‑301 (MQ=255)
                         cccAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCACCGACAAAGCAGTAATAGAGCGAGTTGATGCAATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGAGAAACACGATCACCCCACCAAAGGCCCAGACAAACAGAAACACAAGCGCAAACGACAccgcacacatttc                                                                                                                                                                                                                                                                     >  1:384922/1‑291 (MQ=255)
                             gACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCggg                                                                                                                                                                                                                                                                 <  2:634549/301‑1 (MQ=255)
                                              tCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCACCATTTCCGGGAACACCTTGTTGGCCtt                                                                                                                                                                                                                                                >  1:210064/1‑301 (MQ=255)
                                                 aCGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTGGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATc                                                                                                                                                                                                                                             <  2:86196/301‑1 (MQ=255)
                                                      tAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATgg                                                                                                                                                                                                                                        <  2:272889/301‑1 (MQ=255)
                                                        gAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCAACGCACCCATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGCACACTTTATTGGCGTAACCCATCGAc                                                                                                                                                                                                                                      >  1:197083/1‑301 (MQ=255)
                                                         aaGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACt                                                                                                                                                                                                                                     <  2:250299/301‑1 (MQ=255)
                                                               aTGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCaata                                                                                                                                                                                                                               >  2:119851/1‑301 (MQ=255)
                                                                       gagaCGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTAGATGCAATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGAGCAATCGGCATAGACGACATTTCCGGGAACACATTGTAGGACTAAACCATGCACGACAATAACGTCGcc                                                                                                                                                                                                                       >  1:80314/1‑300 (MQ=255)
                                                                        agaCGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAg                                                                                                                                                                                                                      <  2:225414/301‑1 (MQ=255)
                                                                                aCTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACt                                                                                                                                                                                                              <  2:628751/301‑1 (MQ=255)
                                                                                 cTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCAATg                                                                                                                                                                                                             >  2:180807/1‑301 (MQ=255)
                                                                                   ccccAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGACACCAACACTgcc                                                                                                                                                                                                           >  2:250831/1‑301 (MQ=255)
                                                                                    cccAGGGGGACAAGGCCCGGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAAGTGGCGACGCCCAGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTgccg                                                                                                                                                                                                          <  1:431469/301‑1 (MQ=255)
                                                                                    cccAGCGCGACAAAGCCCTTATCAAAGGTAGAGATGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTgccg                                                                                                                                                                                                          <  2:196871/301‑1 (MQ=255)
                                                                                     caagtgggataaggCCCTGGGCAAAGGTGAGGAGGCAATGGTCAAAACAGGAAATGTGTGAGGCCCCGTGAAAGTGGTAATGGGGCGATGTGATGCCAGAGACGAGGATCATCGTGCCGAGCGTGTGAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTgccgc                                                                                                                                                                                                         <  1:350313/288‑1 (MQ=255)
                                                                                          gcgACAAAGCCTTGAGAAAAGTTAGATAAGCCACTGTCAAAAACAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGcc                                                                                                                                                                                                    <  2:385418/301‑1 (MQ=255)
                                                                                               aaaGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTAcc                                                                                                                                                                                               <  2:210307/301‑1 (MQ=255)
                                                                                                    ccTGGGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacg                                                                                                                                                                                          <  1:119727/301‑1 (MQ=255)
                                                                                                    ccTGATCAAAGGTAGAGAAGCCACTGTTAAAACCAGAAATTGGCCACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacg                                                                                                                                                                                          <  2:183449/301‑1 (MQ=255)
                                                                                                      tgtgctcttccgatctaGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacgac                                                                                                                                                                                        <  2:480586/285‑1 (MQ=255)
                                                                                                       ctgttgcccttcacgacggtcttccgatctccAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGAGGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCGGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgacgac                                                                                                                                                                                        <  1:85213/271‑1 (MQ=255)
                                                                                                          aaaaGGTAGGTAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACcagc                                                                                                                                                                                    <  1:597037/300‑1 (MQ=255)
                                                                                                                tAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAg                                                                                                                                                                              <  2:516236/301‑1 (MQ=255)
                                                                                                                  gagaAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAAAt                                                                                                                                                                            >  1:279903/1‑301 (MQ=255)
                                                                                                                      aGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTgccgacagatcggaagagcgtc                                                                                                                                                                        >  1:479985/1‑285 (MQ=255)
                                                                                                                           cTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTACCGCCAGCCCTACCTGACGACCAGCAGAAAAATCAGTcccac                                                                                                                                                                   >  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)
                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                            
CCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGTCGCAGCCACCACTGCCGCCAGCCCTACCTGACGACCAGCAGAAAGATCGGTCCCCTGGGTGACAATTAACCCTGCCACACCGAGCGCGATAATAATACGCACCGATGACTGGGTGAGAATATTACTTAAGTTCAACAGACTTAAAAATGTTGGGTCCTGGAAAATAATAATCGCCAGCAAAACTAAAAGAACGACGTAAATACCGCCCTCTTTCAGG  >  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: