Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,260,920 G→T 6.9% R112L (CGT→CTT)  yjbR → DUF419 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,260,9200GT 6.9% 1111.6 / 39.0 289R112L (CGT→CTT) yjbRDUF419 family protein
Reads supporting (aligned to +/- strand):  ref base G (134/135);  new base T (11/9);  total (145/144)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 291 total aligned reads displayed.

CAGCACAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGCCACC  >  CP009273/4260772‑4261069
                                                                                                                                                    |                                                                                                                                                     
cagcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTCTCAGCAGGCGGTGAATTTACTGCCGGAGGAAAAACGt                                                                                                                                                      >  2:4458907/1‑150 (MQ=255)
cagcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGt                                                                                                                                                      <  2:841991/150‑1 (MQ=255)
 agcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTa                                                                                                                                                     >  2:2671709/1‑150 (MQ=255)
   cacaGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTaaa                                                                                                                                                   <  2:3677213/150‑1 (MQ=255)
     caGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAAtt                                                                                                                                                 >  1:2383089/1‑150 (MQ=255)
       gCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGc                                                                                                                                               <  1:3033356/150‑1 (MQ=255)
         gATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTg                                                                                                                                             <  1:4488852/150‑1 (MQ=255)
         gATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTg                                                                                                                                             >  1:4062798/1‑150 (MQ=255)
         gATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTg                                                                                                                                             >  2:632453/1‑150 (MQ=255)
           tgtgCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGt                                                                                                                                           <  2:2036755/150‑1 (MQ=255)
              gCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCa                                                                                                                                        <  2:2368918/150‑1 (MQ=255)
                  ccAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAActc                                                                                                                                    <  1:3630980/150‑1 (MQ=255)
                  ccAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAActc                                                                                                                                    >  1:1884744/1‑150 (MQ=255)
                   cAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAActct                                                                                                                                   >  2:391770/1‑150 (MQ=255)
                     aGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGa                                                                                                                                 <  1:1804774/150‑1 (MQ=255)
                        cgccATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAg                                                                                                                              >  2:3417046/1‑150 (MQ=255)
                          ccATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGa                                                                                                                            >  1:4407116/1‑150 (MQ=255)
                           cATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaa                                                                                                                           >  2:533800/1‑150 (MQ=255)
                            aTCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaaa                                                                                                                          >  1:3140317/1‑150 (MQ=255)
                             tCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaaaa                                                                                                                         >  1:2476017/1‑150 (MQ=255)
                             tCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaaaa                                                                                                                         <  1:575287/150‑1 (MQ=255)
                                gAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGc                                                                                                                      >  2:1766337/1‑150 (MQ=255)
                                 aaTAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGcc                                                                                                                     <  1:4112197/150‑1 (MQ=255)
                                 aaTAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGcc                                                                                                                     >  1:2242555/1‑150 (MQ=255)
                                  aTAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCg                                                                                                                    <  2:4876175/150‑1 (MQ=255)
                                      aGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAGAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCa                                                                                                                <  2:1361312/150‑1 (MQ=255)
                                       gcgcACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAg                                                                                                               <  2:4590991/150‑1 (MQ=255)
                                         gcACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGaa                                                                                                             >  1:2974050/1‑150 (MQ=255)
                                          cACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGaaa                                                                                                            >  2:1485258/1‑150 (MQ=255)
                                           aCTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAg                                                                                                           <  2:4560371/150‑1 (MQ=255)
                                            cTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGc                                                                                                          <  1:2726430/150‑1 (MQ=255)
                                             tGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCg                                                                                                         >  2:3091376/1‑150 (MQ=255)
                                               gAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGc                                                                                                       <  2:1451193/150‑1 (MQ=255)
                                               gAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGc                                                                                                       <  1:2082218/150‑1 (MQ=255)
                                                    ccGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTaa                                                                                                  <  1:2770154/150‑1 (MQ=255)
                                                     cGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTatc                                                                                                 <  2:125958/150‑3 (MQ=255)
                                                     cGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAAc                                                                                                 <  2:741691/150‑1 (MQ=255)
                                                      gtgtATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACg                                                                                                <  1:1228524/150‑1 (MQ=255)
                                                      gtgtATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACg                                                                                                <  1:2011552/150‑1 (MQ=255)
                                                      gtgtATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACg                                                                                                <  2:2306668/150‑1 (MQ=255)
                                                          aTCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTa                                                                                            >  2:3295711/1‑150 (MQ=255)
                                                          aTCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTa                                                                                            >  2:1941693/1‑150 (MQ=255)
                                                           tctcGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTAc                                                                                           >  2:3525408/1‑150 (MQ=255)
                                                             tcGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCTGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAACTTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAg                                                                                         <  2:3961411/150‑1 (MQ=255)
                                                               gACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGACGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGACAGCGGCCTTAACGATTACAGCa                                                                                       >  2:2081231/1‑150 (MQ=255)
                                                               gACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCa                                                                                       >  2:215207/1‑150 (MQ=255)
                                                               gACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCa                                                                                       <  1:937933/150‑1 (MQ=255)
                                                               gACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCa                                                                                       >  2:2811681/1‑150 (MQ=255)
                                                                aCGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTTCAGCAt                                                                                      <  1:4781150/150‑1 (MQ=255)
                                                                 cGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATc                                                                                     <  2:2824260/150‑1 (MQ=255)
                                                                  ggTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCg                                                                                    <  1:1981977/150‑1 (MQ=255)
                                                                   gTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCgg                                                                                   <  2:2489354/150‑1 (MQ=255)
                                                                    ttCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGc                                                                                  >  2:888121/1‑150 (MQ=255)
                                                                      cGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCtt                                                                                >  2:2033465/1‑150 (MQ=255)
                                                                      cGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCtt                                                                                <  1:141188/150‑1 (MQ=255)
                                                                       gCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTa                                                                               >  2:3235641/1‑150 (MQ=255)
                                                                         tGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAg                                                                             <  1:2618236/150‑1 (MQ=255)
                                                                         tGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAg                                                                             >  1:2199445/1‑150 (MQ=255)
                                                                          gCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAgg                                                                            <  1:1006186/150‑1 (MQ=255)
                                                                           ccAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGa                                                                           >  1:4873139/1‑150 (MQ=255)
                                                                            cAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGaa                                                                          >  2:4530742/1‑150 (MQ=255)
                                                                                 tCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTg                                                                     <  2:1866049/150‑1 (MQ=255)
                                                                                   gCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGcc                                                                   <  1:4646260/150‑1 (MQ=255)
                                                                                     aaaTCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCgt                                                                 >  2:1001063/1‑150 (MQ=255)
                                                                                       aTCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCgtgt                                                               <  1:5037594/150‑1 (MQ=255)
                                                                                       aTCTATTATCGGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCgtgt                                                               <  1:1740334/150‑1 (MQ=255)
                                                                                         cTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCgtgtgt                                                             >  2:4273441/1‑150 (MQ=255)
                                                                                          tattatCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGtgtgtg                                                            <  2:4667941/150‑1 (MQ=255)
                                                                                               tCTGGTGGATGCGACTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGACAGCGGCCTGAACGATTACAGCATCGGCTGAAGGACGCGTGCCGTGTGTGATGCt                                                       >  1:158611/1‑150 (MQ=255)
                                                                                                  ggtggATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCg                                                    >  1:3853500/1‑150 (MQ=255)
                                                                                                   gtggATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGTAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGc                                                   >  1:4244940/1‑150 (MQ=255)
                                                                                                          cGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCcgc                                            >  1:4574129/1‑150 (MQ=255)
                                                                                                          cGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCcgc                                            >  1:5132723/1‑150 (MQ=255)
                                                                                                           gTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCcgcg                                           <  1:3249476/150‑1 (MQ=255)
                                                                                                            tCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGa                                          <  1:3497584/150‑1 (MQ=255)
                                                                                                            tCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGa                                          <  2:4777585/150‑1 (MQ=255)
                                                                                                                aTCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGt                                      <  1:2563261/150‑1 (MQ=255)
                                                                                                                 tCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGtt                                     <  2:1040449/150‑1 (MQ=255)
                                                                                                                    gcagGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCt                                  >  2:2289376/1‑150 (MQ=255)
                                                                                                                     cagGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTg                                 <  2:866953/150‑1 (MQ=255)
                                                                                                                       ggCGGTGAATTTCCTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGc                               <  1:925884/150‑1 (MQ=255)
                                                                                                                       ggCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGc                               <  1:1106126/150‑1 (MQ=255)
                                                                                                                       ggCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGc                               <  1:292540/150‑1 (MQ=255)
                                                                                                                          ggTGAATTTCCTGCCGGACGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATGACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTa                            <  2:1151046/150‑1 (MQ=255)
                                                                                                                          ggTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTa                            >  2:3658609/1‑150 (MQ=255)
                                                                                                                          ggTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCAGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTa                            >  2:1087263/1‑150 (MQ=255)
                                                                                                                           gTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTAc                           >  2:3496524/1‑150 (MQ=255)
                                                                                                                               aTTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGga                       >  1:1296098/1‑150 (MQ=255)
                                                                                                                                 ttACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGgaga                     >  1:4003523/1‑150 (MQ=255)
                                                                                                                                  tACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGAc                    <  2:1562618/150‑1 (MQ=255)
                                                                                                                                      gCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAgg                <  1:4346642/150‑1 (MQ=255)
                                                                                                                                        cGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGAt              >  1:3345050/1‑150 (MQ=255)
                                                                                                                                            agaaAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcg          >  2:4040502/1‑150 (MQ=255)
                                                                                                                                               aaaaCTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcgccg       >  2:4097348/1‑150 (MQ=255)
                                                                                                                                               aaaaCGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcgccg       >  1:3365996/1‑150 (MQ=255)
                                                                                                                                               aaaaCGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcgccg       <  1:2363660/150‑1 (MQ=255)
                                                                                                                                               aaaaCGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcgccg       <  1:1541550/150‑1 (MQ=255)
                                                                                                                                                aaaCGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGAGTGCCGTGTGTGATGCTTCTCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcgccgc      >  2:110077/1‑150 (MQ=255)
                                                                                                                                                 aaCGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCgccgcc     <  2:4596722/150‑1 (MQ=255)
                                                                                                                                                    gTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGccacc  <  2:162397/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCACAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGCCACC  >  CP009273/4260772‑4261069

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: