Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,777,745 A→T 40.7% I312F (ATT→TTT)  Y75_RS18735 → PTS transporter subunit EIIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,777,7450AT40.7% 37.4 / 19.4 70I312F (ATT→TTT) Y75_RS18735PTS transporter subunit EIIC
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base T (4/13);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 7.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCGCAACCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAT  >  NC_007779/3777510‑3777980
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
cTTCATCTACGGGCAGATTATCTCTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCTGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAACTCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTAACTGCCGCACCGGAAAATCGGGGAAAAGTGGCGGGCCTGATGtttt                                                                                                                                                                                                                                                       >  1:1760254/1‑238 (MQ=255)
cTTCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGagt                                                                                                                                                                                                                                                         >  1:1927651/1‑236 (MQ=255)
 ttCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGGttt                                                                                                                                                                                                                                                       >  1:1731443/1‑237 (MQ=255)
 ttCATCTACGGAAAGTTTATCTTTGGTCCGGCAGCTGTTGCAGGCGGCATTCAGATGTACAGGACGCAGCATCTGCAAGAGTTCAGTATGAGCGCCGAGCCGCTGCAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGCGGGCATTTCTTTAGCGATGTACTTCACTGCCGAACCGGACAATTGGGTAAAAGTGGCAGGACTGCTGAGTc                                                                                                                                                                                                                                                      >  2:1400055/1‑238 (MQ=255)
  tccccttcGGAAAGTTTATCTTGTGGCCGGCGTCTGTGGAAGGTGGCTTTCAGATGTACTGGTCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTccc                                                                                                                                                                                                                                                   <  2:1908452/232‑1 (MQ=255)
  tCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTccc                                                                                                                                                                                                                                                   <  1:1953056/239‑1 (MQ=255)
  gcATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTccc                                                                                                                                                                                                                                                   <  1:817661/238‑1 (MQ=255)
   cATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGGATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGGATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCg                                                                                                                                                                                                                                                  <  1:1877069/239‑1 (MQ=255)
   cATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTACCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCCTTCTGGCTccc                                                                                                                                                                                                                                                   >  1:1740184/1‑238 (MQ=255)
    aTCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCGc                                                                                                                                                                                                                                                 <  1:1698618/239‑1 (MQ=255)
     tCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCGTTCTGGtttcacc                                                                                                                                                                                                                                                   >  2:18147/1‑233 (MQ=255)
     tCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCGTGGTGGTTtccccc                                                                                                                                                                                                                                                >  1:324073/1‑236 (MQ=255)
      cTAGGGACGGTTTTGCTTTGGCCCGGCATCGGGTAAAGGCGGCGTTCCGGTGTAGGGGGCGCAGGATCTGAAAGGGTTCGGTTTGACCCCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCGCaa                                                                                                                                                                                                                                               <  2:1815682/239‑1 (MQ=255)
          ggACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGCTGATTCCCGCCAccct                                                                                                                                                                                                                                           >  1:1447431/1‑237 (MQ=255)
           gACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGTTGGTTACCCCCAACCtt                                                                                                                                                                                                                                          >  1:389358/1‑239 (MQ=255)
           gACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGTTGTTTTTCGCCACCCTa                                                                                                                                                                                                                                        >  1:44810/1‑239 (MQ=255)
                  tATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTTGTGTTTCCCCCAACCCTAAACCCCa                                                                                                                                                                                                                              >  1:1716562/1‑239 (MQ=255)
                  tATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTGGGTGACTTCCGCAGCCtccaaagcc                                                                                                                                                                                                                                   >  2:1606761/1‑230 (MQ=255)
                   aTCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGCTGAGTTCCGCAAACTTTACccccca                                                                                                                                                                                                                             >  1:223635/1‑237 (MQ=255)
                        tGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGGTGTTGATTCCCCCCAACTAAACcccccaagtgg                                                                                                                                                                                                                        >  1:159748/1‑232 (MQ=255)
                        tGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTTCTGGTTCCCGCCCCCCTAACCCGCATGCtgg                                                                                                                                                                                                                      >  1:540138/1‑239 (MQ=255)
                          gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGTTTTCTTATTCCCGCAACCCTAACCCCCACGCtggtg                                                                                                                                                                                                                    >  1:1026356/1‑239 (MQ=255)
                          gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGTTGTTTCCCGCCACCCTAACCCCCATGCtgggg                                                                                                                                                                                                                    >  1:280316/1‑237 (MQ=255)
                          gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGCTTATTTCCCCAACCCTAAACGCCCTGCtggtg                                                                                                                                                                                                                    >  1:422289/1‑239 (MQ=255)
                          gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGTTTGTTTCCGCCACCCTAACCCCCCCGCtgggg                                                                                                                                                                                                                    >  1:1945113/1‑237 (MQ=255)
                          gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCCTTCGTATTTCCCCCACCCTAACAGCCaccccggtg                                                                                                                                                                                                                      >  1:2183122/1‑231 (MQ=255)
                          gTCCGGCAGCTGTGGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCGCAACCTTAACCGCCATGCtggt                                                                                                                                                                                                                     <  1:1187568/238‑1 (MQ=255)
                               gCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGCTGGTTTCCCCCACCTTAAACGCCCTTCTGGTGGGGAt                                                                                                                                                                                                               >  1:165701/1‑239 (MQ=255)
                                              ggcATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACGGGAAAATCGGGTAAAAGTGGCGGGCTTGTTGCTTCCCGCCACCTTAACCCCCCTGCCGGGGGGGGTTACCGCACCGCTggg                                                                                                                                                                                             >  1:2170315/1‑238 (MQ=255)
                                                cATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTGTGGGTTTCCCCCACCCTACACCCCACTCTGGGGGGGATTTACCGAAccccgggga                                                                                                                                                                                             >  2:872430/1‑232 (MQ=255)
                                                cATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTTCTTATTCCCCCAACCCTTACCCCCACGCTGGTGGGGATTACCGGAACGCCTGGaa                                                                                                                                                                                            >  1:316163/1‑239 (MQ=255)
                                                cATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGGTTTCCTCCCCCCTAACCGCCACGCTGGGGGGGAGTAGCGAAACGCTGGAAt                                                                                                                                                                                           >  1:1577780/1‑239 (MQ=255)
                                                cATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGATGGTTGATTCCCCCCCCCCTAAACGCGAAACTAGTGGGGATTACAGGAccccctgaga                                                                                                                                                                                            >  1:1105351/1‑231 (MQ=255)
                                                  ttCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCg                                                                                                                                                                                                                                                  <  1:275786/192‑1 (MQ=255)
                                                   tggtgtGTACTGGGGCACGAACTGCAAGGGTTTATGTTGGAGCGCCGAGCCGTTAAAATCGTTGTTCGCGGAGGGGGGTTTTCCCCTGCACGGTAACTAAAAAATCTTTGGGGCGGGGGGCTTTTCTTTAGCATTGTTCTTCACTGCCCCCCCGGAAAATCGGGTAAAGGGGGCGGGCTTGCTGATCCCCGCAACCTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTa                                                                                                                                                                                        <  2:540138/234‑1 (MQ=255)
                                                            cTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTTTCGATTCCCGCACCCCTAACCGCCCTGCTGGGGGGGAGTAACGAAACGCTGGGATTTACCTTCCTGt                                                                                                                                                                               >  1:42228/1‑239 (MQ=255)
                                                            cTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTTCTGGTTTCCGCAACCTTAACCGCCACGCTGGTGGGGAGTAACGAACCGCCGGGATTTATCTTCCCGt                                                                                                                                                                               >  1:425692/1‑239 (MQ=255)
                                                            cTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGGTTATTCCCCCCACCCTCAACGCCATGCCGGTGGGGATTAACCAAACCCTGGAATTTAACTTCCttg                                                                                                                                                                               >  1:1216887/1‑237 (MQ=255)
                                                            cTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCCTGGTGATTCCCGCCACCCTAACCCCCATGCTGGGGGGGATTACCGAACCGCCGGAATTTTCCTTCCTgg                                                                                                                                                                               >  1:1540425/1‑238 (MQ=255)
                                                               ggCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTTCTTGTTTCCCCCCCCTTTACCCCCACGCTGGGGGGGGTTACCCAACCGCCGGGAATTTCCTTCCTGtt                                                                                                                                                                              >  1:702898/1‑237 (MQ=255)
                                                               ggCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCGTTGCGATTCCCGCCACCTTAAACGCCCATCGTGGGGGGAATACCCGAACGCTTGGAATTTCCTTCCTTTTCa                                                                                                                                                                            >  1:2164473/1‑239 (MQ=255)
                                                               ggCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTTTTGGTTCCCGCAACCTTTACCGCCATTCAGGGGGGGGATAACGGACAGCTGGCAGTTATCTTCCTCtttt                                                                                                                                                                            >  1:1952751/1‑237 (MQ=255)
                                                               ggCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGGTTCTCCCCACCTTTAACCCCCACCTGGTGGGGATTACCGCAACCCCGGAAATTAc                                                                                                                                                                                       >  1:1603219/1‑228 (MQ=255)
                                                                 cgcAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTTCTGGTTTCCCCAACCCTTACCGCCATGGCGGGGGGGAGTTCCGGACCGCCGGGACTTTACTTTCTTTTTctt                                                                                                                                                                           >  1:2082841/1‑237 (MQ=255)
                                                                 cgcAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGGTTAATCCCGCCACCTTAACCGCCACGCATGGGGGGATTTCCGAAACCCCTGGAATTAACCTCCCTGtcatt                                                                                                                                                                           >  1:129567/1‑235 (MQ=255)
                                                                      caTCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCGTGGTGATTTCCGCCACCTTAACCGCCCGACCGGGTGGGGATACCGAAACGCTGGGCGTTTCCTTCCTCTTCCTTTCat                                                                                                                                                                     >  2:1313934/1‑237 (MQ=255)
                                                                      caTCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCCTGCTGTTTCCCCCCACCCTAACCGCCCTGCTGGGGGGGATTTACGAACCGCCGGGATTTTCCTTCCTTTTCAtttttc                                                                                                                                                                     >  1:1627163/1‑235 (MQ=255)
                                                                         cTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCCTGCTGATTTCCGCCACCTTAACCGCCCAGCTGGTGGGGATTTCCGAAACCGCGGGAATTTCCTTCCTGTTTATTTTCACCGt                                                                                                                                                                  >  1:183196/1‑239 (MQ=255)
                                                                               gtgTTCAGTGGGGGCGCGGGGCGGTTAAAATGGTTGTCCCGGGAGGGGGGTTTTGCCCGGCCGGTTAACACAAAATTTTTGGGGGCGGTGGGATTTTCTTGGGGGTTGTACTTCCCGCCCGCACGGAAAACTGGGGAAAAGGTGGGGGTCTTGCTGATTCCCGCACCTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGt                                                                                                                                                            <  2:1447431/236‑1 (MQ=255)
                                                                                            aGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTGGCTGGTTTCCGCCACCCTATACGCCCTCCTGGGGGGGGTTAGCCAACAGCTGGGATTTTTATTCCCGTTTCTTTTACCGTTTCTGTTTTCTGTacacc                                                                                                                                                >  2:1500372/1‑237 (MQ=255)
                                                                                                      cgttGAATTCTTTGTCCCGGGAGGGGGTTTTTCCCCTGACGGTTAACTAAAAAGTTTTGGGGCGCGGGGGTTTTTTTTTAGCGTTTTTCTTCACTGCCCCACCGAAAAAGGGGGTAAAGGGGGCGGGTTTGCTGTTCCCCGAACCTTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTggc                                                                                                                                      <  2:2164473/235‑1 (MQ=255)
                                                                                                      ccgttaaaagtttgttCCGCGGAAGGGGTTTTTGCCCTGCAGGGTACCTCAAAATTCTTGGGGCCGGGGGGCTTTTTTTTGGCGATTTCCTCCCCCGCCGCCCCGGAAAATGGGGTAAAGGGGGGGGTTTTGCTTTTCCCCGAAACTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTggc                                                                                                                                      <  2:422289/228‑1 (MQ=255)
                                                                                                                 tttttCCCGGAGGGGGGTTTTGCCCTCCACGGTAACTAAAAAATTTTGGGGGCGGGGGGCTTTTTTTTAGCGATGTCCTTCACTGCCGCACCGAAAAATGGGAAAAAGGTGGGGGGTTTGCTGATTCCCGCACCCTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTc                                                                                                                          <  1:344393/236‑1 (MQ=255)
                                                                                                                 gttgtCCGGGAGGGCGTTTTTCCCCTGCGCGAAACCTAAAAATTTTTGGGTGCGTGGGGCTTTCTTTTGGTGATTTACTTCTCCCCCCCCCGGAAAATTGGGGAAAAGGGGGGGGGCTTGCTGTTCCCCCCACCTTAAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTc                                                                                                                          <  2:1658941/235‑1 (MQ=255)
                                                                                                                   gttCCCGGAGGGCGTTTTTGCCACGCGCATTATCACAAAATTCTTTGGTGCGGGGGGTATTTTTTTGGTGATTTTCTCCCCTGCCCCCCGGAAAACGGGGAAAAAGGGGGGGGGCTTGCGTTTCCCCGCACCTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGa                                                                                                                        <  2:425692/239‑1 (MQ=255)
                                                                                                                     ttgccgtgaGGGGGTTTTTGCCTGGCCGGTAACTTAAAAATTTGGTGGTGCGGGGGGATTTTTTTAAGGTGTTTTCTTCCCCCCCCCACCGAAAAGTGGGATAAAAGTGGGGGGTTTGCTACTCCCGGAACCTTAATCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGa                                                                                                                        <  2:2073202/231‑1 (MQ=255)
                                                                                                                     tCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGGTTTTTGATTTCCCCACACCTAAACGCCCTGCTGGTGGGGATTACCGAACCCCCGGGACTTTCCCTCCTCTTCATTTTACCGTTTCTGTTTGTCGTATGCGCCGCGGTGGCGGCCTccactgcgga                                                                                                                       >  2:1708344/1‑229 (MQ=255)
                                                                                                                       ccgggAGGGCTGTTTTGCCTTGCACGGTAACTAAAAATTCTTGGGTCCGGGGGGCATTTTTTGAGGGATTTATTCCACTCCCGCACCGAAAAATGGGTAAAAAGGGGGGGGTTTGCTGATTCCCCCAACTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCg                                                                                                                     <  1:431034/236‑1 (MQ=255)
                                                                                                                             ggtcGGTTTTCCCTGCGACGTAACATAAAAATTTTTGGGTGCGGGGGGTTTTTTTTGAGGGTTTTCCTCCCCTGCCGCACGGAAAAAGGGGGAAAAGGGGGGGGTTTGGTTACTCCCCGCAACCTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATg                                                                                                                <  2:134994/235‑1 (MQ=255)
                                                                                                                                   ttttGCCCTGCACGGTAAACCAAAAATTTTTGTGGCCGTGTGTCTTTCTTTTGGGGTTTGTTTCCCCGCCGGCCAGGAAAAGTGGAGAAAAGGTGGGGGTTTTGTGTTTCCCCGAACCTTTATCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTAtct                                                                                                           <  2:1585012/238‑1 (MQ=255)
                                                                                                                                                ggTACCTCAAAAATTTTTGGTGCCGTGGGCTTTTTTTTAGCGTTGTACTTCACTCCCGCCCCGAAAAATCGGTTAAAAGTGGCGGGTTTGCTGTTCCCCGAACCTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAAt                                                                                                                             <  1:464059/206‑1 (MQ=255)
                                                                                                                                                       tCAAAAATCTTTGGTGCGGGGGTCTTTCTTTTAGGGTTTTCTTTCCCTGCCCCCCGGAAAATTGGGGTAAAGGTGGCGGTTTTCCTTTTTCCCGCACCTTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACAt                                                                                      <  1:639291/239‑1 (MQ=255)
                                                                                                                                                                   gttgggtgggggggTATTTTTTTGGCGATTTCCCTCACGCCAGCACGTGAAAGGGGGGTAAGGGGGGGGGTCTTGTGGATCCCCGCACCCTAGCCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAgg                                                                             <  2:1995064/228‑1 (MQ=255)
                                                                                                                                                                         gcggggggtttttCTTTAGCGTTGTCCTTCACCGCCCCCCGGAAAAACGGGGTAAAAGGGGCGGGCTTGCTTATTCCCGCAACTTTAGCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGa                                                                        <  1:110014/228‑1 (MQ=255)
                                                                                                                                                                             ggggtttttCTTTAGCGTTGTCCTTCACTCCCCCCCCGAAAATGGGGAAAAATGTGGCGGGTTGGCTGTTCCCCGCAACTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGa                                                                    <  1:551831/232‑1 (MQ=255)
                                                                                                                                                                                                    ctccccTGCCGCCCCGGAAATTCGGGTAAAAGTGGCGGTTTTTCTTTTCCCCGCATCTTTACCCGCCATGCTGGTGGGAATTACCGATCCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGa                                             <  1:1606761/234‑1 (MQ=255)
                                                                                                                                                                                                    cTTCCCTGCCGCACGGGAAACTGGGTTAAAAGTGGGGGTCTTGCTTTTCCCCGCATTCTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGa                                             <  1:1642169/239‑1 (MQ=255)
                                                                                                                                                                                                          tccgggcccGGAAAACGGGGTAAGAGTGGGGGGTTGCCGTTGCCCCGACCCTTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGa                                       <  1:130514/232‑1 (MQ=255)
                                                                                                                                                                                                                     aaaaGCGGGTAAATGGGGCGTGCTTGCTGTTTCCCGCACCTTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAAc                            <  1:772671/239‑1 (MQ=255)
                                                                                                                                                                                                                            gggttAAATGTGGGGGGTTTGTTGTTCCCCCCACCTTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCAcgc                       <  1:1192700/235‑1 (MQ=255)
                                                                                                                                                                                                                            gggatAAAGGGGGGGGGTTTCCTGTTTCCCACACCTTAGCCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACGTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGTCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGGTGTTCAGAAACCAcgc                       <  1:579711/235‑1 (MQ=255)
                                                                                                                                                                                                                                 aaagGTGGCGGTCTTGCTGATTCCCGCAACCTTACCCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGata                  <  1:1287246/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    ctaTTCCCGCCACCTTAAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAt  <  1:324536/237‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
CTTCATCTACGGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTTTGCCCTGCACGGTAACTCAAAAATCTTTGGTGCCGTGGGCATTTCTTTAGCGATGTACTTCACTGCCGCACCGGAAAATCGGGTAAAAGTGGCGGGCTTGCTGATTCCCGCAACCTTAACCGCCATGCTGGTGGGAATTACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAT  >  NC_007779/3777510‑3777980

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