Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,852,854 G→T 100% P801P (CCG→CCT mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,852,8540GT97.9% 743.2 / 5.8 236P801P (CCG→CCTmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (112/119);  total (115/121)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 236 total aligned reads displayed.

GCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTG  >  CP009273/2852709‑2852998
                                                                                                                                                 |                                                                                                                                                
gCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTaa                                                                                                                                                <  1:66552/148‑1 (MQ=255)
 cGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAc                                                                                                                                               >  1:3288255/1‑148 (MQ=255)
        aTTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccgccg                                                                                                                                        >  2:2573836/1‑148 (MQ=255)
        aTTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccgccg                                                                                                                                        >  1:2305648/1‑148 (MQ=255)
          tGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCt                                                                                                                                      >  2:1888802/1‑148 (MQ=255)
            ccTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTAc                                                                                                                                    >  2:3543017/1‑148 (MQ=255)
                  tGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGt                                                                                                                              >  1:1390951/1‑148 (MQ=255)
                   gCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTg                                                                                                                             >  1:1850614/1‑148 (MQ=255)
                      caGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCACGTGGAt                                                                                                                          <  2:1481625/148‑1 (MQ=255)
                      caGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGAt                                                                                                                          <  2:152902/148‑1 (MQ=255)
                      caGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGAt                                                                                                                          <  2:912484/148‑1 (MQ=255)
                       aGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTTCAGCTCTGTCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATg                                                                                                                         <  1:3235588/148‑1 (MQ=255)
                       aGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATg                                                                                                                         <  1:1226586/148‑1 (MQ=255)
                       aGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGAACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATg                                                                                                                         <  1:113569/148‑1 (MQ=255)
                        gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATgg                                                                                                                        <  1:2451021/148‑1 (MQ=255)
                                gATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCaaa                                                                                                                >  1:179660/1‑148 (MQ=255)
                                gATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCaaa                                                                                                                >  2:3269484/1‑148 (MQ=255)
                                  tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATg                                                                                                              <  2:248427/148‑1 (MQ=255)
                                   ggCGCGGCGATCAAAAGCTACGGCCTTGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGt                                                                                                             <  1:2381094/148‑1 (MQ=255)
                                    gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTc                                                                                                            <  1:2684558/148‑1 (MQ=255)
                                    gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTc                                                                                                            >  2:569310/1‑148 (MQ=255)
                                    gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTc                                                                                                            <  1:899803/148‑1 (MQ=255)
                                    gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTc                                                                                                            >  1:2076726/1‑148 (MQ=255)
                                      gcggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCtt                                                                                                          <  2:3027936/148‑1 (MQ=255)
                                      gcggagAGCAAATGCTACGGCATGGCGGTTGCATCTCTGGCAGGCGTGCCAAAAGAGTTTTTTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCtt                                                                                                          <  1:1340106/148‑1 (MQ=255)
                                        ggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTg                                                                                                        >  1:2438095/1‑148 (MQ=255)
                                         gcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAGAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGc                                                                                                       <  1:537656/148‑1 (MQ=255)
                                         gcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGc                                                                                                       <  1:3501102/148‑1 (MQ=255)
                                          cgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCt                                                                                                      <  1:929190/148‑1 (MQ=255)
                                            aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                                                                                    <  2:1984434/148‑1 (MQ=255)
                                                aaaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGt                                                                                                <  2:2954384/148‑1 (MQ=255)
                                                 aaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTa                                                                                               >  1:875316/1‑148 (MQ=255)
                                                  aGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTAc                                                                                              >  2:1031666/1‑148 (MQ=255)
                                                  aGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTAc                                                                                              <  2:2124384/148‑1 (MQ=255)
                                                  aGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTAc                                                                                              <  2:446055/148‑1 (MQ=255)
                                                  aGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTAc                                                                                              <  1:276820/148‑1 (MQ=255)
                                                    cTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCa                                                                                            >  1:2819569/1‑148 (MQ=255)
                                                      aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCaga                                                                                          >  1:3500463/1‑148 (MQ=255)
                                                          ccTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAgaagaa                                                                                      <  2:2209639/148‑1 (MQ=255)
                                                            tGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAAc                                                                                    <  2:1129995/148‑1 (MQ=255)
                                                             ggCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACt                                                                                   >  2:1020520/1‑148 (MQ=255)
                                                             ggCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACt                                                                                   >  2:425500/1‑148 (MQ=255)
                                                              gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACtt                                                                                  <  2:1686900/148‑1 (MQ=255)
                                                              gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACtt                                                                                  >  2:1285662/1‑148 (MQ=255)
                                                                 gTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGc                                                                               <  2:3140981/148‑1 (MQ=255)
                                                                  ttGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGcc                                                                              >  2:786028/1‑148 (MQ=255)
                                                                  ttGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGcc                                                                              >  2:610351/1‑148 (MQ=255)
                                                                      aGCTCTGGCAGGCGTGCCAAAAGAGGTGATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTAGCCTGCg                                                                          >  1:3041194/1‑148 (MQ=255)
                                                                          ctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                      <  1:1652303/148‑1 (MQ=255)
                                                                             gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAg                                                                   >  1:736533/1‑148 (MQ=255)
                                                                             gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAg                                                                   <  1:921843/148‑1 (MQ=255)
                                                                             gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAg                                                                   >  1:2161603/1‑148 (MQ=255)
                                                                                 gCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTg                                                                                                        <  2:342009/107‑1 (MQ=255)
                                                                                 gCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTg                                                                                                        >  1:342009/1‑107 (MQ=255)
                                                                                  cGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTg                                                              <  2:2083922/148‑1 (MQ=255)
                                                                                        aaaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAAt                                                        <  1:3461877/148‑1 (MQ=255)
                                                                                         aaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATc                                                       <  2:582880/148‑1 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                      <  1:995568/148‑1 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                      >  1:1044604/1‑148 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                      >  1:1618668/1‑148 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                      >  1:2827336/1‑148 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                      >  1:2952938/1‑148 (MQ=255)
                                                                                           agagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCtt                                                     >  2:2224191/1‑148 (MQ=255)
                                                                                              ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                  <  2:2227685/148‑1 (MQ=255)
                                                                                              ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                  >  1:292456/1‑148 (MQ=255)
                                                                                              ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                  <  2:1166662/148‑1 (MQ=255)
                                                                                              ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                  >  1:428380/1‑148 (MQ=255)
                                                                                               gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                 >  1:3055606/1‑148 (MQ=255)
                                                                                                ttattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATcc                                                >  2:3315327/1‑148 (MQ=255)
                                                                                                  attaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCgg                                              <  2:662597/148‑1 (MQ=255)
                                                                                                     aaGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAtt                                           <  2:721146/148‑1 (MQ=255)
                                                                                                     aaGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAtt                                           <  1:1575589/148‑1 (MQ=255)
                                                                                                      aGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtc                                          >  1:2503683/1‑148 (MQ=255)
                                                                                                       gcgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtca                                         <  1:3164642/148‑1 (MQ=255)
                                                                                                       gcgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtca                                         >  1:495223/1‑148 (MQ=255)
                                                                                                         gcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcact                                       <  2:495223/148‑1 (MQ=255)
                                                                                                          cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactc                                      <  1:548589/148‑1 (MQ=255)
                                                                                                               ggCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCAcccc                                 >  2:2483788/1‑148 (MQ=255)
                                                                                                                gCACAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcg                                <  2:3481392/148‑1 (MQ=255)
                                                                                                                   aaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGt                             <  2:2158874/148‑1 (MQ=255)
                                                                                                                   aaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGt                             <  1:2045258/148‑1 (MQ=255)
                                                                                                                   aaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGt                             <  2:2746321/148‑1 (MQ=255)
                                                                                                                   aaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGt                             <  2:363866/148‑1 (MQ=255)
                                                                                                                    aaGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                            >  2:2418142/1‑148 (MQ=255)
                                                                                                                     aGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                           >  2:945763/1‑148 (MQ=255)
                                                                                                                         gCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                       <  2:2305608/148‑1 (MQ=255)
                                                                                                                          cGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcgc                      <  1:1084621/148‑1 (MQ=255)
                                                                                                                           gTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCt                     <  1:2573836/148‑1 (MQ=255)
                                                                                                                              aGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGa                  >  2:2587503/1‑148 (MQ=255)
                                                                                                                                cTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGt                >  2:2059593/1‑148 (MQ=255)
                                                                                                                                  ggAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTgg              <  2:2173858/148‑1 (MQ=255)
                                                                                                                                   gAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa             <  2:1319303/148‑1 (MQ=255)
                                                                                                                                   gAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa             <  2:1318739/148‑1 (MQ=255)
                                                                                                                                      aGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt          <  2:2774014/148‑1 (MQ=255)
                                                                                                                                      aGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt          <  1:2562582/148‑1 (MQ=255)
                                                                                                                                        cATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAACCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTAt        <  2:1772704/148‑1 (MQ=255)
                                                                                                                                           ttCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGc     <  1:2483970/148‑1 (MQ=255)
                                                                                                                                           ttCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAGTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGc     <  1:3557404/148‑1 (MQ=255)
                                                                                                                                              gCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTg  >  2:729968/1‑148 (MQ=255)
                                                                                                                                              gCCTAACGCCACCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTg  >  2:3047652/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTG  >  CP009273/2852709‑2852998

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