Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,869,881 T→A 18.4% K144N (AAA→AAT yidA ← sugar phosphate phosphatase; substrates include erythrose 4‑P and mannose 1‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,869,8810TA18.4% 811.0 / 171.5 305K144N (AAA→AATyidAsugar phosphate phosphatase; substrates include erythrose 4‑P and mannose 1‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base T (140/109);  new base A (26/30);  total (166/139)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 311 total aligned reads displayed.

CGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAA  >  CP009273/3869732‑3870030
                                                                                                                                                     |                                                                                                                                                     
cgTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCAtt                                                                                                                                                       <  2:2890409/150‑1 (MQ=255)
 gTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCAatt                                                                                                                                                      >  2:1320208/1‑150 (MQ=255)
   ttATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATtc                                                                                                                                                    >  1:284625/1‑150 (MQ=255)
    tATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTtct                                                                                                                                                   >  1:991686/1‑150 (MQ=255)
    tATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTtct                                                                                                                                                   >  2:4178203/1‑150 (MQ=255)
     aTCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTtctc                                                                                                                                                  >  2:4827782/1‑150 (MQ=255)
               tCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAg                                                                                                                                        <  1:2016306/150‑1 (MQ=255)
                cGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGa                                                                                                                                       >  1:1896425/1‑150 (MQ=255)
                  aGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAac                                                                                                                                     <  2:1271669/150‑1 (MQ=255)
                  aGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAac                                                                                                                                     <  2:4398701/150‑1 (MQ=255)
                     aaGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACAcc                                                                                                                                  >  2:3698927/1‑150 (MQ=255)
                        tACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGc                                                                                                                               <  2:3338814/150‑1 (MQ=255)
                         aCGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCg                                                                                                                              >  2:4092020/1‑150 (MQ=255)
                           ggCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCACCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGa                                                                                                                            >  2:316977/1‑150 (MQ=255)
                            gcgcACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGaa                                                                                                                           <  1:1895053/150‑1 (MQ=255)
                            gcgcACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGaa                                                                                                                           >  2:1835860/1‑150 (MQ=255)
                             cgcACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAAt                                                                                                                          <  2:1935028/150‑1 (MQ=255)
                                aCTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAAtggt                                                                                                                       >  2:1445664/1‑150 (MQ=255)
                                aCTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAAtggt                                                                                                                       <  2:2981227/150‑1 (MQ=255)
                                   tttCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGc                                                                                                                    >  2:3663859/1‑150 (MQ=255)
                                   tttCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGc                                                                                                                    >  2:2514061/1‑150 (MQ=255)
                                   tttCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGc                                                                                                                    <  1:4583858/150‑1 (MQ=255)
                                    ttCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCa                                                                                                                   <  1:4416151/150‑1 (MQ=255)
                                           cGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAgg                                                                                                            >  1:2231838/1‑150 (MQ=255)
                                             gTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGAt                                                                                                          >  2:4012810/1‑150 (MQ=255)
                                              tataTTTCTCTTTCACTTCCTGCGGAATACGCGCGTTTGCCTGGTCGAGGTTGGCTGGTTCATCCATCATCATCTCTTTGAGGAACTGGGTATTGGGGTCCTTTTTCTCCGCTTCGCAGATCACCAGCGGAATGGTGGCAACGAAGGAtt                                                                                                         <  2:2458197/150‑1 (MQ=255)
                                              tataTTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGAtt                                                                                                         <  1:1707487/150‑1 (MQ=255)
                                                tatTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCa                                                                                                       <  2:2557034/150‑1 (MQ=255)
                                                tatTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGAGTCa                                                                                                       >  2:3481869/1‑150 (MQ=255)
                                                 atTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCAt                                                                                                      <  1:4343688/150‑1 (MQ=255)
                                                 atTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCAt                                                                                                      <  1:1958329/150‑1 (MQ=255)
                                                 atTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCAt                                                                                                      >  2:4593642/1‑150 (MQ=255)
                                                  tttCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATg                                                                                                     >  2:3268072/1‑150 (MQ=255)
                                                  tttCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCGTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATg                                                                                                     <  2:2163769/150‑1 (MQ=255)
                                                    tctcTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCa                                                                                                   <  2:746605/150‑1 (MQ=255)
                                                      tctTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCAcc                                                                                                 >  1:3305104/1‑150 (MQ=255)
                                                       ctTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCg                                                                                                <  1:1356612/150‑1 (MQ=255)
                                                       ctTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCg                                                                                                >  1:3018406/1‑150 (MQ=255)
                                                        tttCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCgt                                                                                               >  1:4282228/1‑150 (MQ=255)
                                                           cACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTgta                                                                                            >  1:2374612/1‑150 (MQ=255)
                                                            acgtCCAGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTgtag                                                                                           <  2:4176376/147‑1 (MQ=255)
                                                            aCTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTgtag                                                                                           >  1:1309879/1‑150 (MQ=255)
                                                              ttCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTgtagta                                                                                         <  1:962912/150‑1 (MQ=255)
                                                               tCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGtagtag                                                                                        >  2:1474228/1‑150 (MQ=255)
                                                                ccTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGc                                                                                       >  1:1477805/1‑150 (MQ=255)
                                                                 cTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCt                                                                                      >  2:2263467/1‑150 (MQ=255)
                                                                       aaTACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATc                                                                                >  1:3806684/1‑150 (MQ=255)
                                                                             cgcgATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGtt                                                                          <  2:1888915/150‑1 (MQ=255)
                                                                                      cTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTa                                                                 <  2:1371464/150‑1 (MQ=255)
                                                                                        ggTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACa                                                               >  1:3665526/1‑150 (MQ=255)
                                                                                         gTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAg                                                              <  2:161589/150‑1 (MQ=255)
                                                                                           cGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCg                                                            >  1:305393/1‑150 (MQ=255)
                                                                                             aGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCgtg                                                          <  2:622251/150‑1 (MQ=255)
                                                                                                aTGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCgtggtg                                                       <  1:4091294/150‑1 (MQ=255)
                                                                                                   gCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCgg                                                    >  2:297943/1‑150 (MQ=255)
                                                                                                    cGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGATTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGt                                                   >  2:1750875/1‑150 (MQ=255)
                                                                                                      ggTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTcc                                                 >  1:835704/1‑150 (MQ=255)
                                                                                                      ggTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTcc                                                 <  2:424809/150‑1 (MQ=255)
                                                                                                         tcatcaATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAgg                                              <  1:4123093/150‑1 (MQ=255)
                                                                                                           atcaatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGc                                            >  1:826605/1‑150 (MQ=255)
                                                                                                           atcaatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGc                                            >  1:2450070/1‑150 (MQ=255)
                                                                                                              aatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTg                                         >  2:3380317/1‑150 (MQ=255)
                                                                                                              aatcaTCATCACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTg                                         >  1:2920402/1‑150 (MQ=255)
                                                                                                               atcatcatcACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTgg                                        <  1:2219604/150‑1 (MQ=255)
                                                                                                               atcatcatcACTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTgg                                        <  1:2868953/150‑1 (MQ=255)
                                                                                                                tcatcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGa                                       >  1:4689121/1‑150 (MQ=255)
                                                                                                                   tcatcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAAt                                    <  2:4643452/150‑1 (MQ=255)
                                                                                                                    catcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATg                                   <  2:1747589/150‑1 (MQ=255)
                                                                                                                    catcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATg                                   <  1:1630539/150‑1 (MQ=255)
                                                                                                                     atcaCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATga                                  >  1:1005506/1‑150 (MQ=255)
                                                                                                                       caCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATgaga                                >  1:1096082/1‑150 (MQ=255)
                                                                                                                       caCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATgaga                                >  1:1772498/1‑150 (MQ=255)
                                                                                                                       caCTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGATGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATgaga                                >  1:1486748/1‑150 (MQ=255)
                                                                                                                       caCTTTCAGGAACTGGGTATTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATgaga                                <  2:1621775/150‑1 (MQ=255)
                                                                                                                          tttCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAAcc                             <  1:1211796/150‑1 (MQ=255)
                                                                                                                           ttCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCg                            >  2:3833182/1‑150 (MQ=255)
                                                                                                                              aGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACt                         >  1:3347766/1‑150 (MQ=255)
                                                                                                                              aGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACt                         >  1:2076591/1‑150 (MQ=255)
                                                                                                                              aGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACt                         <  1:2288064/150‑1 (MQ=255)
                                                                                                                              aGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGAAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACt                         >  1:3945828/1‑150 (MQ=255)
                                                                                                                               ggAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACtt                        <  1:4499749/150‑1 (MQ=255)
                                                                                                                               ggAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACtt                        >  1:4853869/1‑150 (MQ=255)
                                                                                                                               ggAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACtt                        <  1:22658/150‑1 (MQ=255)
                                                                                                                               ggAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACtt                        <  1:4807228/150‑1 (MQ=255)
                                                                                                                                gAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTc                       <  2:1638148/150‑1 (MQ=255)
                                                                                                                                   cTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTcgcg                    >  1:1412343/1‑150 (MQ=255)
                                                                                                                                     gggTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCgag                  <  1:1751152/150‑1 (MQ=255)
                                                                                                                                     gggTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCgag                  <  1:2634823/150‑1 (MQ=255)
                                                                                                                                        tATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCgagaga               >  2:1066097/1‑150 (MQ=255)
                                                                                                                                         aTTGGGGTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGagagag              >  2:4716021/1‑150 (MQ=255)
                                                                                                                                          ttGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGt             <  1:2465057/150‑1 (MQ=255)
                                                                                                                                            ggggTCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGttt           <  1:3132359/150‑1 (MQ=255)
                                                                                                                                             gggTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGtttt          <  1:2677899/150‑1 (MQ=255)
                                                                                                                                             gggTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGtttt          >  2:3037301/1‑150 (MQ=255)
                                                                                                                                                tCCATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTcc       >  2:263980/1‑150 (MQ=255)
                                                                                                                                                 ccATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCa      >  2:4506695/1‑150 (MQ=255)
                                                                                                                                                 ccATATTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCa      >  1:3249790/1‑150 (MQ=255)
                                                                                                                                                  cATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAg     <  2:148996/150‑1 (MQ=255)
                                                                                                                                                   aTTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAgg    <  1:3396961/150‑1 (MQ=255)
                                                                                                                                                    ttttCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGa   >  1:4746732/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATACGCGCGATAGCCTGGTCGAGGATGGCGGGTTCATCAATCATCATCACTTTCAGGAACTGGGTATTGGGGTCCATTTTCTCCGCTTCGCAGAACACCAGCGGAATGGTGGCAACGAAGGATTCATGCACCGTGTAGTAGCTGATATCACGGTTGGCGGTGTACAGCGTGGTGCGGTCCAGGGCGTGGAAATGAGAACCGACTTCGCGAGAGAGTTTTTCCAGGAA  >  CP009273/3869732‑3870030

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: