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.0% 621.5 / 8.7 201P801P (CCG→CCTmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (101/94);  total (104/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 201 total aligned reads displayed.

GCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCT  >  CP009273/2852709‑2852996
                                                                                                                                                 |                                                                                                                                              
gCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTaa                                                                                                                                              <  2:1285549/148‑1 (MQ=255)
 cGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAAc                                                                                                                                             >  1:1532213/1‑148 (MQ=255)
   acacCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgc                                                                                                                                           >  2:1687699/1‑148 (MQ=255)
     acCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccg                                                                                                                                         <  1:2785703/148‑1 (MQ=255)
      ccATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccgc                                                                                                                                        >  1:1052020/1‑148 (MQ=255)
      ccATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccgc                                                                                                                                        >  2:2594464/1‑148 (MQ=255)
        aTTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccgccg                                                                                                                                      >  1:787984/1‑148 (MQ=255)
           gCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTa                                                                                                                                   <  1:2072517/148‑1 (MQ=255)
             cTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACg                                                                                                                                 <  1:2292836/148‑1 (MQ=255)
                tATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCaa                                                                                                                              <  2:1793770/148‑1 (MQ=255)
                tATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCaa                                                                                                                              <  1:170897/148‑1 (MQ=255)
                    cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTgg                                                                                                                          <  1:600682/148‑1 (MQ=255)
                    cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTgg                                                                                                                          <  2:1174832/148‑1 (MQ=255)
                    cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTgg                                                                                                                          >  2:3318497/1‑148 (MQ=255)
                     acaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGa                                                                                                                         <  1:2121963/148‑1 (MQ=255)
                     acaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGa                                                                                                                         <  1:3151501/148‑1 (MQ=255)
                      caGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCACGTGGAt                                                                                                                        >  1:2847296/1‑148 (MQ=255)
                      caGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGAt                                                                                                                        <  2:7886/148‑1 (MQ=255)
                        gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATgg                                                                                                                      <  2:3093511/148‑1 (MQ=255)
                          gTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTa                                                                                                                    <  2:2260772/148‑1 (MQ=255)
                           tGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTAc                                                                                                                   >  2:2628731/1‑148 (MQ=255)
                            gCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACg                                                                                                                  <  1:1573324/148‑1 (MQ=255)
                               ggATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCaa                                                                                                               >  2:1433748/1‑148 (MQ=255)
                                gATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCaaa                                                                                                              <  2:922520/148‑1 (MQ=255)
                                    gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTc                                                                                                          >  2:1835079/1‑148 (MQ=255)
                                     cgcgGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCt                                                                                                         <  2:1980279/148‑1 (MQ=255)
                                     cgcgGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCt                                                                                                         <  1:1956694/148‑1 (MQ=255)
                                       cggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCttt                                                                                                       <  2:81369/148‑1 (MQ=255)
                                            aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                                                                                  >  2:132671/1‑148 (MQ=255)
                                            aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                                                                                  <  2:850278/148‑1 (MQ=255)
                                            aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                                                                                  <  2:850312/148‑1 (MQ=255)
                                             gCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTc                                                                                                 <  2:1508655/148‑1 (MQ=255)
                                              cAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCa                                                                                                >  2:2585369/1‑148 (MQ=255)
                                               aaaaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAg                                                                                               >  2:3451687/1‑148 (MQ=255)
                                                aaaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGt                                                                                              >  1:1246940/1‑148 (MQ=255)
                                                 aaGCTACGGCCTGGCGGTTGCAGCTCAGGAAGGCGTGCAAAAAGAGCTTATTAAGCGCGCACGGCAAAAGATGCGTGAGCTGGAAAGCATATAGCCTAACGCAGCCGCTACGCAACTGGATGGTACGCAAATGTCTATGCAGTCAGTa                                                                                             >  2:1139434/1‑148 (MQ=255)
                                                      aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCaga                                                                                        >  2:2688326/1‑148 (MQ=255)
                                                       cGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaa                                                                                       >  2:2931969/1‑148 (MQ=255)
                                                          ccTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAgaagaa                                                                                    >  1:2501481/1‑148 (MQ=255)
                                                          ccTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAgaagaa                                                                                    <  1:2386431/148‑1 (MQ=255)
                                                           cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGaaa                                                                                   <  1:1369442/148‑1 (MQ=255)
                                                            tGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAAc                                                                                  <  2:931807/148‑1 (MQ=255)
                                                            tGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAAc                                                                                  >  1:3210944/1‑148 (MQ=255)
                                                              gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACtt                                                                                <  2:2020332/148‑1 (MQ=255)
                                                                 gTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGc                                                                             <  2:654092/148‑1 (MQ=255)
                                                                 gTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGc                                                                             <  2:651648/148‑1 (MQ=255)
                                                                   tGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCt                                                                           >  1:3147095/1‑148 (MQ=255)
                                                                     cAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGc                                                                         >  2:1320063/1‑148 (MQ=255)
                                                                        ctctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                      >  2:2817401/1‑148 (MQ=255)
                                                                        ctctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                      <  1:668345/148‑1 (MQ=255)
                                                                        ctctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                      >  1:2458152/1‑148 (MQ=255)
                                                                             gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAg                                                                 >  2:1909715/1‑148 (MQ=255)
                                                                                ggCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctc                                                              >  1:3093678/1‑148 (MQ=255)
                                                                                   gTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTgg                                                           <  1:2587016/148‑1 (MQ=255)
                                                                                        aaaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAAt                                                      <  2:2043887/148‑1 (MQ=255)
                                                                                        aaaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAAt                                                      <  2:914232/148‑1 (MQ=255)
                                                                                         aaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATc                                                     <  1:566250/148‑1 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGTAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                    >  2:1702407/1‑148 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                    <  1:109770/148‑1 (MQ=255)
                                                                                            gagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTg                                                  <  2:278469/148‑1 (MQ=255)
                                                                                             agGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGa                                                 <  1:2229253/148‑1 (MQ=255)
                                                                                             agGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGa                                                 >  1:2907086/1‑148 (MQ=255)
                                                                                             agGTTATTAAGCGCGCACGACACAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGa                                                 <  2:2230119/148‑1 (MQ=255)
                                                                                              ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                >  2:2794566/1‑148 (MQ=255)
                                                                                                  attaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCgg                                            <  1:1867259/148‑1 (MQ=255)
                                                                                                     aaGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAtt                                         <  2:3362544/148‑1 (MQ=255)
                                                                                                      aGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtc                                        <  2:2360990/148‑1 (MQ=255)
                                                                                                      aGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtc                                        <  1:2001797/148‑1 (MQ=255)
                                                                                                         gcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcact                                     <  2:1710289/148‑1 (MQ=255)
                                                                                                          cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactc                                    <  1:1301614/148‑1 (MQ=255)
                                                                                                            ctcGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactcac                                  <  2:3558025/146‑1 (MQ=255)
                                                                                                                   aaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGt                           >  1:739658/1‑148 (MQ=255)
                                                                                                                   aaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGt                           <  1:2533769/148‑1 (MQ=255)
                                                                                                                    aaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                          >  2:1801255/1‑148 (MQ=255)
                                                                                                                    aaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                          >  2:1687301/1‑148 (MQ=255)
                                                                                                                     aGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                         <  1:2594464/148‑1 (MQ=255)
                                                                                                                     aGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                         >  1:335251/1‑148 (MQ=255)
                                                                                                                     aGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                         >  1:334901/1‑148 (MQ=255)
                                                                                                                       cTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAgg                       >  1:2439460/1‑148 (MQ=255)
                                                                                                                           gTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCt                   <  2:2458152/148‑1 (MQ=255)
                                                                                                                            tGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTg                  >  1:3007269/1‑148 (MQ=255)
                                                                                                                            tGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTg                  >  1:3007314/1‑148 (MQ=255)
                                                                                                                             gAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTgg                 >  1:3160284/1‑148 (MQ=255)
                                                                                                                               gCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAg               >  2:1638022/1‑148 (MQ=255)
                                                                                                                                cTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGt              >  1:525792/1‑148 (MQ=255)
                                                                                                                                 tGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTg             >  2:2098663/1‑148 (MQ=255)
                                                                                                                                  ggAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTgg            <  2:3508529/148‑1 (MQ=255)
                                                                                                                                  ggAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTgg            >  1:2431592/1‑148 (MQ=255)
                                                                                                                                   gAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa           >  1:2190781/1‑148 (MQ=255)
                                                                                                                                      aGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt        <  1:147672/148‑1 (MQ=255)
                                                                                                                                       gCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTa       <  1:2063049/148‑1 (MQ=255)
                                                                                                                                       gCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTa       >  2:754205/1‑148 (MQ=255)
                                                                                                                                       gCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTa       >  2:813827/1‑148 (MQ=255)
                                                                                                                                          tttCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCg    >  2:3523126/1‑148 (MQ=255)
                                                                                                                                          tttCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCg    >  2:1764902/1‑148 (MQ=255)
                                                                                                                                          tttCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCg    <  1:1632750/148‑1 (MQ=255)
                                                                                                                                          tttCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCg    <  1:444346/148‑1 (MQ=255)
                                                                                                                                            tCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCt  <  1:2655152/148‑1 (MQ=255)
                                                                                                                                            tCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCt  <  2:104568/148‑1 (MQ=255)
                                                                                                                                            tCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCt  <  2:104418/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCT  >  CP009273/2852709‑2852996

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: