Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,572,227 G→T 9.8% I213I (ATC→ATA hsdM ← DNA methyltransferase M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,572,2270GT 9.8% 956.4 / 47.0 267I213I (ATC→ATAhsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base G (103/137);  new base T (11/16);  total (114/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

AATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTT  >  CP009273/4572080‑4572375
                                                                                                                                                   |                                                                                                                                                    
aaTATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGa                                                                                                                                                     >  2:4045784/1‑149 (MQ=255)
  tatGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATa                                                                                                                                                   >  1:13685/1‑149 (MQ=255)
   atGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATaa                                                                                                                                                  >  1:3417917/1‑149 (MQ=255)
         cTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAcgcgc                                                                                                                                            <  1:829398/149‑1 (MQ=255)
         cTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAcgcgc                                                                                                                                            <  1:1249142/149‑1 (MQ=255)
         cTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAcgcgc                                                                                                                                            <  2:588443/149‑1 (MQ=255)
            cGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGt                                                                                                                                         <  2:4313112/149‑1 (MQ=255)
            cGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGt                                                                                                                                         <  1:2867665/149‑1 (MQ=255)
            cGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGt                                                                                                                                         >  2:2841734/1‑149 (MQ=255)
             ggCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTg                                                                                                                                        <  1:1083501/149‑1 (MQ=255)
              gCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTgg                                                                                                                                       >  1:3131346/1‑149 (MQ=255)
                aGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGAt                                                                                                                                     <  2:4256258/149‑1 (MQ=255)
                  gTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCt                                                                                                                                   <  2:1472860/149‑1 (MQ=255)
                   ttttCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTg                                                                                                                                  >  1:2074075/1‑149 (MQ=255)
                   ttttCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTg                                                                                                                                  <  2:2041193/149‑1 (MQ=255)
                       cACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGaaa                                                                                                                              <  1:2846794/149‑1 (MQ=255)
                       cACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGaaa                                                                                                                              <  2:4777455/149‑1 (MQ=255)
                        aCCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAAt                                                                                                                             >  2:750180/1‑149 (MQ=255)
                          cGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAACCGCGCGGTGGATCTGGAAATcc                                                                                                                           >  2:556033/1‑149 (MQ=255)
                          cGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATcc                                                                                                                           <  1:2422457/149‑1 (MQ=255)
                           gTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCt                                                                                                                          <  2:279378/149‑1 (MQ=255)
                           gTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCt                                                                                                                          >  1:3660725/1‑149 (MQ=255)
                           gTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGGTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCt                                                                                                                          <  2:4521658/149‑1 (MQ=255)
                                 aCCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgtgt                                                                                                                    <  1:2617624/149‑1 (MQ=255)
                                  cccAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTc                                                                                                                   >  2:1038309/1‑149 (MQ=255)
                                  cccAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTc                                                                                                                   <  1:743273/149‑1 (MQ=255)
                                  cccAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTc                                                                                                                   <  1:3985933/149‑1 (MQ=255)
                                    cAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGc                                                                                                                 <  1:3301612/149‑1 (MQ=255)
                                    cAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAATACAGCGTGTCGcc                                                                                                                >  2:476138/1‑149 (MQ=255)
                                    cAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGc                                                                                                                 >  2:3109090/1‑149 (MQ=255)
                                    cAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGc                                                                                                                 >  1:688132/1‑149 (MQ=255)
                                     agagTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGcc                                                                                                                <  2:4064090/149‑1 (MQ=255)
                                         tgtTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCa                                                                                                            >  2:969327/1‑149 (MQ=255)
                                           ttGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAg                                                                                                          <  1:3444970/149‑1 (MQ=255)
                                             gCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgt                                                                                                        <  2:3720715/149‑1 (MQ=255)
                                                 aGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtcgtc                                                                                                    <  1:4572835/149‑1 (MQ=255)
                                                  gACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTcc                                                                                                   >  2:782949/1‑149 (MQ=255)
                                                   aCGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCa                                                                                                  <  1:4468838/149‑1 (MQ=255)
                                                     ggATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGa                                                                                                <  1:2356000/149‑1 (MQ=255)
                                                     ggATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGa                                                                                                <  1:3296306/149‑1 (MQ=255)
                                                       aTTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGTGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATc                                                                                              <  1:3461221/149‑1 (MQ=255)
                                                       aTTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATc                                                                                              <  1:2175929/149‑1 (MQ=255)
                                                       aTTGCGCCGCCGTGGTCGAGGGTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATc                                                                                              <  2:2524991/149‑1 (MQ=255)
                                                         tGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCAt                                                                                            <  2:3260282/149‑1 (MQ=255)
                                                          gcgcCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCAtt                                                                                           <  1:1754649/149‑1 (MQ=255)
                                                               gccgTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGttt                                                                                      >  1:380676/1‑149 (MQ=255)
                                                                   tGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGa                                                                                  <  1:2138133/149‑1 (MQ=255)
                                                                     gTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACt                                                                                <  2:2395482/149‑1 (MQ=255)
                                                                       cGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTa                                                                              <  2:3210499/149‑1 (MQ=255)
                                                                        gAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTaa                                                                             <  1:3963746/149‑1 (MQ=255)
                                                                         aGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAAc                                                                            >  1:1264865/1‑149 (MQ=255)
                                                                         aGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAAc                                                                            >  1:3710772/1‑149 (MQ=255)
                                                                           gTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACAt                                                                          >  1:4323487/1‑149 (MQ=255)
                                                                              gCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGc                                                                       >  2:2467027/1‑149 (MQ=255)
                                                                               ccTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCg                                                                      <  2:2149481/149‑1 (MQ=255)
                                                                                cTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCATGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCgg                                                                     >  1:866680/1‑149 (MQ=255)
                                                                                 ttCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGt                                                                    <  2:2960956/149‑1 (MQ=255)
                                                                                 tccAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGt                                                                    <  1:4710130/147‑1 (MQ=255)
                                                                                  tCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTc                                                                   <  1:133104/149‑1 (MQ=255)
                                                                                  tCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTc                                                                   <  1:1487965/149‑1 (MQ=255)
                                                                                   cAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCTAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCg                                                                  <  2:4637207/149‑1 (MQ=255)
                                                                                      tatCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCt                                                               >  1:2486241/1‑149 (MQ=255)
                                                                                      tatCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCt                                                               <  2:4363476/149‑1 (MQ=255)
                                                                                       atCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCtt                                                              >  2:3771662/1‑149 (MQ=255)
                                                                                        tCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtc                                                             >  2:1994669/1‑149 (MQ=255)
                                                                                          gTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaa                                                           <  1:3539774/149‑1 (MQ=255)
                                                                                           tGCATCAGGCAGTTCATAAGTGCCAGACTACGGGTGCCGGGCAACATTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTACAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaat                                                          >  2:1939691/1‑149 (MQ=255)
                                                                                            gcagcaGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaatc                                                         <  1:4532093/149‑1 (MQ=255)
                                                                                              agcagGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaatcaa                                                       >  1:2244224/1‑149 (MQ=255)
                                                                                               gcaggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCaaa                                                      >  2:4056878/1‑149 (MQ=255)
                                                                                                 aggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCaaaaa                                                    <  2:928329/149‑1 (MQ=255)
                                                                                                 aggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCaaaaa                                                    <  2:3952087/149‑1 (MQ=255)
                                                                                                 aggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCaaaaa                                                    <  2:3294866/149‑1 (MQ=255)
                                                                                                    cagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGcc                                                 <  1:4093916/149‑1 (MQ=255)
                                                                                                     agTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGccc                                                <  1:3850094/149‑1 (MQ=255)
                                                                                                         catcaGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCccgccg                                            <  2:3479224/149‑1 (MQ=255)
                                                                                                            caGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTAc                                         <  1:2028417/149‑1 (MQ=255)
                                                                                                              gTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCt                                       >  1:4516673/1‑149 (MQ=255)
                                                                                                               tGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTg                                      <  1:2731644/149‑1 (MQ=255)
                                                                                                                  cAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCACGCCGTACCTgccg                                   >  2:4456201/1‑149 (MQ=255)
                                                                                                                   aGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgccgc                                  >  1:1963607/1‑149 (MQ=255)
                                                                                                                    gacgacGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgccgcc                                 >  2:1431273/1‑149 (MQ=255)
                                                                                                                    gacgacGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgccgcc                                 <  2:4187193/149‑1 (MQ=255)
                                                                                                                       gacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCggg                              <  1:1900927/149‑1 (MQ=255)
                                                                                                                       gacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCggg                              >  2:273406/1‑149 (MQ=255)
                                                                                                                        acgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGt                             >  2:361243/1‑149 (MQ=255)
                                                                                                                          gggTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTcc                           <  2:1518457/149‑1 (MQ=255)
                                                                                                                                cGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGcacc                     <  2:3613264/149‑1 (MQ=255)
                                                                                                                                 gggCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAATGCCCGCCGTACCTGCCGCCGGGTCCTGcacca                    >  2:814194/1‑149 (MQ=255)
                                                                                                                                  ggCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGcaccac                   >  1:4306417/1‑149 (MQ=255)
                                                                                                                                   gCACCAGTTCGAGGCCTATAAACGCGCGTTGGATCTGTAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACtt                 <  1:476138/149‑1 (MQ=255)
                                                                                                                                    cACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACtt                 <  2:501338/149‑1 (MQ=255)
                                                                                                                                     aCCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTc                <  1:1114621/149‑1 (MQ=255)
                                                                                                                                      ccAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCa               >  2:2601408/1‑149 (MQ=255)
                                                                                                                                       cAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCAc              >  1:1629710/1‑149 (MQ=255)
                                                                                                                                       cAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCAc              >  1:4248862/1‑149 (MQ=255)
                                                                                                                                       cAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCAc              <  1:3625671/149‑1 (MQ=255)
                                                                                                                                            cGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGc         <  2:1211047/149‑1 (MQ=255)
                                                                                                                                             gAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATAAAAAAGCCCTCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCt        <  1:3161093/149‑1 (MQ=255)
                                                                                                                                                 ccGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCgg    >  1:604614/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTT  >  CP009273/4572080‑4572375

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: