Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,850,790 C→A 100% I113I (ATC→ATA)  mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,850,7900CA96.0% 883.0 / 25.4 275I113I (ATC→ATA) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (8/3);  new base A (127/137);  total (135/140)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 277 total aligned reads displayed.

GCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGA  >  CP009273/2850643‑2850936
                                                                                                                                                   |                                                                                                                                                  
gCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGt                                                                                                                                                  >  2:1734876/1‑150 (MQ=255)
  cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                                >  1:2922117/1‑150 (MQ=255)
  cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                                >  2:4395426/1‑150 (MQ=255)
      cccGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGcc                                                                                                                                            <  1:3655538/150‑1 (MQ=255)
       ccGATGGCGGGGATTCACTACCATGCGGTGGAAAACTATCTCGCCAAACTGTTGAATCAGTTAGAGTCCGTTGCCATCTTCGAACAAATTGTCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCa                                                                                                                                           <  1:1407299/150‑1 (MQ=255)
         gATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAgg                                                                                                                                         <  2:3697139/150‑1 (MQ=255)
          aTGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGc                                                                                                                                        >  1:272741/1‑150 (MQ=255)
           tGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCa                                                                                                                                       <  2:4477953/150‑1 (MQ=255)
            ggCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCAc                                                                                                                                      <  1:2527189/150‑1 (MQ=255)
            ggCGGGGATTCCCGACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCAc                                                                                                                                      <  2:1469128/150‑1 (MQ=255)
              cGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCa                                                                                                                                    >  1:1209460/1‑150 (MQ=255)
                gggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATc                                                                                                                                  >  2:2223526/1‑150 (MQ=255)
                  gATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAg                                                                                                                                <  2:2395283/150‑1 (MQ=255)
                          aCCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAg                                                                                                                        <  1:4516293/150‑1 (MQ=255)
                           ccATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGc                                                                                                                       >  1:2870694/1‑150 (MQ=255)
                            cATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGcc                                                                                                                      >  1:135123/1‑150 (MQ=255)
                             aTGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGccc                                                                                                                     <  1:608541/150‑1 (MQ=255)
                             aTGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGccc                                                                                                                     <  1:3228306/150‑1 (MQ=255)
                             aTGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGccc                                                                                                                     <  2:607590/150‑1 (MQ=255)
                                 ggTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCAAGGCACCATCAGCGATGAAGCCCTGtt                                                                                                                 <  1:4786/150‑1 (MQ=255)
                                  gTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTg                                                                                                                >  2:3074616/1‑150 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGGGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                              >  1:4248169/1‑150 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                              <  1:1988246/150‑1 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                              <  1:3586454/150‑1 (MQ=255)
                                      aaaaCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAgg                                                                                                            >  1:2612809/1‑150 (MQ=255)
                                         aCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGc                                                                                                         <  1:3119511/150‑1 (MQ=255)
                                         aCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGc                                                                                                         <  1:1180501/150‑1 (MQ=255)
                                          cTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCg                                                                                                        <  2:183523/150‑1 (MQ=255)
                                            aTCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTc                                                                                                      <  1:3961386/150‑1 (MQ=255)
                                            aTCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTc                                                                                                      <  1:3961275/150‑1 (MQ=255)
                                                    aaaCTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAAc                                                                                              <  1:189528/150‑1 (MQ=255)
                                                       cTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctg                                                                                           <  1:2417667/150‑1 (MQ=255)
                                                       cTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctg                                                                                           <  1:2417301/150‑1 (MQ=255)
                                                       cTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGCCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctg                                                                                           <  2:2522443/150‑1 (MQ=255)
                                                        tGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctgc                                                                                          >  2:4301820/1‑150 (MQ=255)
                                                           tGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTgg                                                                                       <  1:1439465/150‑1 (MQ=255)
                                                             aaTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcg                                                                                     >  1:4133957/1‑150 (MQ=255)
                                                              aTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcgg                                                                                    >  1:3771255/1‑150 (MQ=255)
                                                               tCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcggc                                                                                   >  2:3391679/1‑150 (MQ=255)
                                                                   ggAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATc                                                                               <  2:792969/150‑1 (MQ=255)
                                                                    gagaGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCt                                                                              <  2:4167547/150‑1 (MQ=255)
                                                                     agagTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTg                                                                             <  1:665625/150‑1 (MQ=255)
                                                                     agagTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTg                                                                             >  1:1872567/1‑150 (MQ=255)
                                                                       agTCCGTTGCCATCTACGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGc                                                                           >  2:3437029/1‑150 (MQ=255)
                                                                        gTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCa                                                                          >  2:1457962/1‑150 (MQ=255)
                                                                           cGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGa                                                                       <  2:2798853/150‑1 (MQ=255)
                                                                            gTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGAc                                                                      <  1:2721230/150‑1 (MQ=255)
                                                                             ttGCCATCTGCGAACAAATTGGCGCTCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGAcc                                                                     <  1:2976633/150‑2 (MQ=255)
                                                                              tGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAg                                                                    <  1:3460694/150‑1 (MQ=255)
                                                                              tGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAg                                                                    >  1:3561701/1‑150 (MQ=255)
                                                                              tGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAg                                                                    >  1:2027511/1‑150 (MQ=255)
                                                                               gCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGc                                                                   <  2:4418983/150‑1 (MQ=255)
                                                                                ccATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCa                                                                  >  2:2193189/1‑150 (MQ=255)
                                                                                 cATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTAGAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCaa                                                                 >  2:3789693/1‑150 (MQ=255)
                                                                                   tCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAg                                                               <  1:3079511/150‑1 (MQ=255)
                                                                                   tCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAg                                                               >  1:3696412/1‑150 (MQ=255)
                                                                                    cTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAgg                                                              >  2:3967153/1‑150 (MQ=255)
                                                                                     tGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGt                                                             <  1:6027/150‑1 (MQ=255)
                                                                                       cGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTGCGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGttt                                                           <  1:390625/150‑1 (MQ=255)
                                                                                        gAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTc                                                          >  1:239613/1‑150 (MQ=255)
                                                                                         aaCAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCg                                                         <  2:3328327/150‑1 (MQ=255)
                                                                                          aCAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCgg                                                        >  1:387458/1‑150 (MQ=255)
                                                                                          aCAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCgg                                                        <  2:4605426/150‑1 (MQ=255)
                                                                                           cAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGc                                                       >  1:2466121/1‑150 (MQ=255)
                                                                                            aaaTTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCt                                                      <  1:3301111/150‑1 (MQ=255)
                                                                                              aTTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAc                                                    <  2:1325996/150‑1 (MQ=255)
                                                                                              aTTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAc                                                    <  2:3906735/150‑1 (MQ=255)
                                                                                               ttGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAcg                                                   >  2:3605683/1‑150 (MQ=255)
                                                                                                tGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAcgc                                                  >  2:1572773/1‑150 (MQ=255)
                                                                                                 ggCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAcgcg                                                 <  2:3916304/150‑1 (MQ=255)
                                                                                                   cGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGAc                                               >  2:59783/1‑150 (MQ=255)
                                                                                                     aTCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGc                                             >  1:3564707/1‑150 (MQ=255)
                                                                                                     aTCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGc                                             <  1:3470567/150‑1 (MQ=255)
                                                                                                     aTCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGc                                             >  1:4063265/1‑150 (MQ=255)
                                                                                                       ccGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTg                                           <  1:2362055/150‑1 (MQ=255)
                                                                                                       ccGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTg                                           <  1:1460267/150‑1 (MQ=255)
                                                                                                        cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                          <  1:202330/150‑1 (MQ=255)
                                                                                                        cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                          >  1:1013092/1‑150 (MQ=255)
                                                                                                          gCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGat                                        <  2:1994434/150‑1 (MQ=255)
                                                                                                          gCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGat                                        <  2:1262444/150‑1 (MQ=255)
                                                                                                           cGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGata                                       >  2:1076174/1‑150 (MQ=255)
                                                                                                           cGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGata                                       <  1:2660771/150‑1 (MQ=255)
                                                                                                           cGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGata                                       >  2:2172823/1‑150 (MQ=255)
                                                                                                           cGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGata                                       <  2:3675333/150‑1 (MQ=255)
                                                                                                            gACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGatat                                      <  1:2003998/150‑1 (MQ=255)
                                                                                                              ccAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCa                                    >  2:4334067/1‑150 (MQ=255)
                                                                                                               cAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAg                                   >  1:4452481/1‑150 (MQ=255)
                                                                                                                        tCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCg                          >  2:1442738/1‑150 (MQ=255)
                                                                                                                           ggTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGttt                       >  2:1505964/1‑150 (MQ=255)
                                                                                                                            gTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGtttt                      <  2:3269167/150‑1 (MQ=255)
                                                                                                                            gTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGtttt                      >  2:1180941/1‑150 (MQ=255)
                                                                                                                              tGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCg                    <  2:2810661/150‑1 (MQ=255)
                                                                                                                                  cgcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTg                <  1:2155511/150‑1 (MQ=255)
                                                                                                                                  cgcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTg                >  1:2416076/1‑150 (MQ=255)
                                                                                                                                   gcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGa               >  2:4424066/1‑150 (MQ=255)
                                                                                                                                   gcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCAGGATATCAGTTCCGGGCGATTTCGCCTGa               >  1:155144/1‑150 (MQ=255)
                                                                                                                                        gTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGaa          <  2:2894836/150‑1 (MQ=255)
                                                                                                                                             gCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGc     <  1:779727/150‑1 (MQ=255)
                                                                                                                                             gCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGc     <  2:1013133/150‑1 (MQ=255)
                                                                                                                                                tataGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGa  >  2:4636554/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGA  >  CP009273/2850643‑2850936

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: ‑