Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,850,790 C→A 91.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,7900CA91.4% 942.1 / 76.4 315I113I (ATC→ATAmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (17/10);  new base A (149/139);  total (166/149)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 317 total aligned reads displayed.

GAGCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGAC  >  CP009273/2850641‑2850937
                                                                                                                                                     |                                                                                                                                                   
gagCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATc                                                                                                                                                     <  1:4108005/150‑1 (MQ=255)
   ccGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGtt                                                                                                                                                  <  2:1736108/150‑1 (MQ=255)
   ccGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGtt                                                                                                                                                  <  2:2096589/150‑1 (MQ=255)
    cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                                 >  1:2951588/1‑150 (MQ=255)
    cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                                 >  1:2546333/1‑150 (MQ=255)
    cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                                 <  1:750651/150‑1 (MQ=255)
     gATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTAc                                                                                                                                                <  2:1982539/150‑1 (MQ=255)
         ccGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCa                                                                                                                                            >  2:2383817/1‑150 (MQ=255)
          cGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCACAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAg                                                                                                                                           >  2:3767601/1‑150 (MQ=255)
          cGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAg                                                                                                                                           >  1:2243910/1‑150 (MQ=255)
             tGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCa                                                                                                                                        >  1:1976446/1‑150 (MQ=255)
                cGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCa                                                                                                                                     >  2:3235887/1‑150 (MQ=255)
                  gggATTCCCTACCATGCGGTGGAAAACTATCTAGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATc                                                                                                                                   <  2:2408382/150‑1 (MQ=255)
                   ggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCa                                                                                                                                  <  2:4317670/150‑1 (MQ=255)
                   ggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCa                                                                                                                                  <  2:849237/150‑1 (MQ=255)
                    gATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAg                                                                                                                                 >  2:1909095/1‑150 (MQ=255)
                    gATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAg                                                                                                                                 <  2:391225/150‑1 (MQ=255)
                    gATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAg                                                                                                                                 <  1:3969846/150‑1 (MQ=255)
                    gATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAg                                                                                                                                 >  1:3201780/1‑150 (MQ=255)
                     aTTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGc                                                                                                                                >  2:100317/1‑150 (MQ=255)
                     aTTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGc                                                                                                                                >  2:4184090/1‑150 (MQ=255)
                       tCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGa                                                                                                                              <  1:1059483/150‑1 (MQ=255)
                         ccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATg                                                                                                                            >  1:3857241/1‑150 (MQ=255)
                            aCCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAg                                                                                                                         >  1:3101929/1‑150 (MQ=255)
                             ccATGCGGTGGAAAACTATCTCGCCAACCTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGc                                                                                                                        <  2:687878/150‑1 (MQ=255)
                              cATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGcc                                                                                                                       >  2:2139728/1‑150 (MQ=255)
                                tGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGACCAAATTGGCGATCCGGCGACCCGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTCTAGTTACGCCAGGCACCATCAGCGATGAAGCCCt                                                                                                                     >  2:1858535/1‑150 (MQ=255)
                                 gCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTg                                                                                                                    <  1:1154291/150‑1 (MQ=255)
                                 gCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTg                                                                                                                    <  2:25414/150‑1 (MQ=255)
                                  cGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGt                                                                                                                   >  2:330146/1‑150 (MQ=255)
                                   ggTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGtt                                                                                                                  >  1:3562808/1‑150 (MQ=255)
                                       gAAATCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAg                                                                                                              <  2:3585430/150‑1 (MQ=255)
                                        aaaaCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAgg                                                                                                             <  1:179223/150‑1 (MQ=255)
                                         aaaCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGa                                                                                                            >  2:2033314/1‑150 (MQ=255)
                                          aaCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAg                                                                                                           >  2:1054909/1‑150 (MQ=255)
                                          aaCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAg                                                                                                           >  2:4463273/1‑150 (MQ=255)
                                            cTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCg                                                                                                         >  1:3316092/1‑150 (MQ=255)
                                             tATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGACATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGt                                                                                                        >  2:730471/1‑150 (MQ=255)
                                                          tGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctgc                                                                                           >  2:1650337/1‑150 (MQ=255)
                                                           ggTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctgct                                                                                          <  2:1092810/150‑1 (MQ=255)
                                                            gTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctgctg                                                                                         <  1:398510/150‑1 (MQ=255)
                                                             tGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTgg                                                                                        >  1:2834194/1‑150 (MQ=255)
                                                                aTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcgg                                                                                     >  2:2263797/1‑150 (MQ=255)
                                                                aTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcgg                                                                                     >  2:3795677/1‑150 (MQ=255)
                                                                aTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcgg                                                                                     >  1:2286569/1‑150 (MQ=255)
                                                                   aGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTa                                                                                  >  1:2033577/1‑150 (MQ=255)
                                                                    gggAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTAt                                                                                 <  2:4560072/150‑1 (MQ=255)
                                                                     ggAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATc                                                                                >  2:4524833/1‑150 (MQ=255)
                                                                      gagaGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCt                                                                               <  1:813250/150‑1 (MQ=255)
                                                                      gagaGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCt                                                                               <  1:4473006/150‑1 (MQ=255)
                                                                       agagTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTg                                                                              <  1:2240709/150‑1 (MQ=255)
                                                                        gagTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTgg                                                                             <  1:4254258/150‑1 (MQ=255)
                                                                            ccGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAgg                                                                         >  2:1610862/1‑150 (MQ=255)
                                                                              gTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGAc                                                                       >  2:4711752/1‑150 (MQ=255)
                                                                                tGCCATCTGCGAACACAGTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAg                                                                     <  2:650870/150‑1 (MQ=255)
                                                                                     tCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAg                                                                >  1:1315799/1‑150 (MQ=255)
                                                                                     tCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAg                                                                <  1:3063663/150‑1 (MQ=255)
                                                                                      cTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAgg                                                               <  1:2773754/150‑1 (MQ=255)
                                                                                      cTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAgg                                                               <  2:3201780/150‑1 (MQ=255)
                                                                                      cTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAgg                                                               <  2:1845274/150‑1 (MQ=255)
                                                                                      cTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAgg                                                               <  2:1163573/150‑1 (MQ=255)
                                                                                         cGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGttt                                                            >  1:3723106/1‑150 (MQ=255)
                                                                                          gAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTc                                                           <  2:462399/150‑1 (MQ=255)
                                                                                          gAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTc                                                           <  1:4033264/150‑1 (MQ=255)
                                                                                            aCAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCgg                                                         <  1:884986/150‑1 (MQ=255)
                                                                                             cAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGc                                                        >  2:144237/1‑150 (MQ=255)
                                                                                              aaaTTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCt                                                       <  1:1823123/150‑1 (MQ=255)
                                                                                                  tGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAcgc                                                   >  1:1006868/1‑150 (MQ=255)
                                                                                                  tGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAcgc                                                   >  1:1905905/1‑150 (MQ=255)
                                                                                                     cGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGAc                                                <  1:441762/150‑1 (MQ=255)
                                                                                                      gATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACg                                               <  2:2233315/150‑1 (MQ=255)
                                                                                                      gATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACg                                               >  2:1821854/1‑150 (MQ=255)
                                                                                                        tCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCt                                             <  1:2595541/150‑1 (MQ=255)
                                                                                                          cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGCCAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                           <  1:1487826/150‑1 (MQ=255)
                                                                                                          cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                           >  2:1060917/1‑150 (MQ=255)
                                                                                                            gCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGat                                         <  2:2600651/150‑1 (MQ=255)
                                                                                                              gACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGatat                                       >  2:3700748/1‑150 (MQ=255)
                                                                                                              gACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGatat                                       >  2:1378114/1‑150 (MQ=255)
                                                                                                                ccAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCa                                     <  1:1738133/150‑1 (MQ=255)
                                                                                                                ccAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCa                                     <  1:638448/150‑1 (MQ=255)
                                                                                                                  aGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGt                                   >  1:1439617/1‑150 (MQ=255)
                                                                                                                     aaaGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTcc                                <  1:3900297/150‑1 (MQ=255)
                                                                                                                      aaGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCg                               >  1:68791/1‑150 (MQ=255)
                                                                                                                       aGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCgg                              >  1:1645603/1‑150 (MQ=255)
                                                                                                                         gTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGc                            >  2:2825995/1‑150 (MQ=255)
                                                                                                                             ggTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGttt                        >  1:549141/1‑150 (MQ=255)
                                                                                                                             ggTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGttt                        <  1:4508408/150‑1 (MQ=255)
                                                                                                                               ttGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTc                      <  2:2951588/150‑1 (MQ=255)
                                                                                                                                 gAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGc                    >  1:660645/1‑150 (MQ=255)
                                                                                                                                  aGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGcc                   >  2:3152243/1‑150 (MQ=255)
                                                                                                                                    cgcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTg                 <  2:2650142/150‑1 (MQ=255)
                                                                                                                                      cAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAg               >  1:2994040/1‑150 (MQ=255)
                                                                                                                                        aaGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCg             >  2:1943115/1‑150 (MQ=255)
                                                                                                                                        aaGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGATCCGGGCGTGTTCGCCTGAGCg             >  1:4579514/1‑150 (MQ=255)
                                                                                                                                            tgtgCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAAcc         >  2:128398/1‑150 (MQ=255)
                                                                                                                                            tgtgCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAAcc         <  1:281906/150‑1 (MQ=255)
                                                                                                                                            tgtgCGTATAGTGACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAAcc         <  2:3946230/150‑1 (MQ=255)
                                                                                                                                             gtgCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCg        <  1:2139728/150‑1 (MQ=255)
                                                                                                                                               gCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGc      <  2:2834194/150‑1 (MQ=255)
                                                                                                                                                 gTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTg    >  1:4301354/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAGCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTGAC  >  CP009273/2850641‑2850937

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: