Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 2,850,790 C→A 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,7900CA99.6% 744.9 / ‑4.9 226I113I (ATC→ATAmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (114/111);  minor base T (0/1);  total (114/112)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 226 total aligned reads displayed.

GCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTG  >  CP009273/2850643‑2850935
                                                                                                                                                   |                                                                                                                                                 
 ccGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGtt                                                                                                                                                <  2:3000535/150‑1 (MQ=255)
  cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                               >  1:2577230/1‑150 (MQ=255)
  cGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTa                                                                                                                                               >  2:868691/1‑150 (MQ=255)
    aTCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACg                                                                                                                                             >  2:915427/1‑150 (MQ=255)
        cGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAg                                                                                                                                         <  2:2592376/150‑1 (MQ=255)
         gATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAgg                                                                                                                                        <  2:1780819/150‑1 (MQ=255)
          aTGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGc                                                                                                                                       <  1:2915353/150‑1 (MQ=255)
            ggCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCAc                                                                                                                                     >  1:850463/1‑150 (MQ=255)
              cGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCa                                                                                                                                   <  2:1552492/150‑1 (MQ=255)
              cGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCa                                                                                                                                   <  1:3405974/150‑1 (MQ=255)
               ggggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCAt                                                                                                                                  >  1:2826457/1‑150 (MQ=255)
               ggggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGAGTATAGTTACGCCAGGCACCAt                                                                                                                                  >  2:1476028/1‑150 (MQ=255)
                 ggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCa                                                                                                                                <  1:1246841/150‑1 (MQ=255)
                 ggATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCa                                                                                                                                <  1:2688195/150‑1 (MQ=255)
                  gATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAg                                                                                                                               <  1:493390/150‑1 (MQ=255)
                  gATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAg                                                                                                                               <  2:2247316/150‑1 (MQ=255)
                   aTTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGc                                                                                                                              >  1:2013646/1‑150 (MQ=255)
                   aTTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGc                                                                                                                              >  1:101281/1‑150 (MQ=255)
                      cccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGTGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGAt                                                                                                                           <  2:2518728/150‑1 (MQ=255)
                      cccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGAt                                                                                                                           >  1:553908/1‑150 (MQ=255)
                      cccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGAt                                                                                                                           >  2:1736893/1‑150 (MQ=255)
                      cccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGAt                                                                                                                           >  2:54364/1‑150 (MQ=255)
                      cccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGAt                                                                                                                           <  1:1334435/150‑1 (MQ=255)
                       ccTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATg                                                                                                                          <  2:1517824/150‑1 (MQ=255)
                        cTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGGTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGa                                                                                                                         >  2:2540027/1‑150 (MQ=255)
                          aCCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAg                                                                                                                       >  1:1279821/1‑150 (MQ=255)
                           ccATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGc                                                                                                                      >  1:892783/1‑150 (MQ=255)
                           ccATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGGATAGTTACGCCAGGCACCATCAGCGATGACGc                                                                                                                      >  2:3291729/1‑150 (MQ=255)
                            cATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGcc                                                                                                                     <  2:32464/150‑1 (MQ=255)
                              tGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCt                                                                                                                   <  2:1372660/150‑1 (MQ=255)
                              tGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGGATAGTTACGCCAGGCACCATCAGCGATGAAGCCCt                                                                                                                   >  2:2497298/1‑150 (MQ=255)
                                   tGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGc                                                                                                              >  1:3638788/1‑150 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                             <  2:1585121/150‑1 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                             >  2:128951/1‑150 (MQ=255)
                                    ggAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGGATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCa                                                                                                             >  2:1842547/1‑150 (MQ=255)
                                       aaaCTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGa                                                                                                          >  2:1128484/1‑150 (MQ=255)
                                           tATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGt                                                                                                      >  1:3718687/1‑150 (MQ=255)
                                           tATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGt                                                                                                      <  2:1670630/150‑1 (MQ=255)
                                           tATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGt                                                                                                      >  1:979876/1‑150 (MQ=255)
                                            aTCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTc                                                                                                     <  1:1162351/150‑1 (MQ=255)
                                            aTCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTc                                                                                                     >  1:3431118/1‑150 (MQ=255)
                                                  ccAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACa                                                                                               <  2:3205757/150‑1 (MQ=255)
                                                   cAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACaa                                                                                              >  2:3131546/1‑150 (MQ=255)
                                                     aaCTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAAcc                                                                                            <  2:1051245/150‑1 (MQ=255)
                                                       cTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctg                                                                                          >  2:143748/1‑150 (MQ=255)
                                                        tGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACctgc                                                                                         >  2:3154523/1‑150 (MQ=255)
                                                          gTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCACTGTTGCAGGAGCGTCAGGACAACctgctg                                                                                       <  1:1877225/150‑1 (MQ=255)
                                                            gAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggc                                                                                     <  2:3573419/150‑1 (MQ=255)
                                                             aaTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcg                                                                                    >  1:2346544/1‑150 (MQ=255)
                                                             aaTCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTggcg                                                                                    >  1:3383144/1‑150 (MQ=255)
                                                                 aGGGAGATACCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGTTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGTAGCGTCAGGACAACCTGCTGGCGGCTa                                                                                >  1:2651043/1‑150 (MQ=255)
                                                                 aGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTa                                                                                >  1:745867/1‑150 (MQ=255)
                                                                  gggAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTAt                                                                               >  2:3674538/1‑150 (MQ=255)
                                                                  gggAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTAt                                                                               >  1:683681/1‑150 (MQ=255)
                                                                    gagaGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCt                                                                             <  2:1946242/150‑1 (MQ=255)
                                                                     agagTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTg                                                                            >  1:1808398/1‑150 (MQ=255)
                                                                        gTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCa                                                                         >  1:3630577/1‑150 (MQ=255)
                                                                         tCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAg                                                                        >  1:2187946/1‑150 (MQ=255)
                                                                          tcGTTGCCATCTGCGAACTAATTGGCGATCCGGCGATCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAgg                                                                       <  1:3688792/149‑1 (MQ=255)
                                                                           cGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGa                                                                      >  2:3785776/1‑150 (MQ=255)
                                                                            gTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGAc                                                                     >  2:3213087/1‑150 (MQ=255)
                                                                             ttGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACa                                                                    <  1:1784090/150‑1 (MQ=255)
                                                                                   tCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAg                                                              <  1:3497108/150‑1 (MQ=255)
                                                                                        gAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTc                                                         <  1:1377459/150‑1 (MQ=255)
                                                                                         aaCAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCg                                                        >  2:684033/1‑150 (MQ=255)
                                                                                            aaaTTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCt                                                     >  2:2773581/1‑150 (MQ=255)
                                                                                            aaaTTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCt                                                     <  2:1107595/150‑1 (MQ=255)
                                                                                             aaTTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTa                                                    >  1:2446213/1‑150 (MQ=255)
                                                                                               ttGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTAcg                                                  >  2:597620/1‑150 (MQ=255)
                                                                                                   cGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGAc                                              <  2:2103880/150‑1 (MQ=255)
                                                                                                    gATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATTGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACg                                             <  1:2607279/150‑1 (MQ=255)
                                                                                                      tCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCt                                           <  2:1502135/150‑1 (MQ=255)
                                                                                                        cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                         >  2:3529059/1‑150 (MQ=255)
                                                                                                        cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                         >  1:2349792/1‑150 (MQ=255)
                                                                                                        cGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTgg                                         >  1:2350724/1‑150 (MQ=255)
                                                                                                         ggCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGa                                        >  2:357147/1‑150 (MQ=255)
                                                                                                            gACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGatat                                     <  1:2993622/150‑1 (MQ=255)
                                                                                                            gACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGatat                                     >  2:1033865/1‑150 (MQ=255)
                                                                                                              ccAGCAACGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGGTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCa                                   <  1:3256219/150‑1 (MQ=255)
                                                                                                              ccAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCa                                   <  2:3798217/150‑1 (MQ=255)
                                                                                                               cAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAg                                  <  2:589273/150‑1 (MQ=255)
                                                                                                                aGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGt                                 <  1:155411/150‑1 (MQ=255)
                                                                                                                aGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGt                                 <  2:2397284/150‑1 (MQ=255)
                                                                                                                 gCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGtt                                <  2:3663455/150‑1 (MQ=255)
                                                                                                                 gCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGtt                                <  2:2468941/150‑1 (MQ=255)
                                                                                                                    aaGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCg                             >  2:2087083/1‑150 (MQ=255)
                                                                                                                     aGGTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCgg                            <  2:849387/150‑1 (MQ=255)
                                                                                                                       gTCCGGTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGc                          >  2:347181/1‑150 (MQ=255)
                                                                                                                            gTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGtttt                     >  2:2337771/1‑150 (MQ=255)
                                                                                                                            gTTGAGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGtttt                     >  2:2081826/1‑150 (MQ=255)
                                                                                                                                aGCGCAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGcc                 >  1:2268441/1‑150 (MQ=255)
                                                                                                                                 gcgcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCt                >  1:2726234/1‑150 (MQ=255)
                                                                                                                                  cgcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTg               >  1:1647857/1‑150 (MQ=255)
                                                                                                                                  cgcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTg               >  2:2825233/1‑150 (MQ=255)
                                                                                                                                   gcAAAGTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGa              <  1:413848/150‑1 (MQ=255)
                                                                                                                                        gTTGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGaa         >  2:3041462/1‑150 (MQ=255)
                                                                                                                                         ttGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAAc        <  1:1476028/150‑1 (MQ=255)
                                                                                                                                         ttGTGCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAAc        <  1:2070675/150‑1 (MQ=255)
                                                                                                                                             gCGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGc    >  2:2882841/1‑150 (MQ=255)
                                                                                                                                              cGTATAGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCt   <  1:3039168/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCGATCCCGATGGCGGGGATTCCCTACCATGCGGTGGAAAACTATCTCGCCAAACTGGTGAATCAGGGAGAGTCCGTTGCCATCTGCGAACAAATTGGCGATCCGGCGACCAGCAAAGGTCCGGTTGAGCGCAAAGTTGTGCGTATCGTTACGCCAGGCACCATCAGCGATGAAGCCCTGTTGCAGGAGCGTCAGGACAACCTGCTGGCGGCTATCTGGCAGGACAGCAAAGGTTTCGGCTACGCGACGCTGGATATCAGTTCCGGGCGTTTTCGCCTGAGCGAACCGGCTG  >  CP009273/2850643‑2850935

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: