Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,097,262 C→T 88.3% F297F (TTC→TTT)  mutY → adenine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,097,2620CT88.3% 698.7 / 71.0 240F297F (TTC→TTT) mutYadenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (19/9);  new base T (107/105);  total (126/114)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

CGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAC  >  CP009273/3097116‑3097408
                                                                                                                                                  |                                                                                                                                                  
  tCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACtt                                                                                                                                               <  1:1003772/150‑1 (MQ=255)
  tCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACtt                                                                                                                                               <  2:1535304/150‑1 (MQ=255)
    cgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAg                                                                                                                                             <  1:2180945/150‑1 (MQ=255)
     gccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGa                                                                                                                                            >  2:4615313/1‑150 (MQ=255)
      ccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATATCCTGACGCAACTGACCGCGTTTCGGCATTCCTTCAGCCATTTTCACTTAgtt                                                                                                                                           <  2:181334/150‑3 (MQ=255)
       cgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGata                                                                                                                                          <  2:3169951/150‑1 (MQ=255)
         aGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATGCCTTCAGCCATTTTCACTTAGATAtt                                                                                                                                        <  2:2538251/150‑1 (MQ=255)
         aGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATAtt                                                                                                                                        <  1:644332/150‑1 (MQ=255)
         aGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATAtt                                                                                                                                        <  2:3733719/150‑1 (MQ=255)
         aGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATAtt                                                                                                                                        <  2:2499668/150‑1 (MQ=255)
            ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                     <  1:4394103/150‑1 (MQ=255)
            ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                     <  1:1938345/150‑1 (MQ=255)
            ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                     <  1:691613/150‑1 (MQ=255)
            ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                     >  2:1770838/1‑150 (MQ=255)
            ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                     <  2:3003123/150‑1 (MQ=255)
              aTTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGcc                                                                                                                                   >  1:4041614/1‑150 (MQ=255)
              aTTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGcc                                                                                                                                   <  1:256224/150‑1 (MQ=255)
               ttGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCt                                                                                                                                  >  1:1821495/1‑150 (MQ=255)
                   ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtgt                                                                                                                              <  1:433189/150‑1 (MQ=255)
                   ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtgt                                                                                                                              <  2:2612816/150‑1 (MQ=255)
                   ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTAtgt                                                                                                                              <  2:4442908/150‑1 (MQ=255)
                     ggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTgg                                                                                                                            >  1:2416245/1‑150 (MQ=255)
                       cGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCt                                                                                                                          >  1:2901518/1‑150 (MQ=255)
                       cGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCt                                                                                                                          >  1:3747661/1‑150 (MQ=255)
                       cGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCt                                                                                                                          <  1:3978101/150‑1 (MQ=255)
                         gTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTc                                                                                                                        <  1:2025734/150‑1 (MQ=255)
                          tttATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTcc                                                                                                                       <  1:3706171/150‑1 (MQ=255)
                           ttATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTccc                                                                                                                      >  2:2582445/1‑150 (MQ=255)
                               aCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCgtgt                                                                                                                  >  2:2760244/1‑150 (MQ=255)
                                cTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtc                                                                                                                 <  2:3918829/150‑1 (MQ=255)
                                atGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTTACGCAACTGACCGCGTTTCGGCTTACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTTGCTTCCCGTgtc                                                                                                                 <  2:668035/149‑1 (MQ=255)
                                    ttCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCa                                                                                                             >  2:3486491/1‑150 (MQ=255)
                                     tCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCAt                                                                                                            >  1:1482863/1‑150 (MQ=255)
                                      cccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCAtt                                                                                                           >  2:3094590/1‑150 (MQ=255)
                                      cccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCAtt                                                                                                           >  2:3094559/1‑150 (MQ=255)
                                       ccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTc                                                                                                          >  1:804433/1‑150 (MQ=255)
                                         gCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCAc                                                                                                        >  2:2974316/1‑150 (MQ=255)
                                            gTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCgg                                                                                                     >  1:496093/1‑150 (MQ=255)
                                             tttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGc                                                                                                    <  2:3265307/150‑1 (MQ=255)
                                             tttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGc                                                                                                    >  2:3443204/1‑150 (MQ=255)
                                              ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCt                                                                                                   <  2:368472/150‑1 (MQ=255)
                                                  cgacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                                               >  2:3606234/1‑150 (MQ=255)
                                                  cgacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                                               >  1:1560851/1‑150 (MQ=255)
                                                  cgacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                                               >  1:1102378/1‑150 (MQ=255)
                                                  cgacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGAGTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                                               >  2:4552429/1‑150 (MQ=255)
                                                     cgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGa                                                                                            >  2:3817483/1‑150 (MQ=255)
                                                      gaagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGAt                                                                                           >  1:4112833/1‑150 (MQ=255)
                                                       aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                                                          <  1:4559776/150‑1 (MQ=255)
                                                       aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                                                          >  1:1260285/1‑150 (MQ=255)
                                                       aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                                                          <  1:2491531/150‑1 (MQ=255)
                                                       aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                                                          >  1:1310625/1‑150 (MQ=255)
                                                       aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                                                          >  2:1189877/1‑150 (MQ=255)
                                                        agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                         >  2:159742/1‑150 (MQ=255)
                                                        agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                         >  2:2324892/1‑150 (MQ=255)
                                                           aaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAgg                                                                                      <  2:4371728/150‑1 (MQ=255)
                                                            aGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGc                                                                                     <  1:2700361/150‑1 (MQ=255)
                                                             gTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCa                                                                                    <  2:3191998/150‑1 (MQ=255)
                                                               ttGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAAt                                                                                  <  1:2815555/150‑1 (MQ=255)
                                                                tGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATg                                                                                 >  2:4122575/1‑150 (MQ=255)
                                                                    gCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGct                                                                             >  1:2569204/1‑150 (MQ=255)
                                                                    gCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGct                                                                             >  1:3161908/1‑150 (MQ=255)
                                                                      aGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGctct                                                                           <  1:3080917/150‑1 (MQ=255)
                                                                       gTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTg                                                                          >  2:1788636/1‑150 (MQ=255)
                                                                        tggctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTgg                                                                         <  2:2219904/150‑1 (MQ=255)
                                                                                gcAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                                                 >  1:1366317/1‑150 (MQ=255)
                                                                                  aaCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAg                                                               <  1:3648076/150‑1 (MQ=255)
                                                                                  aaCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAg                                                               >  1:2505895/1‑150 (MQ=255)
                                                                                   aCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgc                                                              <  2:2870085/150‑1 (MQ=255)
                                                                                    cGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgcg                                                             >  1:1195899/1‑150 (MQ=255)
                                                                                     ggCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgcgc                                                            <  2:144350/150‑1 (MQ=255)
                                                                                       cAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCaa                                                          <  1:2449185/150‑1 (MQ=255)
                                                                                         gATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAcc                                                        >  2:4299975/1‑150 (MQ=255)
                                                                                            tgctgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgcc                                                     >  2:708516/1‑150 (MQ=255)
                                                                                            tgctgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgcc                                                     <  1:1681639/150‑1 (MQ=255)
                                                                                             gctgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgccg                                                    <  2:3516208/150‑1 (MQ=255)
                                                                                               tgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTc                                                  <  1:250094/150‑1 (MQ=255)
                                                                                                 ccGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAg                                                <  1:3534928/150‑1 (MQ=255)
                                                                                                   gATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                                              >  1:1168736/1‑150 (MQ=255)
                                                                                                   gATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                                              >  2:3339834/1‑150 (MQ=255)
                                                                                                      aaCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGc                                           <  2:1149399/150‑1 (MQ=255)
                                                                                                         cTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTa                                        <  1:75006/150‑1 (MQ=255)
                                                                                                           gACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGc                                      >  1:3766777/1‑150 (MQ=255)
                                                                                                            aCGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCg                                     <  2:3412089/150‑1 (MQ=255)
                                                                                                             cGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCgg                                    >  2:4631764/1‑150 (MQ=255)
                                                                                                             cGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCgg                                    >  1:2049606/1‑150 (MQ=255)
                                                                                                                 aCTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTcc                                >  2:3587708/1‑150 (MQ=255)
                                                                                                                  cTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTccc                               <  1:530301/150‑1 (MQ=255)
                                                                                                                  cTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTccc                               <  1:1084768/150‑1 (MQ=255)
                                                                                                                   tGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCACTGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCg                              >  1:479466/1‑150 (MQ=255)
                                                                                                                     aCCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTg                            >  1:2053354/1‑150 (MQ=255)
                                                                                                                             tCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGttt                    >  1:2980836/1‑150 (MQ=255)
                                                                                                                              cGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTg                   >  2:398701/1‑150 (MQ=255)
                                                                                                                              cGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTg                   >  2:399116/1‑150 (MQ=255)
                                                                                                                                   tACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAca              >  1:1884921/1‑150 (MQ=255)
                                                                                                                                        tCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGt         >  1:1725714/1‑150 (MQ=255)
                                                                                                                                        tCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCACTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGt         <  2:1988042/150‑1 (MQ=255)
                                                                                                                                        tCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGt         >  1:3511405/1‑150 (MQ=255)
                                                                                                                                         cAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGtt        <  2:2363723/150‑1 (MQ=255)
                                                                                                                                               ttttCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAc  >  2:1065798/1‑150 (MQ=255)
                                                                                                                                               tttCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAc  >  2:3509480/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAC  >  CP009273/3097116‑3097408

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