Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,852,854 G→T 87.3% 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,8540GT87.3% 461.7 / 71.6 197P801P (CCG→CCTmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base T (83/89);  total (94/103)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 197 total aligned reads displayed.

GCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGA  >  CP009273/2852709‑2852999
                                                                                                                                                 |                                                                                                                                                 
gCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTaa                                                                                                                                                 <  2:3291741/148‑1 (MQ=255)
gCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTaa                                                                                                                                                 <  2:1567797/148‑1 (MQ=255)
gCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTaa                                                                                                                                                 <  1:1922994/148‑1 (MQ=255)
    cacCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgcc                                                                                                                                             >  1:184736/1‑148 (MQ=255)
    cacCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgcc                                                                                                                                             >  1:184605/1‑148 (MQ=255)
    cacCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAAcgcc                                                                                                                                             >  2:2584417/1‑148 (MQ=255)
     acCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGTCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccg                                                                                                                                            >  1:852447/1‑148 (MQ=255)
     acCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGTCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccg                                                                                                                                            >  1:2542152/1‑148 (MQ=255)
     acCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTCAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccg                                                                                                                                            <  1:206527/148‑1 (MQ=255)
          tGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCt                                                                                                                                       >  2:2045559/1‑148 (MQ=255)
          tGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCt                                                                                                                                       >  2:1156228/1‑148 (MQ=255)
            ccTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTAc                                                                                                                                     >  1:97092/1‑148 (MQ=255)
             cTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACg                                                                                                                                    <  2:685734/148‑1 (MQ=255)
              tttATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGc                                                                                                                                   >  2:1909456/1‑148 (MQ=255)
              tttATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGc                                                                                                                                   >  2:1411811/1‑148 (MQ=255)
                      caGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGAt                                                                                                                           >  2:1988627/1‑148 (MQ=255)
                         cGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGt                                                                                                                        >  1:1642514/1‑148 (MQ=255)
                          gTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTa                                                                                                                       <  1:2925830/148‑1 (MQ=255)
                            gCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACg                                                                                                                     <  1:1214236/148‑1 (MQ=255)
                               ggATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCaa                                                                                                                  >  1:2147358/1‑148 (MQ=255)
                                gATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCaaa                                                                                                                 <  1:2587929/148‑1 (MQ=255)
                                gATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCaaa                                                                                                                 >  2:1345278/1‑148 (MQ=255)
                                  tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATg                                                                                                               <  2:616770/148‑1 (MQ=255)
                                   ggCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGt                                                                                                              <  2:287051/148‑1 (MQ=255)
                                         gcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGc                                                                                                        <  1:1347242/148‑1 (MQ=255)
                                          cgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCACATGTCTTTGCt                                                                                                       >  1:26888/1‑148 (MQ=255)
                                          cgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCt                                                                                                       <  1:2884590/148‑1 (MQ=255)
                                          cgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCt                                                                                                       >  1:921225/1‑148 (MQ=255)
                                            aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTtgttgt                                                                                                     <  1:1911907/148‑1 (MQ=255)
                                            aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                                                                                     <  2:1235352/148‑1 (MQ=255)
                                            aGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGt                                                                                                     <  1:1443242/148‑1 (MQ=255)
                                              cAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCa                                                                                                   >  1:3130270/1‑148 (MQ=255)
                                               aaaaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAg                                                                                                  >  1:2952838/1‑148 (MQ=255)
                                                aaaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGt                                                                                                 >  2:43040/1‑148 (MQ=255)
                                                    cTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCa                                                                                             >  1:1450318/1‑148 (MQ=255)
                                                      aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCaga                                                                                           <  2:2184436/148‑1 (MQ=255)
                                                      aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCaga                                                                                           <  2:520718/148‑1 (MQ=255)
                                                           cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAACGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGaaa                                                                                      <  2:841308/148‑1 (MQ=255)
                                                           cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGaaa                                                                                      <  1:2780171/148‑1 (MQ=255)
                                                           cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGaaa                                                                                      <  1:1985634/148‑1 (MQ=255)
                                                            tGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAAc                                                                                     <  1:3310761/148‑1 (MQ=255)
                                                              gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACtt                                                                                   <  1:1831151/148‑1 (MQ=255)
                                                                ggTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCg                                                                                 <  1:1164074/148‑1 (MQ=255)
                                                                 gTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGc                                                                                <  1:534916/148‑1 (MQ=255)
                                                                  ttGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGcc                                                                               <  2:3243442/148‑1 (MQ=255)
                                                                     cAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGc                                                                            >  1:124015/1‑148 (MQ=255)
                                                                     cAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGc                                                                            >  2:27628/1‑148 (MQ=255)
                                                                        ctctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                         <  2:2645918/148‑1 (MQ=255)
                                                                        ctctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                         <  1:2051529/148‑1 (MQ=255)
                                                                        ctctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                         <  1:2472926/148‑1 (MQ=255)
                                                                          ctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                       >  1:2979809/1‑148 (MQ=255)
                                                                          ctGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                       >  1:2916228/1‑148 (MQ=255)
                                                                            ggCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGaa                                                                     >  2:1391377/1‑148 (MQ=255)
                                                                             gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAg                                                                    >  2:2972549/1‑148 (MQ=255)
                                                                               aGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                  <  2:167131/148‑1 (MQ=255)
                                                                               aGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                  <  2:167172/148‑1 (MQ=255)
                                                                                 gCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctct                                                                >  1:1320080/1‑148 (MQ=255)
                                                                                    tGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTCCGTGAGCTGGAAAGCATATCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTAGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGa                                                             >  2:1726565/1‑148 (MQ=255)
                                                                                        aaaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAAt                                                         <  1:354955/148‑1 (MQ=255)
                                                                                        aaaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAAt                                                         >  1:1703043/1‑148 (MQ=255)
                                                                                          aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACTCATGTGGATGGTAAGCAAATGTATTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                       <  2:2910570/148‑1 (MQ=255)
                                                                                            gagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTg                                                     <  1:729415/148‑1 (MQ=255)
                                                                                             agGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGa                                                    <  1:3231540/148‑1 (MQ=255)
                                                                                                ttattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATcc                                                 >  2:100584/1‑148 (MQ=255)
                                                                                                 tattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                <  2:1711721/148‑1 (MQ=255)
                                                                                                 gaagcaGCCCGCACGCCACAAGCTCCGTGACCTGGCAAGCCGTCCGCCTAACGCCGCCGCTACGCACGTGGATGGTCCGCAAATGTGTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                <  2:1888510/143‑1 (MQ=255)
                                                                                                  attaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCgg                                               <  2:2372856/148‑1 (MQ=255)
                                                                                                   ttaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACTCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                              <  2:2829504/148‑1 (MQ=255)
                                                                                                     aaGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAtt                                            <  1:1447813/148‑1 (MQ=255)
                                                                                                         gcgcACGGCATAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGTTGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcact                                        <  2:1816759/148‑1 (MQ=255)
                                                                                                         gcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcact                                        <  2:97092/148‑1 (MQ=255)
                                                                                                           gcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactca                                      <  1:1876446/148‑1 (MQ=255)
                                                                                                           gcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactca                                      <  1:2379483/148‑1 (MQ=255)
                                                                                                              cGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCAccc                                   >  2:780705/1‑148 (MQ=255)
                                                                                                               ggCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCAcccc                                  >  1:1092706/1‑148 (MQ=255)
                                                                                                                  aaaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgcg                               >  1:457630/1‑148 (MQ=255)
                                                                                                                    aaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                             >  2:790683/1‑148 (MQ=255)
                                                                                                                     aGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                            <  1:2772415/148‑1 (MQ=255)
                                                                                                                     aGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                            <  1:1036165/148‑1 (MQ=255)
                                                                                                                         gCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                        <  2:1360826/148‑1 (MQ=255)
                                                                                                                         gCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                        <  2:2500165/148‑1 (MQ=255)
                                                                                                                         gCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                        <  1:1909456/148‑1 (MQ=255)
                                                                                                                         gCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                        <  1:1135715/148‑1 (MQ=255)
                                                                                                                          cGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcgc                       >  2:1639086/1‑148 (MQ=255)
                                                                                                                          cGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcgc                       >  1:672387/1‑148 (MQ=255)
                                                                                                                           gTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCt                      <  2:825872/148‑1 (MQ=255)
                                                                                                                             gAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTgg                    >  1:584034/1‑148 (MQ=255)
                                                                                                                             gAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTgg                    >  1:2247920/1‑148 (MQ=255)
                                                                                                                                   gAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa              >  1:259223/1‑148 (MQ=255)
                                                                                                                                   gAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa              >  2:3149647/1‑148 (MQ=255)
                                                                                                                                    aaaGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt             <  1:139317/148‑1 (MQ=255)
                                                                                                                                     aaGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAtt            >  2:252627/1‑148 (MQ=255)
                                                                                                                                      aGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt           <  1:599889/148‑1 (MQ=255)
                                                                                                                                      aGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt           <  1:795800/148‑1 (MQ=255)
                                                                                                                                      aGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCAACCCGCGTCAGGCGCTGGAGTGGAttt           >  2:1304487/1‑148 (MQ=255)
                                                                                                                                        cATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTAt         <  2:873171/148‑1 (MQ=255)
                                                                                                                                          tttCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCg       <  1:2584417/148‑1 (MQ=255)
                                                                                                                                           ttCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAATAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGc      <  2:1442238/148‑1 (MQ=255)
                                                                                                                                           gtCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGc      <  2:1819792/147‑1 (MQ=255)
                                                                                                                                              gCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTg   >  1:2189067/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGA  >  CP009273/2852709‑2852999

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: