Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,852,854 G→T 90.5% 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,8540GT90.5% 896.5 / 87.2 327P801P (CCG→CCTmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (12/19);  new base T (145/151);  total (157/170)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 327 total aligned reads displayed.

GGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAA  >  CP009273/2852708‑2853000
                                                                                                                                                  |                                                                                                                                                  
ggCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTaa                                                                                                                                                  >  2:3068991/1‑149 (MQ=255)
 gCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAc                                                                                                                                                 >  1:3156074/1‑149 (MQ=255)
  cGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcg                                                                                                                                                >  2:2249532/1‑149 (MQ=255)
   gACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgc                                                                                                                                               <  2:2092635/149‑1 (MQ=255)
    acacCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgcc                                                                                                                                              >  2:448726/1‑149 (MQ=255)
    acacCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgcc                                                                                                                                              >  2:2816820/1‑149 (MQ=255)
         aTTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAAcgccgccgc                                                                                                                                         >  2:4189843/1‑149 (MQ=255)
          ttGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCt                                                                                                                                        >  1:3628071/1‑149 (MQ=255)
            gCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGTCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTAc                                                                                                                                      <  2:2398011/149‑1 (MQ=255)
             ccTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACg                                                                                                                                     >  2:680898/1‑149 (MQ=255)
              cTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGc                                                                                                                                    >  1:3829387/1‑149 (MQ=255)
              cTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGc                                                                                                                                    >  1:1336546/1‑149 (MQ=255)
               tttATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCa                                                                                                                                   >  2:4416017/1‑149 (MQ=255)
               tttATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCa                                                                                                                                   <  1:2637087/149‑1 (MQ=255)
                 tATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAg                                                                                                                                 <  1:2957900/149‑1 (MQ=255)
                     cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGa                                                                                                                             >  2:1650076/1‑149 (MQ=255)
                      acaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGAt                                                                                                                            >  1:551245/1‑149 (MQ=255)
                           gTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTAc                                                                                                                       <  2:4497236/149‑1 (MQ=255)
                                  aTGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATg                                                                                                                <  2:1984604/149‑1 (MQ=255)
                                   tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGt                                                                                                               >  1:1215399/1‑149 (MQ=255)
                                      cgcgGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCtt                                                                                                            <  2:4001825/149‑1 (MQ=255)
                                         ggcgAGCAAAATCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGc                                                                                                         <  1:338577/149‑1 (MQ=255)
                                         ggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGc                                                                                                         <  2:978002/149‑1 (MQ=255)
                                          gcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCt                                                                                                        >  1:4416767/1‑149 (MQ=255)
                                           cgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTg                                                                                                       <  1:231030/149‑1 (MQ=255)
                                           cgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTg                                                                                                       <  1:4582680/149‑1 (MQ=255)
                                              gCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAACTGTCTTTGCTGTCa                                                                                                    <  1:3926037/149‑1 (MQ=255)
                                              gCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCa                                                                                                    <  1:1492908/149‑1 (MQ=255)
                                               cAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGAGGATGGTACGCAAATGTCTTTGCTGTCAg                                                                                                   >  1:4092302/1‑149 (MQ=255)
                                                  aaGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTAc                                                                                                >  2:4210161/1‑149 (MQ=255)
                                                    gCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCa                                                                                              <  1:157833/149‑1 (MQ=255)
                                                    gCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCa                                                                                              >  1:4755633/1‑149 (MQ=255)
                                                      tACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCaga                                                                                            <  1:1816815/149‑1 (MQ=255)
                                                       aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaa                                                                                           <  2:1500305/149‑1 (MQ=255)
                                                        cGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaag                                                                                          <  2:4379634/149‑1 (MQ=255)
                                                         ggCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaaga                                                                                         <  2:3812789/149‑1 (MQ=255)
                                                         ggCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCagaaga                                                                                         >  2:2717036/1‑149 (MQ=255)
                                                          gCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAgaagaa                                                                                        >  1:3845508/1‑149 (MQ=255)
                                                          gCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAgaagaa                                                                                        <  2:4049367/149‑1 (MQ=255)
                                                           ccTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATTGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGaaa                                                                                       <  2:2897589/149‑1 (MQ=255)
                                                            cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAAc                                                                                      <  2:1813614/149‑1 (MQ=255)
                                                            cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAAc                                                                                      <  1:4656516/149‑1 (MQ=255)
                                                            cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAAc                                                                                      >  2:1243602/1‑149 (MQ=255)
                                                             tGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACt                                                                                     >  2:4763334/1‑149 (MQ=255)
                                                               gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTc                                                                                   >  1:1989746/1‑149 (MQ=255)
                                                               gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTc                                                                                   >  1:3075065/1‑149 (MQ=255)
                                                                 ggTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGc                                                                                 >  1:3698589/1‑149 (MQ=255)
                                                                  gTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGcc                                                                                >  1:2840081/1‑149 (MQ=255)
                                                                       aGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCgg                                                                           >  1:4361563/1‑149 (MQ=255)
                                                                        gCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                          >  2:3339661/1‑149 (MQ=255)
                                                                              gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGACAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCCAATGTCTTTGCTGTCAGTACCAGAAGAACCTTCGCCTGCGGTCGAAGc                                                                    <  2:2969107/149‑1 (MQ=255)
                                                                              gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGc                                                                    <  2:2161021/149‑1 (MQ=255)
                                                                                aGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctc                                                                  >  1:4494778/1‑149 (MQ=255)
                                                                                 ggCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctct                                                                 <  2:2036941/149‑1 (MQ=255)
                                                                                 ggCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctct                                                                 >  1:286744/1‑149 (MQ=255)
                                                                                 ggCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGctct                                                                 <  1:1039282/149‑1 (MQ=255)
                                                                                     tGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaa                                                             >  1:1644907/1‑149 (MQ=255)
                                                                                         aaaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATc                                                         <  2:1780457/149‑1 (MQ=255)
                                                                                          aaaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                        <  1:2773156/149‑1 (MQ=255)
                                                                                           aaGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCtt                                                       >  1:1465800/1‑149 (MQ=255)
                                                                                            agagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTg                                                      >  1:974526/1‑149 (MQ=255)
                                                                                            agagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTg                                                      >  2:3243006/1‑149 (MQ=255)
                                                                                             gagGTTATTCAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGa                                                     <  1:728181/149‑1 (MQ=255)
                                                                                               ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                   <  2:4269116/149‑1 (MQ=255)
                                                                                               ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                   >  2:1205290/1‑149 (MQ=255)
                                                                                               ggTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                   >  1:4694366/1‑149 (MQ=255)
                                                                                                 ttattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                 >  2:285676/1‑149 (MQ=255)
                                                                                                  tattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCgg                                                >  2:3553207/1‑149 (MQ=255)
                                                                                                  tattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCgg                                                >  2:3504151/1‑149 (MQ=255)
                                                                                                   attaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                               >  1:2168101/1‑149 (MQ=255)
                                                                                                   attaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                               <  1:2946975/149‑1 (MQ=255)
                                                                                                      aaGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtc                                            <  2:2337845/149‑1 (MQ=255)
                                                                                                       aGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtca                                           <  2:416016/149‑1 (MQ=255)
                                                                                                           cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactca                                       <  2:3422486/149‑1 (MQ=255)
                                                                                                            gcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcactcac                                      >  1:3711381/1‑149 (MQ=255)
                                                                                                                 gCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgc                                 >  2:189991/1‑149 (MQ=255)
                                                                                                                    aaaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                              <  2:793146/149‑1 (MQ=255)
                                                                                                                     aaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                             >  2:2545183/1‑149 (MQ=255)
                                                                                                                     aaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                             >  2:3208087/1‑149 (MQ=255)
                                                                                                                     aaGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCa                             >  2:1469140/1‑149 (MQ=255)
                                                                                                                      aGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAg                            <  1:2493635/149‑1 (MQ=255)
                                                                                                                      aGCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAg                            >  2:3951709/1‑149 (MQ=255)
                                                                                                                       gCTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAgg                           <  1:1884853/149‑1 (MQ=255)
                                                                                                                        cTGCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgc                          <  2:919348/149‑1 (MQ=255)
                                                                                                                          gCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcgc                        <  2:2440641/149‑1 (MQ=255)
                                                                                                                          gCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcgc                        <  1:2694507/149‑1 (MQ=255)
                                                                                                                          gCGTGAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcgc                        >  2:4057314/1‑149 (MQ=255)
                                                                                                                              gAGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGa                    >  1:3474248/1‑149 (MQ=255)
                                                                                                                               aGCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAg                   <  1:3517412/149‑1 (MQ=255)
                                                                                                                                gCTGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGt                  <  2:3627504/149‑1 (MQ=255)
                                                                                                                                gCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGt                  >  1:1706854/1‑149 (MQ=255)
                                                                                                                                  tGGAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTgg                <  1:3649682/149‑1 (MQ=255)
                                                                                                                                   ggAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa               <  2:1439397/149‑1 (MQ=255)
                                                                                                                                   ggAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa               <  1:81640/149‑1 (MQ=255)
                                                                                                                                   ggAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCGGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGa               <  2:3829387/149‑1 (MQ=255)
                                                                                                                                    gAAAGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt              >  2:460245/1‑149 (MQ=255)
                                                                                                                                      aaGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt            <  2:188443/149‑1 (MQ=255)
                                                                                                                                       aGCATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTa           >  2:1354164/1‑149 (MQ=255)
                                                                                                                                       aGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTa           <  2:531184/149‑1 (MQ=255)
                                                                                                                                         cATTTCGCCTAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATc         <  2:372745/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGCGACACCATTGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAA  >  CP009273/2852708‑2853000

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: