Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,850,790 C→A 73.4% 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,7900CA73.4% 458.1 / 247.2 282I113I (ATC→ATAmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base A (103/104);  total (139/143)
Fisher's exact test for biased strand distribution p-value = 8.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 284 total aligned reads displayed.

GCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCG  >  CP009273/2850643‑2850939
                                                                                                                                                   |                                                                                                                                                     
 ccGAGCCCGATGGCGGGGATTCCCTACCCTGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGtt                                                                                                                                                    <  1:4943271/150‑1 (MQ=255)
  cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTa                                                                                                                                                   <  2:3293918/150‑1 (MQ=255)
  cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                                   <  2:4808306/150‑1 (MQ=255)
   gATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTAc                                                                                                                                                  >  1:1935383/1‑150 (MQ=255)
   gATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTAc                                                                                                                                                  <  1:4577825/150‑1 (MQ=255)
    aTCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACg                                                                                                                                                 <  2:2039949/150‑1 (MQ=255)
         gttGGCGGTGATTCCCTACGATGCGGTGGAAAACTATCTCGCCAAACTGGATAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAgg                                                                                                                                            <  1:2676195/148‑1 (MQ=255)
            ggCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCAc                                                                                                                                         <  1:3328794/150‑1 (MQ=255)
             gCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCAcc                                                                                                                                        >  1:1629944/1‑150 (MQ=255)
              cGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCa                                                                                                                                       >  2:1363461/1‑150 (MQ=255)
              cGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCa                                                                                                                                       <  1:3897860/150‑1 (MQ=255)
                 ggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCa                                                                                                                                    <  2:4206580/150‑1 (MQ=255)
                 ggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCa                                                                                                                                    <  2:360261/150‑1 (MQ=255)
                   aTTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGc                                                                                                                                  >  2:343964/1‑150 (MQ=255)
                     tCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGa                                                                                                                                >  2:313063/1‑150 (MQ=255)
                      cccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGAt                                                                                                                               >  1:2404406/1‑150 (MQ=255)
                      cccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGAt                                                                                                                               >  1:4642447/1‑150 (MQ=255)
                          aCCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAg                                                                                                                           <  2:5043881/150‑1 (MQ=255)
                           ccATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGc                                                                                                                          <  1:5048793/150‑1 (MQ=255)
                            cATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGcc                                                                                                                         >  1:1954217/1‑150 (MQ=255)
                            cATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGcc                                                                                                                         >  2:2597896/1‑150 (MQ=255)
                              tGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCt                                                                                                                       <  2:4010957/150‑1 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                                 >  1:5145932/1‑150 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                                 >  2:1342195/1‑150 (MQ=255)
                                        aaCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAg                                                                                                             >  1:1356878/1‑150 (MQ=255)
                                         aCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGc                                                                                                            >  1:3671685/1‑150 (MQ=255)
                                         aCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGc                                                                                                            <  2:3586375/150‑1 (MQ=255)
                                         aCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGc                                                                                                            <  2:1220297/150‑1 (MQ=255)
                                         aCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGc                                                                                                            <  2:2248596/150‑1 (MQ=255)
                                            aTCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTc                                                                                                         <  2:427631/150‑1 (MQ=255)
                                            aTCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTc                                                                                                         >  1:533184/1‑150 (MQ=255)
                                            aTCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTc                                                                                                         <  1:1556868/150‑1 (MQ=255)
                                             tctcGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCa                                                                                                        >  2:4335916/1‑150 (MQ=255)
                                                cGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGa                                                                                                     <  1:143504/150‑1 (MQ=255)
                                                 gCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGAc                                                                                                    <  2:398971/150‑1 (MQ=255)
                                                 gCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGAc                                                                                                    <  2:1982621/150‑1 (MQ=255)
                                                    aaaCTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAAc                                                                                                 >  2:4386490/1‑150 (MQ=255)
                                                    aaaCTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAAc                                                                                                 >  1:1808986/1‑150 (MQ=255)
                                                     aaCTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAAcc                                                                                                >  2:1827657/1‑150 (MQ=255)
                                                        tGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctgc                                                                                             <  1:3782023/150‑1 (MQ=255)
                                                         ggTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctgct                                                                                            <  1:1238931/150‑1 (MQ=255)
                                                           tGAATCAGGGAGAGTCCGGTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTgg                                                                                          <  2:1331784/150‑1 (MQ=255)
                                                            gAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTGCGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggc                                                                                         <  2:1048858/150‑1 (MQ=255)
                                                              aTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcgg                                                                                       >  2:1631218/1‑150 (MQ=255)
                                                                cAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCt                                                                                     >  1:343329/1‑150 (MQ=255)
                                                                cAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCt                                                                                     <  1:2591080/150‑1 (MQ=255)
                                                                cAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCt                                                                                     <  1:1771769/150‑1 (MQ=255)
                                                                cAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCt                                                                                     <  2:4602517/150‑1 (MQ=255)
                                                                  gggAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTAt                                                                                   >  2:4956356/1‑150 (MQ=255)
                                                                  gggAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCGGGCGGCTAt                                                                                   >  2:744321/1‑150 (MQ=255)
                                                                    gagaGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCCGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCGGGCGGCTATCt                                                                                 >  2:1260216/1‑150 (MQ=255)
                                                                    gagaGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCt                                                                                 >  2:1260395/1‑150 (MQ=255)
                                                                     agagTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTg                                                                                >  1:4538189/1‑150 (MQ=255)
                                                                       agTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGc                                                                              >  1:3577477/1‑150 (MQ=255)
                                                                          ccGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAgg                                                                           >  2:956705/1‑150 (MQ=255)
                                                                          ccGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAgg                                                                           <  1:42465/150‑1 (MQ=255)
                                                                            gTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGAc                                                                         >  1:704469/1‑150 (MQ=255)
                                                                             ttGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACa                                                                        >  2:131238/1‑150 (MQ=255)
                                                                              tGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAg                                                                       <  1:2298522/150‑1 (MQ=255)
                                                                               gCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGc                                                                      <  1:1403050/150‑1 (MQ=255)
                                                                               gCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGc                                                                      <  1:2109912/150‑1 (MQ=255)
                                                                               gCCAACTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACACCCTGCTGGCGGCTATCTGGCAGGACAGc                                                                      <  2:3162004/150‑1 (MQ=255)
                                                                                ccATCTGCGAACAAATTGGTGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCa                                                                     <  1:1845221/150‑1 (MQ=255)
                                                                                 cATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCaa                                                                    >  1:198334/1‑150 (MQ=255)
                                                                                 cATCTGCGAACAAAGTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCaa                                                                    <  1:4724980/150‑1 (MQ=255)
                                                                                     tGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGt                                                                <  1:4283309/150‑1 (MQ=255)
                                                                                      gCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGtt                                                               >  1:1207171/1‑150 (MQ=255)
                                                                                      gCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGtt                                                               <  2:1577001/150‑1 (MQ=255)
                                                                                        gAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTc                                                             <  1:2086492/150‑1 (MQ=255)
                                                                                         aaCAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCg                                                            <  1:3013344/150‑1 (MQ=255)
                                                                                                tGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAcgc                                                     >  2:2264733/1‑150 (MQ=255)
                                                                                                     aTCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGc                                                >  2:2827753/1‑150 (MQ=255)
                                                                                                      tCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCt                                               >  1:3283347/1‑150 (MQ=255)
                                                                                                      tCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCt                                               >  1:4073591/1‑150 (MQ=255)
                                                                                                      tCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCt                                               >  1:4190653/1‑150 (MQ=255)
                                                                                                        cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                             <  2:3753468/150‑1 (MQ=255)
                                                                                                             aCCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATc                                        <  1:1900275/150‑1 (MQ=255)
                                                                                                               cAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAg                                      >  2:3271717/1‑150 (MQ=255)
                                                                                                               cAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAg                                      >  2:651926/1‑150 (MQ=255)
                                                                                                                aGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGt                                     >  2:2350368/1‑150 (MQ=255)
                                                                                                                aGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGt                                     >  2:1862674/1‑150 (MQ=255)
                                                                                                                   aaaGGTCCGGTTGAGCTCAAAGTTGTGAGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAAACTGCTGGCGGCTATCTGGCACGACAGCAAAGGTTTCGGCTACGCGACGCTGGATAACAGTTcc                                  >  2:4422127/1‑150 (MQ=255)
                                                                                                                   aaaGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTcc                                  >  1:280118/1‑150 (MQ=255)
                                                                                                                       gTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGc                              >  2:4489559/1‑150 (MQ=255)
                                                                                                                        tCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCg                             <  1:1096571/150‑1 (MQ=255)
                                                                                                                            gTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGCTGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGtttt                         <  2:1801672/150‑1 (MQ=255)
                                                                                                                            gTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGtttt                         <  1:1022509/150‑1 (MQ=255)
                                                                                                                             ttGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTc                        <  2:2775871/150‑1 (MQ=255)
                                                                                                                             ttGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTc                        >  1:1414450/1‑150 (MQ=255)
                                                                                                                              tGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACACCCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTACGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCg                       >  2:2718897/1‑150 (MQ=255)
                                                                                                                                  cgcAATGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTg                   >  2:629038/1‑150 (MQ=255)
                                                                                                                                   gcAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGa                  <  1:3110082/150‑1 (MQ=255)
                                                                                                                                    cAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAg                 <  2:3836539/150‑1 (MQ=255)
                                                                                                                                    cAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAg                 <  2:1207211/150‑1 (MQ=255)
                                                                                                                                       aGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGTCTACGCGACGCTGGATATCAGTTCCGGGCGGTTTCGCCTGAGCGa              >  2:1481616/1‑150 (MQ=255)
                                                                                                                                       aGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGa              >  2:786965/1‑150 (MQ=255)
                                                                                                                                        gTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGaa             <  1:89710/150‑1 (MQ=255)
                                                                                                                                          tgtgCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAAcc           <  1:1491643/150‑1 (MQ=255)
                                                                                                                                            tgCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCgg         >  1:2359069/1‑150 (MQ=255)
                                                                                                                                             gCGTATAGTGACGCCAGGCCCCATCCGCGATGAAGCCCTGGTGCAGGAGCGTCAGGACAACCGGCTGGCGGCGATCTGGCAGGACAGCAACGGTTTCGGCTACGCGACGCTGGATATCCGTTCCGGGCGTTTTCGCCTGAGCGAACCGGc        <  2:3644090/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGACCG  >  CP009273/2850643‑2850939

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: