Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,411,867 G→T 5.7% E115* (GAA→TAA)  ulaC → L‑ascorbate‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,411,8670GT 5.7% 1327.3 / 33.3 331E115* (GAA→TAA) ulaCL‑ascorbate‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (160/152);  new base T (11/8);  total (171/160)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

CAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGAT  >  CP009273/4411718‑4412016
                                                                                                                                                     |                                                                                                                                                     
cAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAgg                                                                                                                                                       <  2:2490090/150‑1 (MQ=255)
  aTGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGTaa                                                                                                                                                     <  1:4590650/150‑1 (MQ=255)
    gCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGt                                                                                                                                                   <  1:4779440/150‑1 (MQ=255)
       gCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTggg                                                                                                                                                <  1:4516551/150‑1 (MQ=255)
        cACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGc                                                                                                                                               <  2:2897609/150‑1 (MQ=255)
          cGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGcat                                                                                                                                             <  2:4164207/150‑1 (MQ=255)
          cGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGGTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGcat                                                                                                                                             <  2:4533712/150‑1 (MQ=255)
               gTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGc                                                                                                                                        >  1:1244156/1‑150 (MQ=255)
                tCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCa                                                                                                                                       <  2:2453773/150‑1 (MQ=255)
                   ggAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGAt                                                                                                                                    <  1:1358364/150‑1 (MQ=255)
                    gAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATc                                                                                                                                   >  2:3068212/1‑150 (MQ=255)
                    gAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATc                                                                                                                                   >  1:3455775/1‑150 (MQ=255)
                    gAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATc                                                                                                                                   >  2:3108074/1‑150 (MQ=255)
                     aaGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCg                                                                                                                                  >  1:3859796/1‑150 (MQ=255)
                        agGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCa                                                                                                                               <  2:2945993/150‑1 (MQ=255)
                           gCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAAcc                                                                                                                            >  1:4030470/1‑150 (MQ=255)
                            cGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCt                                                                                                                           >  2:1443231/1‑150 (MQ=255)
                            cGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCt                                                                                                                           <  1:1095130/150‑1 (MQ=255)
                              ttAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGt                                                                                                                         >  2:2234520/1‑150 (MQ=255)
                               tAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGtt                                                                                                                        >  1:4078111/1‑150 (MQ=255)
                                   aaaaCCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaa                                                                                                                    >  2:2067211/1‑150 (MQ=255)
                                   aaaaCCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGAGGGCATCATGCAGATCGTCAACCTGTTtgaa                                                                                                                    >  2:3356915/1‑150 (MQ=255)
                                    aaaCCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaag                                                                                                                   <  1:2850746/150‑1 (MQ=255)
                                    aaaCCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaag                                                                                                                   >  2:2916040/1‑150 (MQ=255)
                                    aaaCCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaag                                                                                                                   <  2:684938/150‑1 (MQ=255)
                                           tttCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaagatgaaga                                                                                                            >  2:3515208/1‑150 (MQ=255)
                                              ctctctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGaa                                                                                                         >  1:1801803/1‑150 (MQ=255)
                                              ctctctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGaa                                                                                                         <  1:2654410/150‑1 (MQ=255)
                                                ctctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAAtt                                                                                                       <  1:4387168/150‑1 (MQ=255)
                                                ctctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAAtt                                                                                                       >  1:1552350/1‑150 (MQ=255)
                                                 tctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAAttt                                                                                                      >  2:483367/1‑150 (MQ=255)
                                                 tctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAAttt                                                                                                      <  2:3837499/150‑1 (MQ=255)
                                                  ctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAAtttt                                                                                                     >  2:5130091/1‑150 (MQ=255)
                                                   tGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTg                                                                                                    >  1:3054540/1‑150 (MQ=255)
                                                      tAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGAcc                                                                                                 <  1:109775/150‑1 (MQ=255)
                                                       aaCGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCg                                                                                                >  1:2285816/1‑150 (MQ=255)
                                                        aCGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGt                                                                                               >  2:1355059/1‑150 (MQ=255)
                                                         cGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGtt                                                                                              >  1:3499504/1‑150 (MQ=255)
                                                         cGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGtt                                                                                              >  1:491420/1‑150 (MQ=255)
                                                          gCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGttt                                                                                             <  1:1702193/150‑1 (MQ=255)
                                                          gCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGttt                                                                                             <  1:3581552/150‑1 (MQ=255)
                                                             gAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTAcg                                                                                          <  1:4931378/150‑1 (MQ=255)
                                                              aaaaaGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTAcgc                                                                                         <  1:2144020/150‑1 (MQ=255)
                                                              aaaaaGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTAcgc                                                                                         <  1:3355133/150‑1 (MQ=255)
                                                               aaaaGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTAcgcg                                                                                        <  2:4996596/150‑1 (MQ=255)
                                                                 aaGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGcc                                                                                      >  1:2858336/1‑150 (MQ=255)
                                                                  aGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCt                                                                                     >  2:2359220/1‑150 (MQ=255)
                                                                      ncTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCg                                                                                 >  1:2561293/2‑150 (MQ=255)
                                                                      gCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCg                                                                                 >  1:5103555/1‑150 (MQ=255)
                                                                       cTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGt                                                                                <  2:763104/150‑1 (MQ=255)
                                                                         ggAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTAc                                                                              >  1:1915948/1‑150 (MQ=255)
                                                                          gAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTAcc                                                                             >  1:2597356/1‑150 (MQ=255)
                                                                           aGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCg                                                                            >  2:868748/1‑150 (MQ=255)
                                                                            gTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGa                                                                           <  2:3611482/150‑1 (MQ=255)
                                                                             ttCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAg                                                                          >  2:1533386/1‑150 (MQ=255)
                                                                               cAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACGTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCa                                                                        >  1:2018815/1‑150 (MQ=255)
                                                                               cAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCa                                                                        >  2:3212750/1‑150 (MQ=255)
                                                                               cAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCa                                                                        >  1:4157157/1‑150 (MQ=255)
                                                                                aaCCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTATGACCGTTTACGCGCCTGCCGAACCGAGCAg                                                                       >  1:4722128/1‑150 (MQ=255)
                                                                                 aCCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAgg                                                                      >  1:2131852/1‑150 (MQ=255)
                                                                                  ccACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGa                                                                     <  2:1204219/150‑1 (MQ=255)
                                                                                  ccACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGa                                                                     >  1:4842627/1‑150 (MQ=255)
                                                                                     cGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGt                                                                  <  1:97543/150‑1 (MQ=255)
                                                                                       aTGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTAc                                                                >  2:1613066/1‑150 (MQ=255)
                                                                                        tGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACt                                                               >  2:55518/1‑150 (MQ=255)
                                                                                            aaCGATCCGGTGGATATCCTCATAACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGAt                                                           >  2:3748747/1‑150 (MQ=255)
                                                                                               gATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGAtctc                                                        <  1:197523/150‑1 (MQ=255)
                                                                                               gATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGAtctc                                                        >  1:4219976/1‑150 (MQ=255)
                                                                                                aTCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCa                                                       >  1:2898189/1‑150 (MQ=255)
                                                                                                     gTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGTAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGAc                                                  >  2:3445071/1‑150 (MQ=255)
                                                                                                       ggATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCg                                                <  2:3162949/150‑1 (MQ=255)
                                                                                                         atatCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCa                                              >  2:4009447/1‑150 (MQ=255)
                                                                                                          tatCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGTAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCAc                                             <  1:2186544/150‑1 (MQ=255)
                                                                                                          tatCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCAc                                             >  2:635440/1‑150 (MQ=255)
                                                                                                           atCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCAcc                                            >  2:4213746/1‑150 (MQ=255)
                                                                                                           atCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCAcc                                            >  2:1708220/1‑150 (MQ=255)
                                                                                                           atCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCAcc                                            >  1:4902863/1‑150 (MQ=255)
                                                                                                               tcatcaCCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAAcg                                        <  2:2601151/150‑1 (MQ=255)
                                                                                                                catcaCCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAAcgc                                       <  1:633178/150‑1 (MQ=255)
                                                                                                                catcaCCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAAcgc                                       <  1:3212750/150‑1 (MQ=255)
                                                                                                                         ggcggcGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTa                              <  1:1134237/150‑1 (MQ=255)
                                                                                                                         ggcggcGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTa                              <  1:761246/150‑1 (MQ=255)
                                                                                                                            ggcggTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTaaga                           >  2:1941507/1‑150 (MQ=255)
                                                                                                                               ggTCGATGCCAATACTCACCAGTAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAgg                        >  1:1215441/1‑150 (MQ=255)
                                                                                                                               ggTCGATGCCAATACTCACCAGTAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAgg                        >  1:2807567/1‑150 (MQ=255)
                                                                                                                               ggTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAgg                        >  2:4682723/1‑150 (MQ=255)
                                                                                                                                 tCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGaa                      <  1:3983392/150‑1 (MQ=255)
                                                                                                                                 tCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGaa                      <  1:1815688/150‑1 (MQ=255)
                                                                                                                                 tCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGaa                      <  2:4722011/150‑1 (MQ=255)
                                                                                                                                 tCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGaa                      >  1:4524586/1‑150 (MQ=255)
                                                                                                                                  cGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTGGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTGAAGAAGGAAt                     >  2:3954215/1‑150 (MQ=255)
                                                                                                                                   gATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAAtt                    >  1:550798/1‑150 (MQ=255)
                                                                                                                                    aTGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTg                   <  2:2400856/150‑1 (MQ=255)
                                                                                                                                    aTGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTg                   <  1:4264318/150‑1 (MQ=255)
                                                                                                                                        cAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACa               <  2:4251049/150‑1 (MQ=255)
                                                                                                                                          aTACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATg             <  1:3072690/150‑1 (MQ=255)
                                                                                                                                              tCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCAt         >  2:2877537/1‑150 (MQ=255)
                                                                                                                                                  cAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCATTAcc     >  2:1141240/1‑150 (MQ=255)
                                                                                                                                                   aGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCg    >  2:2664891/1‑150 (MQ=255)
                                                                                                                                                    ggAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCGTAAGAAGGAAGTGAACATGTCATTACCGa   >  1:1330352/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGAT  >  CP009273/4411718‑4412016

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: