Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,423,366 C→A 10.0% E118* (GAA→TAA)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,423,3660CA10.0% 654.7 / 2.6 242E118* (GAA→TAA) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (59/136);  new base A (31/3);  total (103/139)
Fisher's exact test for biased strand distribution p-value = 1.35e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

CGCGCCGCGGATTAAGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCA  >  CP009273/2423217‑2423515
                                                                                                                                                     |                                                                                                                                                     
cgcgCCGCGGATTAAGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTttc                                                                                                                                                       >  1:1455630/1‑150 (MQ=255)
  cgccgcGGATGAAGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTttctt                                                                                                                                                     <  2:2864895/150‑1 (MQ=255)
        ggATTAAGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACa                                                                                                                                               >  2:1523065/1‑150 (MQ=255)
        ggATTAAGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACa                                                                                                                                               <  2:5174268/150‑1 (MQ=255)
             aaGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGTCGTTTTTCTTCAAACAGtttt                                                                                                                                          <  2:1686177/150‑1 (MQ=255)
             aaGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTTTTCAAAAATtttt                                                                                                                                          >  1:514064/1‑150 (MQ=255)
             aaGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAAAAGtttt                                                                                                                                          >  2:553028/1‑150 (MQ=255)
                 ggTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTaa                                                                                                                                      >  2:5243603/1‑149 (MQ=255)
                 ggTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTGTTTATTAAAAAAATTTTTTaa                                                                                                                                      >  1:5028430/1‑149 (MQ=255)
                    ttGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCa                                                                                                                                   <  1:1027118/150‑1 (MQ=255)
                      gtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAAc                                                                                                                                 <  1:4617489/150‑1 (MQ=255)
                       tgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACt                                                                                                                                <  2:2706384/150‑1 (MQ=255)
                          cagcaATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCCAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGt                                                                                                                             >  2:4714340/1‑150 (MQ=255)
                              aaTGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGc                                                                                                                         <  2:3201978/150‑1 (MQ=255)
                                 ggCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGtt                                                                                                                      <  1:4997079/150‑1 (MQ=255)
                                 ggCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGtt                                                                                                                      <  1:830388/150‑1 (MQ=255)
                                 ggCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTATTAAAAAAGTTTTTTAAGCAACTCGTGAGCGtt                                                                                                                      >  1:4656546/1‑150 (MQ=255)
                                 ggCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTGTTTGTTAAAAAAATTTTTTAAGCAAAACGTGAGCGtt                                                                                                                      >  2:2651662/1‑150 (MQ=255)
                                  gCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTg                                                                                                                     <  1:3538617/150‑1 (MQ=255)
                                    aGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGt                                                                                                                   >  2:1201211/1‑150 (MQ=255)
                                      cGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCa                                                                                                                 <  2:4510582/150‑1 (MQ=255)
                                      cGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCa                                                                                                                 <  1:684511/150‑1 (MQ=255)
                                      cGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCa                                                                                                                 <  2:4474502/150‑1 (MQ=255)
                                      cGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCa                                                                                                                 <  1:2517286/150‑1 (MQ=255)
                                           atatTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTg                                                                                                            <  2:3598574/150‑1 (MQ=255)
                                            tatTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTAGTGAGCGTTGGTCAGATTGc                                                                                                           >  2:4812254/1‑150 (MQ=255)
                                            tatTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTATTCAAACAATTTTTTACGCAAAACGTGAGCGTTGGTCAGAtttc                                                                                                           >  1:418472/1‑148 (MQ=255)
                                             atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGcc                                                                                                          <  2:3868463/150‑1 (MQ=255)
                                              ttGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCg                                                                                                         <  2:1057364/150‑1 (MQ=255)
                                              ttGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCg                                                                                                         <  1:4720711/150‑1 (MQ=255)
                                               tGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTAAAAAAATTTTTTAAACAAAAAGTTATAGGTGGTGAGATTGCCGt                                                                                                        >  1:748698/1‑150 (MQ=255)
                                                 tCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTATTCAAACAGTTTTTTACGCAACACGTGAGCGTTGGTCAGATTGCCGGtg                                                                                                      >  2:3678697/1‑150 (MQ=255)
                                                      ttCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATATCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTATTAAAACAGTTTTTTAAGCAACACGTGAGCGTTGGTCAGATTTCCGGTGTGGGc                                                                                                 >  1:424374/1‑150 (MQ=255)
                                                        ccAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTATTAAAAAATTTTTTTACGCCAAAAGTGAGCGGTGGTCAGATTGCCCTTGTGGGcca                                                                                               >  1:1556879/1‑148 (MQ=255)
                                                         cAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAg                                                                                              <  2:4216293/150‑1 (MQ=255)
                                                          aGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACACGTGAGCGTTGGTCAGATTTCCGTTTTGGGCAAGc                                                                                             >  1:2970004/1‑150 (MQ=255)
                                                           gCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTTCCGTTGTGGGCAAGCg                                                                                            >  2:4284198/1‑150 (MQ=255)
                                                           gCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCg                                                                                            <  1:3193099/150‑1 (MQ=255)
                                                            cGGGTAGTGGAGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTTTTAAAAAAATTTTTTTAGCAAAACGAGAGCGTTGGGCAGATATCAGTTGTGGGCAAGCGt                                                                                           >  1:2005303/1‑150 (MQ=255)
                                                             gggTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTa                                                                                          >  2:2930215/1‑150 (MQ=255)
                                                             gggTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTa                                                                                          <  1:2199383/150‑1 (MQ=255)
                                                              ggTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTaa                                                                                         <  2:2325391/150‑1 (MQ=255)
                                                               gTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAAt                                                                                        <  2:3060407/150‑1 (MQ=255)
                                                                tAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAAAAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAAtt                                                                                       >  2:1974696/1‑149 (MQ=255)
                                                                  gTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGcc                                                                                     <  2:2753675/150‑1 (MQ=255)
                                                                   tGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCa                                                                                    <  1:4304986/150‑1 (MQ=255)
                                                                      cGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATAc                                                                                 <  1:1327194/150‑1 (MQ=255)
                                                                       ggAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACg                                                                                <  1:3640038/150‑1 (MQ=255)
                                                                       ggAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGGTTTTTACGCAACTCGTGAGCGGTGGTCAGATTGCCGGTGTGGGCAAGCGTTATGCCCTACg                                                                                >  2:692128/1‑150 (MQ=255)
                                                                         aaGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACggg                                                                              <  2:4427959/150‑1 (MQ=255)
                                                                            ttGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAg                                                                           <  1:481484/150‑1 (MQ=255)
                                                                              gTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGtt                                                                         <  1:1687025/150‑1 (MQ=255)
                                                                               tCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTa                                                                        <  2:1455630/150‑1 (MQ=255)
                                                                               tCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTa                                                                        <  2:3782436/150‑1 (MQ=255)
                                                                                   gCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGt                                                                    <  2:3210971/150‑1 (MQ=255)
                                                                                   gCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGt                                                                    <  2:1740383/150‑1 (MQ=255)
                                                                                    cTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTa                                                                   <  1:2804985/150‑1 (MQ=255)
                                                                                     tCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCCTCAAACAGTTTTTTACGCAACACGTGAGCGTTGGTCAGAATGCCGTTGTGGGCAAGCGTAATGCCATACGGGGCGTTAACGGaa                                                                  >  2:112414/1‑150 (MQ=255)
                                                                                         tGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTTCGCAACTCGTGAGCGTTGGTAAGATTTCCGTTGTGGGCAAGCGTAATTCCATACGGGGAGTTAACGTAAAACg                                                              >  1:2620603/1‑150 (MQ=255)
                                                                                         tGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAAAAGTTTTTTAAGCAACTACTGAGCGTTGGTCACATTTCCCGTGTGGGCAAGCGTAATTCCCAACGGGGAGTTAACGTAAAACg                                                              >  2:1931552/1‑150 (MQ=255)
                                                                                         gggCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACGCGAGAGCGTTGGTCCGATTGCCGGGGTGGGCCAGCGTAATGCCAGACCGGGCGTGAACGTAAAACg                                                              >  1:3901464/2‑150 (MQ=255)
                                                                                            cGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCt                                                           <  1:1767427/150‑1 (MQ=255)
                                                                                            cGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAATTTTTTAAGCAACACGTGAGCGTTGGTCAGAATTCCGTTGTGGGCCAGCGTAATGCCAAACGGGGAGTTAACCTAAAAAGGcg                                                           >  1:856234/1‑149 (MQ=255)
                                                                                             gAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTg                                                          <  2:1964147/150‑1 (MQ=255)
                                                                                                 atatTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTATTAAAAAAGTTTTTTACGCAACTCGTGATCGTTGGTCAGATTTCCGGTGTGGGCCAGCGGAATGCCATAAGGGGGGGTTACGTAAAAAGGCTGCGCt                                                      >  2:1376698/1‑150 (MQ=255)
                                                                                                 atatTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTGTTTGTTCAAACAGTTTTTTTAGCAAAAAGGGAGCGTTTGTCAGATTGCCGTTTTGGGCAAGCGGAATGCCAAAACGGGAGTTAACGGAAAAACGCTGCGCt                                                      >  2:4015754/1‑150 (MQ=255)
                                                                                                   atTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTAAAACATTTTTTTAAGCCACACGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATTCCATACGGGGAGGTAACGTAAAAACGCTGCGCTTc                                                    >  2:4391370/1‑150 (MQ=255)
                                                                                                    ttAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAAAAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCa                                                   >  2:4997310/1‑150 (MQ=255)
                                                                                                      aaGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGGTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCaga                                                 >  2:3411739/1‑150 (MQ=255)
                                                                                                         cAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGgc                                              <  2:831120/150‑1 (MQ=255)
                                                                                                            aaTTTCCGTGTCGTATGTGGTTTTGTTGTGGCGGCTTTTTTCTTCAAACTTTTTTTTACGAAACTCGTGAGCGTTGGTCAGTTTGCGGTTGTGGGCAAGCTTAATGCAATAGGGGGTGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCt                                           <  1:906292/150‑1 (MQ=255)
                                                                                                             aTTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTGTTCAAACAGTTTTTTACGCAACTAGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCAAACGGGGAGTTAACGTAAAACGGCTGCGCTTAAGAGGCGCTg                                          >  2:2897619/1‑150 (MQ=255)
                                                                                                               ttCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGa                                        <  1:1189213/150‑1 (MQ=255)
                                                                                                                  cGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCt                                     <  2:4139786/150‑1 (MQ=255)
                                                                                                                      tcgtatgtgttgTGGTTGTGGCGGCGTTTTCCTTAAAACAGTTTTTTACGCAACTCGTGTGCGTTGGTCAGATTGCCGTTGTGGGCAAGCTTAATGCCATAGGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCa                                 <  1:4417687/144‑1 (MQ=255)
                                                                                                                        ggATGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCATCTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGc                               <  1:2030538/150‑1 (MQ=255)
                                                                                                                        ggAAGTGGTGTTGATGTGGCGGCGTTTTTCATCAAACAATTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATTCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGTGCTGCCAGc                               >  2:674698/1‑150 (MQ=255)
                                                                                                                        ggAAGTGGTGTTGATGTGGCGGCGTTTTTATTCAAACAGGTTTTTTAGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATTCCATACGGGGTGTTATCGTAAAACGGCTGCGCGTCAGAGGCGCAGGAGCTGCCAGc                               >  2:1063822/1‑150 (MQ=255)
                                                                                                                         gAAGTGGTGTTGATGTGGCGGCGTTTTTATTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTTCCGTTGTGGGCAAGCGTACTGCCATAAGGGGATTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGcc                              >  2:3892162/1‑150 (MQ=255)
                                                                                                                          aaGTGGTTTTGATGTGGCGGCGTTTTTATTCAAACAGTTTTTAACGCAAATCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCTTCCGGGGATTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCg                             <  1:5089311/150‑1 (MQ=255)
                                                                                                                          aaGTGGTGTTGATGTGGCGGCGTTTTTATTCAAAAATTTTTTTTAGCAACACGTGAGCGTTTGTAAGATTGCCGGTGTGGGCAAACCTAATTCCATACCGGGAGTTAAAATAAAACCGCAGCGCGTAAAAGGCGCAGGAGCTGCCAcccc                             >  1:380643/1‑149 (MQ=255)
                                                                                                                          aaGTGGTGTTGAGGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGAGTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCg                             <  2:4964268/150‑1 (MQ=255)
                                                                                                                            gtggtgTTGATGTGGCGGCGTGTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGGTGTGGGCCAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTg                           >  2:3375188/1‑150 (MQ=255)
                                                                                                                                 tttgttgTGGCGTCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCACGCGTAATGCAATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTa                      <  2:4375761/149‑1 (MQ=255)
                                                                                                                                   tGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAAc                    <  2:248116/150‑1 (MQ=255)
                                                                                                                                   tGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAAc                    <  1:3478647/150‑1 (MQ=255)
                                                                                                                                    gATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATCCGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACg                   <  2:4549358/150‑1 (MQ=255)
                                                                                                                                    gATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACg                   <  1:4922895/150‑1 (MQ=255)
                                                                                                                                    gATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACg                   <  1:3404334/150‑1 (MQ=255)
                                                                                                                                           cggcgTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCCAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGa            >  2:3457370/1‑150 (MQ=255)
                                                                                                                                            ggcgTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGAc           <  1:4056566/150‑1 (MQ=255)
                                                                                                                                            ggcgTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGAc           >  1:371197/1‑150 (MQ=255)
                                                                                                                                            ggcgTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGAc           >  2:1902471/1‑150 (MQ=255)
                                                                                                                                             gcgTTTTTCTTCAAACAGTTTTTTCCGCAACTCGTGAGCGTTGGTCAGTTTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGAcc          <  2:2213601/150‑1 (MQ=255)
                                                                                                                                               ttttttATTAAAAAAGTTTTTTACGCAAGTCGTGTGGTTTGGTCATTTTGCGGTTGTGGGCAAGCTTAATGCCATAGGGGGTGTTAACTTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCaa        <  1:2651662/149‑1 (MQ=255)
                                                                                                                                                  tttCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGc     >  1:4813758/1‑150 (MQ=255)
                                                                                                                                                   ttcttcAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGcc    >  2:1911629/1‑150 (MQ=255)
                                                                                                                                                    tcttcAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGccc   <  1:246628/150‑1 (MQ=255)
                                                                                                                                                     cttcAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCa  >  1:454383/1‑150 (MQ=255)
                                                                                                                                                     cttcAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCa  >  1:534229/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCCGCGGATTAAGCGGTTTGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTTTTTCTTCAAACAGTTTTTTACGCAACTCGTGAGCGTTGGTCAGATTGCCGTTGTGGGCAAGCGTAATGCCATACGGGGAGTTAACGTAAAACGGCTGCGCTTCAGAGGCGCTGGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCA  >  CP009273/2423217‑2423515

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑