Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,186,393 G→A 23.3% I350I (ATC→ATT)  purB ← adenylosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,186,3930GA23.3% 659.8 / 201.8 305I350I (ATC→ATT) purBadenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base G (117/117);  new base A (32/39);  total (149/156)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

TTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCGGGTTAACTTT  >  CP009273/1186244‑1186542
                                                                                                                                                     |                                                                                                                                                     
ttttCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCg                                                                                                                                                       >  2:4673306/1‑150 (MQ=255)
ttttCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCg                                                                                                                                                       >  1:2677589/1‑150 (MQ=255)
       tGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCa                                                                                                                                                >  1:4591656/1‑150 (MQ=255)
        gCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCAc                                                                                                                                               <  2:402282/150‑1 (MQ=255)
        gCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCAc                                                                                                                                               <  2:3035709/150‑1 (MQ=255)
            tAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCg                                                                                                                                           >  1:2938782/1‑150 (MQ=255)
            tAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCg                                                                                                                                           <  1:2024932/150‑1 (MQ=255)
            tAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCg                                                                                                                                           >  1:4857697/1‑150 (MQ=255)
            tAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCg                                                                                                                                           >  1:3430895/1‑150 (MQ=255)
               cGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAgg                                                                                                                                        >  1:4495952/1‑150 (MQ=255)
                gACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGt                                                                                                                                       >  1:1369809/1‑150 (MQ=255)
                gACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGt                                                                                                                                       <  2:4539334/150‑1 (MQ=255)
                  cGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTa                                                                                                                                     >  2:2838900/1‑150 (MQ=255)
                  cGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTAAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTa                                                                                                                                     <  2:709153/150‑1 (MQ=255)
                   gCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTAc                                                                                                                                    <  2:1193378/150‑1 (MQ=255)
                   gCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTAc                                                                                                                                    >  1:3224224/1‑150 (MQ=255)
                     aTAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACgc                                                                                                                                  >  2:3273956/1‑150 (MQ=255)
                      tAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACgca                                                                                                                                 <  2:1330113/150‑1 (MQ=255)
                      tAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACgca                                                                                                                                 <  2:1155050/150‑1 (MQ=255)
                      tAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACgca                                                                                                                                 >  2:4475033/1‑150 (MQ=255)
                              tGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTa                                                                                                                         <  2:3501787/150‑1 (MQ=255)
                              tGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTa                                                                                                                         <  2:917549/150‑1 (MQ=255)
                              tGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTa                                                                                                                         <  1:389712/150‑1 (MQ=255)
                                  ttGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGt                                                                                                                     >  1:1254457/1‑150 (MQ=255)
                                   tGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTc                                                                                                                    >  1:2897340/1‑150 (MQ=255)
                                       tCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCggtc                                                                                                                <  1:754824/150‑1 (MQ=255)
                                        cagccagcACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCggtca                                                                                                               <  1:3603328/150‑1 (MQ=255)
                                        cagccagcACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCggtca                                                                                                               >  1:1621518/1‑150 (MQ=255)
                                        cagccagcACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCggtca                                                                                                               >  1:400829/1‑150 (MQ=255)
                                                aCTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGc                                                                                                       >  2:1612018/1‑150 (MQ=255)
                                                 cTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCt                                                                                                      >  1:1759256/1‑150 (MQ=255)
                                                       aGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGc                                                                                                <  2:3163508/150‑1 (MQ=255)
                                                       aGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGc                                                                                                >  2:2912042/1‑150 (MQ=255)
                                                         ttGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCgg                                                                                              <  2:2208665/150‑1 (MQ=255)
                                                          tgtgATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCggg                                                                                             <  2:3435161/150‑1 (MQ=255)
                                                             gATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGaaa                                                                                          <  2:1294348/150‑1 (MQ=255)
                                                                ccAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCg                                                                                       >  1:1089156/1‑150 (MQ=255)
                                                                ccAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCg                                                                                       >  1:672471/1‑150 (MQ=255)
                                                                       aTCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAGCCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGttt                                                                                >  2:3142103/1‑150 (MQ=255)
                                                                       aTCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGttt                                                                                <  1:363786/150‑1 (MQ=255)
                                                                       aTCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGttt                                                                                <  2:2391538/150‑1 (MQ=255)
                                                                        tCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTttg                                                                               >  2:1810315/1‑150 (MQ=255)
                                                                           agcagATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTttgctt                                                                            <  2:166849/150‑1 (MQ=255)
                                                                             cagATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTttgcttgc                                                                          >  1:928593/1‑150 (MQ=255)
                                                                             cagATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTttgcttgc                                                                          >  2:1085221/1‑150 (MQ=255)
                                                                                aTGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAg                                                                       <  1:2304392/150‑1 (MQ=255)
                                                                                aTGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAg                                                                       <  2:1473875/150‑1 (MQ=255)
                                                                                 tGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGa                                                                      >  1:4851405/1‑150 (MQ=255)
                                                                                   gTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGAtg                                                                    <  2:782089/150‑1 (MQ=255)
                                                                                     cACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGAtgct                                                                  <  1:4537354/150‑1 (MQ=255)
                                                                                          ttCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAAt                                                             <  2:152000/150‑1 (MQ=255)
                                                                                            cACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGTAGTTTGCTTGCCAGATGCTGCAATAc                                                           <  2:4303876/150‑1 (MQ=255)
                                                                                             aCTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATAcc                                                          >  1:4774810/1‑150 (MQ=255)
                                                                                             aCTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATAcc                                                          >  1:492618/1‑150 (MQ=255)
                                                                                               ttCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACcgc                                                        <  1:208281/150‑1 (MQ=255)
                                                                                               ttCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGAGACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACcgc                                                        <  1:2638757/150‑1 (MQ=255)
                                                                                                 ccAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGt                                                      >  1:4758953/1‑150 (MQ=255)
                                                                                                    gTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTTCTTGCCAGATGCTGCAATACCGCGTTgg                                                   >  2:2518396/1‑150 (MQ=255)
                                                                                                    gTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTgg                                                   <  1:50906/150‑1 (MQ=255)
                                                                                                      ttGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGaa                                                 <  1:3036534/150‑1 (MQ=255)
                                                                                                      ttGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGaa                                                 <  1:3036181/150‑1 (MQ=255)
                                                                                                       tGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGaaa                                                <  1:627902/150‑1 (MQ=255)
                                                                                                       tGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGaaa                                                >  2:2740596/1‑150 (MQ=255)
                                                                                                        gCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAg                                               <  2:4417989/150‑1 (MQ=255)
                                                                                                        gCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAg                                               <  1:4735992/150‑1 (MQ=255)
                                                                                                         cTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCCGCGGGAAACCGGCCGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAgg                                              >  1:4712584/1‑150 (MQ=255)
                                                                                                         cTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAgg                                              >  2:3130770/1‑150 (MQ=255)
                                                                                                         cTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAgg                                              >  1:2259799/1‑150 (MQ=255)
                                                                                                         cTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAgg                                              <  2:482358/150‑1 (MQ=255)
                                                                                                           cACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGcc                                            <  2:1645727/150‑1 (MQ=255)
                                                                                                              gCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAg                                         <  2:4391428/150‑1 (MQ=255)
                                                                                                                  ttCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTc                                     <  2:2299846/150‑1 (MQ=255)
                                                                                                                   tCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTcc                                    >  1:2533018/1‑150 (MQ=255)
                                                                                                                     aGGGTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCt                                  >  1:1855916/1‑150 (MQ=255)
                                                                                                                       ggTGGATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTc                                <  1:2521004/150‑1 (MQ=255)
                                                                                                                          tgATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGa                             >  1:1471303/2‑150 (MQ=255)
                                                                                                                          ggATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGa                             <  2:353957/150‑1 (MQ=255)
                                                                                                                          ggATTGATATGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGa                             >  1:979673/1‑150 (MQ=255)
                                                                                                                           gATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAg                            <  1:3280209/150‑1 (MQ=255)
                                                                                                                           gATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAg                            <  1:4666896/150‑1 (MQ=255)
                                                                                                                              tGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGttt                         <  2:2669440/150‑1 (MQ=255)
                                                                                                                                 tatGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCCGGTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCg                      >  2:543163/1‑150 (MQ=255)
                                                                                                                                  atGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCGGCAATACCGCGGTGGAAAGGCCCAGATTCCCTTCGGGGTTTTCGa                     >  1:1474421/1‑150 (MQ=255)
                                                                                                                                  atGCAATCAAGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGa                     <  1:4161249/150‑1 (MQ=255)
                                                                                                                                      aaTCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCGTCGGAGTTTTAGAAGTc                 >  1:1011662/1‑150 (MQ=255)
                                                                                                                                       aTCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCCGCGGGAAACCGGCAGTTTGCTTGCCAGATGCGGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCg                >  2:3932714/1‑150 (MQ=255)
                                                                                                                                       aTCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCg                >  1:4412966/1‑150 (MQ=255)
                                                                                                                                       aTCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCg                <  2:4851078/150‑1 (MQ=255)
                                                                                                                                       aTCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCg                >  1:2595657/1‑150 (MQ=255)
                                                                                                                                          aaGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATc             <  1:2839443/150‑1 (MQ=255)
                                                                                                                                          aaGGCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATc             <  1:4383848/150‑1 (MQ=255)
                                                                                                                                            ggCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCgg           >  2:377712/1‑150 (MQ=255)
                                                                                                                                             gCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCGTGCCAGATGCTGCAATACCGCGGTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGCTCggg          >  1:2921555/1‑150 (MQ=255)
                                                                                                                                             gCATAACCAATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCggg          <  2:3440824/150‑1 (MQ=255)
                                                                                                                                              cATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCGGGt         <  2:662316/150‑1 (MQ=255)
                                                                                                                                               aTAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCGGGtt        >  2:4205415/1‑150 (MQ=255)
                                                                                                                                               aTAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCGGGtt        <  1:352117/150‑1 (MQ=255)
                                                                                                                                                 aaCCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCGGGTTaa      >  1:4053641/1‑150 (MQ=255)
                                                                                                                                                    cGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCGGGTTAACtt   >  2:590577/1‑150 (MQ=255)
                                                                                                                                                     gATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGCTCGGGTTAACtgt  >  2:1860522/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTTCGATGCCATAGCGACGCATAACTGTCTGGATTGGTTCAGCCAGCACTTCCCAGTTGTGATCCAGTTCATCCAGCAGATGGTCACGGTTCACTTCCAGTTTGCTCACGCCTTTCAGGGTGGATTGATATGCAATCAAGGCATAACCGATACCCACGCCGAGGTTACGCAGCACGGTAGAGTCGGTCAGGTCACGCTGCCAGCGGGAAACCGGCAGTTTGCTTGCCAGATGCTGCAATACCGCGTTGGAAAGGCCCAGATTCCCTTCGGAGTTTTCGAAGTCGATCGGGTTAACTTT  >  CP009273/1186244‑1186542

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