Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,097,262 C→T 64.6% 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,2620CT64.6% 272.2 / 273.3 243F297F (TTC→TTT) mutYadenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (42/44);  new base T (67/90);  total (109/134)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 247 total aligned reads displayed.

AGCGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAC  >  CP009273/3097114‑3097408
                                                                                                                                                    |                                                                                                                                                  
  cGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCAc                                                                                                                                                 <  1:365003/150‑1 (MQ=255)
  cGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCAc                                                                                                                                                 >  1:328175/1‑150 (MQ=255)
  cGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCAc                                                                                                                                                 >  2:3081991/1‑150 (MQ=255)
     ccgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTa                                                                                                                                              <  2:3219873/150‑1 (MQ=255)
      cgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAg                                                                                                                                             >  1:2295474/1‑150 (MQ=255)
      cgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAg                                                                                                                                             <  1:217405/150‑1 (MQ=255)
      cgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAg                                                                                                                                             <  2:2655243/150‑1 (MQ=255)
       gccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGa                                                                                                                                            >  1:189317/1‑150 (MQ=255)
       gccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGa                                                                                                                                            >  1:191601/1‑150 (MQ=255)
          gAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGatat                                                                                                                                         <  1:340516/150‑1 (MQ=255)
           aGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATAtt                                                                                                                                        <  1:2883597/150‑1 (MQ=255)
             cGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgt                                                                                                                                      >  2:607459/1‑150 (MQ=255)
              ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                     <  1:709609/150‑1 (MQ=255)
              ggATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAAAGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                     <  2:2692832/150‑1 (MQ=255)
                aTTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACAGACCGAGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGcc                                                                                                                                   >  2:244126/1‑150 (MQ=255)
                  tgtgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTa                                                                                                                                 <  2:4324998/150‑1 (MQ=255)
                    tgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtg                                                                                                                               <  1:2746260/150‑1 (MQ=255)
                    tgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAAAGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtg                                                                                                                               >  1:2577205/1‑150 (MQ=255)
                     ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAtgt                                                                                                                              <  1:2091576/150‑1 (MQ=255)
                           gTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTc                                                                                                                        <  2:2237264/150‑1 (MQ=255)
                              tataCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCg                                                                                                                     <  2:2256834/150‑1 (MQ=255)
                              tataCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCg                                                                                                                     >  1:849349/1‑150 (MQ=255)
                              tataCTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCg                                                                                                                     >  1:849289/1‑150 (MQ=255)
                                  cTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtc                                                                                                                 <  1:1513429/150‑1 (MQ=255)
                                    gTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtcgt                                                                                                               >  1:1627515/1‑150 (MQ=255)
                                      ttCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAAAGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCa                                                                                                             >  2:3327238/1‑150 (MQ=255)
                                         ccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTc                                                                                                          <  1:554688/150‑1 (MQ=255)
                                         ccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTc                                                                                                          <  1:793034/150‑1 (MQ=255)
                                         ccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAAAGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTc                                                                                                          <  2:1698710/150‑1 (MQ=255)
                                         ccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAAAGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTc                                                                                                          <  2:1698153/150‑1 (MQ=255)
                                          cGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCa                                                                                                         <  2:809352/150‑1 (MQ=255)
                                            cAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCAcc                                                                                                       <  2:1751804/150‑1 (MQ=255)
                                             aGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCg                                                                                                      >  2:726537/1‑150 (MQ=255)
                                                ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCt                                                                                                   >  1:1080835/1‑150 (MQ=255)
                                                ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCt                                                                                                   <  1:2526431/150‑1 (MQ=255)
                                                 tGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTg                                                                                                  >  2:4050409/1‑150 (MQ=255)
                                                    cgacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                                               >  2:524029/1‑150 (MQ=255)
                                                    cgacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                                               <  1:4513610/150‑1 (MQ=255)
                                                     gacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATg                                                                                              <  2:1890573/150‑1 (MQ=255)
                                                     gacgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATg                                                                                              <  2:1890042/150‑1 (MQ=255)
                                                      acgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATgg                                                                                             <  1:3602026/150‑1 (MQ=255)
                                                      acgaAGAAAGTTTGCGGCAGTGGCTGGCGCAAAGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATgg                                                                                             <  1:913368/150‑1 (MQ=255)
                                                        gaagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGAt                                                                                           >  1:3380240/1‑150 (MQ=255)
                                                         aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                                                          >  2:782760/1‑150 (MQ=255)
                                                          agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                         <  2:3106883/150‑1 (MQ=255)
                                                          agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                         <  2:3067665/150‑1 (MQ=255)
                                                          agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                         <  2:174368/150‑1 (MQ=255)
                                                          agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                         >  2:1966870/1‑150 (MQ=255)
                                                            aaaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAg                                                                                       <  1:513196/150‑1 (MQ=255)
                                                            aaaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAg                                                                                       <  2:3334669/150‑1 (MQ=255)
                                                            aaaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAg                                                                                       <  2:602389/150‑1 (MQ=255)
                                                             aaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAgg                                                                                      >  1:57917/1‑150 (MQ=255)
                                                              aGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGGGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGc                                                                                     <  1:4476305/150‑1 (MQ=255)
                                                                    cGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCAGTCACCGGCTGCATGGATGAAGGCAATgcg                                                                               >  2:2947773/1‑150 (MQ=255)
                                                                      gCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATCTTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGct                                                                             <  2:4130330/150‑1 (MQ=255)
                                                                       cAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGctc                                                                            >  1:974598/1‑150 (MQ=255)
                                                                          tggctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTgg                                                                         >  2:3681947/1‑150 (MQ=255)
                                                                           ggctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGt                                                                        <  2:4540352/150‑1 (MQ=255)
                                                                                    aaCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAg                                                               <  2:1241271/150‑1 (MQ=255)
                                                                                         cAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCaa                                                          <  1:2032667/150‑1 (MQ=255)
                                                                                          aGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAc                                                         >  1:493008/1‑150 (MQ=255)
                                                                                          aGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAc                                                         >  1:2018308/1‑150 (MQ=255)
                                                                                          aGATTGCTGCCGATAACCTGACGAAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAc                                                         >  1:492668/1‑150 (MQ=255)
                                                                                                ctgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGt                                                   <  2:2734064/150‑1 (MQ=255)
                                                                                                 tgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTc                                                  <  1:389873/150‑1 (MQ=255)
                                                                                                  gcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCa                                                 <  1:3480035/150‑1 (MQ=255)
                                                                                                  gcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCa                                                 <  1:3479980/150‑1 (MQ=255)
                                                                                                  gcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCa                                                 <  1:3080152/150‑1 (MQ=255)
                                                                                                   ccGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAg                                                <  2:3530230/150‑1 (MQ=255)
                                                                                                   ccGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAg                                                <  2:2652481/150‑1 (MQ=255)
                                                                                                   ccGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAg                                                <  2:270207/150‑1 (MQ=255)
                                                                                                    cGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGt                                               >  2:1343220/1‑150 (MQ=255)
                                                                                                     gATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                                              >  1:4388687/1‑150 (MQ=255)
                                                                                                     gATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                                              >  1:796123/1‑150 (MQ=255)
                                                                                                      aTAACCTGACGCATCTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTg                                             <  2:1536734/150‑1 (MQ=255)
                                                                                                      aTAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTg                                             >  1:3150767/1‑150 (MQ=255)
                                                                                                        aaCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGc                                           >  2:3983325/1‑150 (MQ=255)
                                                                                                         aCCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGcc                                          >  2:2225916/1‑150 (MQ=255)
                                                                                                          ccTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCt                                         <  2:518853/150‑1 (MQ=255)
                                                                                                                gCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGc                                   >  2:1422744/1‑150 (MQ=255)
                                                                                                                 cAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCt                                  >  1:3597599/1‑150 (MQ=255)
                                                                                                                   aCTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTcc                                <  2:1799624/150‑1 (MQ=255)
                                                                                                                      gACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGt                             >  1:416631/1‑150 (MQ=255)
                                                                                                                      gACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGt                             <  2:646963/150‑1 (MQ=255)
                                                                                                                         cgcgTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGa                          <  2:2300034/150‑1 (MQ=255)
                                                                                                                            gTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCg                       >  1:597092/1‑150 (MQ=255)
                                                                                                                            gTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCg                       >  1:2117603/1‑150 (MQ=255)
                                                                                                                            gTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGAATGGATTAACGCAATGAGCTCTGGTATAACATATCGCAACCGCCGTCAGTTGGCCTATCGGCTCCCGTGGAGCg                       >  1:2940434/1‑150 (MQ=255)
                                                                                                                            gTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGACATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCg                       >  1:597405/1‑150 (MQ=255)
                                                                                                                            gTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCg                       >  2:3447168/1‑150 (MQ=255)
                                                                                                                               tCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGttt                    <  1:1395076/150‑1 (MQ=255)
                                                                                                                                  gCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCGGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTAGGCCTAGCGGCTCCCGTGGAGCGTTTGtt                 <  1:4486440/150‑1 (MQ=255)
                                                                                                                                    aTACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAc               <  1:4004180/150‑1 (MQ=255)
                                                                                                                                     tACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAca              <  1:3687945/150‑1 (MQ=255)
                                                                                                                                       ccTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAcagc            >  1:2254186/1‑150 (MQ=255)
                                                                                                                                       ccTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAcagc            >  1:3248141/1‑150 (MQ=255)
                                                                                                                                           cAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGtt        <  1:100816/150‑1 (MQ=255)
                                                                                                                                            aGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTa       >  1:500615/1‑150 (MQ=255)
                                                                                                                                              ccATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACg     <  1:2430260/150‑1 (MQ=255)
                                                                                                                                                 ttttCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAc  >  2:4314964/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCGTCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCAC  >  CP009273/3097114‑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: ‑