Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,110,346 G→A 100% T88I (ACC→ATC)  wcaL ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,110,3460GA95.2% 524.8 / 13.3 167T88I (ACC→ATC) wcaLputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (6/2);  new base A (74/85);  total (80/87)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 167 total aligned reads displayed.

GAGTTTTGCTGCGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGC  >  CP009273/2110201‑2110490
                                                                                                                                                 |                                                                                                                                                
   ttttGCTGCGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATAttt                                                                                                                                             >  1:2338458/1‑148 (MQ=255)
   ttttGCTGCGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCAGAGCGTTTGAGGTTGAGCGCCTGCCAGATAttt                                                                                                                                             >  2:2682390/1‑148 (MQ=255)
     ttGCTGCGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATAttttt                                                                                                                                           <  2:1938247/148‑1 (MQ=255)
         tgcGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGa                                                                                                                                       >  1:798334/1‑148 (MQ=255)
         tgcGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGa                                                                                                                                       >  1:798121/1‑148 (MQ=255)
           cGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATg                                                                                                                                     >  2:617454/1‑148 (MQ=255)
                          gACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGt                                                                                                                      >  1:210535/1‑148 (MQ=255)
                           aCCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTc                                                                                                                     >  2:245323/1‑148 (MQ=255)
                           aCCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTc                                                                                                                     >  2:2233505/1‑148 (MQ=255)
                           aCCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTc                                                                                                                     >  1:1012949/1‑148 (MQ=255)
                           aCCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGGTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTc                                                                                                                     <  1:2500333/148‑1 (MQ=255)
                            ccAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTct                                                                                                                    >  2:172243/1‑148 (MQ=255)
                             cAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTctg                                                                                                                   <  1:1645961/148‑1 (MQ=255)
                              aaaaTGAGCGATGAACACATCGGCGCGAAACGGTGTTTCGACCTGGCCGAAAATGGCAGACAAAATCAGTTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTctgg                                                                                                                  <  1:2492673/148‑1 (MQ=255)
                                aaTGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTctggct                                                                                                                >  1:1908000/1‑148 (MQ=255)
                                  tGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTctggctgg                                                                                                              >  2:2220640/1‑148 (MQ=255)
                                   gAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCtggctggc                                                                                                             >  2:389825/1‑148 (MQ=255)
                                    aGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTggctggct                                                                                                            >  2:832687/1‑148 (MQ=255)
                                        aTGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGt                                                                                                        <  2:2266845/148‑1 (MQ=255)
                                         tGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTg                                                                                                       <  2:462273/148‑1 (MQ=255)
                                          gAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTgg                                                                                                      >  1:1065131/1‑148 (MQ=255)
                                           aaCACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGgc                                                                                                     >  2:1088664/1‑148 (MQ=255)
                                           aaCACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGgc                                                                                                     >  1:95771/1‑148 (MQ=255)
                                             cacaTCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGgcgc                                                                                                   <  2:1489561/148‑1 (MQ=255)
                                               caTCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTGTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAg                                                                                                 >  2:3231490/1‑148 (MQ=255)
                                                aTCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGt                                                                                                <  1:1825857/148‑1 (MQ=255)
                                                 tCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGtt                                                                                               >  2:3352384/1‑148 (MQ=255)
                                                  cGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGttt                                                                                              >  1:1085580/1‑148 (MQ=255)
                                                        gAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCAc                                                                                        <  1:1488317/148‑1 (MQ=255)
                                                         aatcGGTGTTGCGACCTGGCCTCAAATGGCTTACAATATGAGGTTCCGCGACTCGGGAACATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGATTGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGcctgt                                                                                       <  2:3355098/145‑4 (MQ=255)
                                                          aaCGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACtt                                                                                      <  1:1973641/148‑1 (MQ=255)
                                                             ggTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGc                                                                                   <  2:2106492/148‑1 (MQ=255)
                                                                  tGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTa                                                                              <  2:1362404/148‑1 (MQ=255)
                                                                  tGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTa                                                                              >  1:1621978/1‑148 (MQ=255)
                                                                     gACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGt                                                                           <  1:617939/148‑1 (MQ=255)
                                                                     gACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGt                                                                           <  2:278890/148‑1 (MQ=255)
                                                                      aCCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGtt                                                                          >  1:45600/1‑148 (MQ=255)
                                                                        cTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCg                                                                        <  2:1157623/148‑1 (MQ=255)
                                                                         tGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGt                                                                       <  1:1245933/148‑1 (MQ=255)
                                                                          ggCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTc                                                                      <  2:2528810/148‑1 (MQ=255)
                                                                          ggCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTc                                                                      <  1:1542931/148‑1 (MQ=255)
                                                                          ggCCGCAAATGGCAGACAAAATCAGGTACCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTc                                                                      <  1:535958/148‑1 (MQ=255)
                                                                           gCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTcc                                                                     >  2:2308694/1‑148 (MQ=255)
                                                                           gCCGCAAAACCCACACAAAATCACCTTCCGCGACTCGGCACCATAGCGATTGAGGATGAGCGCCTGCCAGAAAATTTTACGATGAAAGCCGCGCAAGGTCTGGCTGGCTCGCTGGCGCAGTTTCGCCACTTTGCCCGAAGGTtcctcc                                                                     <  2:3503912/148‑1 (MQ=255)
                                                                             cGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTg                                                                   >  1:1431561/1‑148 (MQ=255)
                                                                               cAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGTCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTa                                                                 <  2:2057264/148‑1 (MQ=255)
                                                                                aaaTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTaa                                                                >  1:2522055/1‑148 (MQ=255)
                                                                                 aaTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAc                                                               <  1:1204666/148‑1 (MQ=255)
                                                                                 aaTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCAGAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAc                                                               <  1:645046/148‑1 (MQ=255)
                                                                                  aTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAcc                                                              >  2:1458131/1‑148 (MQ=255)
                                                                                  aTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAcc                                                              <  1:3440186/148‑1 (MQ=255)
                                                                                  aTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAcc                                                              >  2:2629852/1‑148 (MQ=255)
                                                                                  aTGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAcc                                                              <  2:2574729/148‑1 (MQ=255)
                                                                                  aTGGCAGACAAAATCAGGGTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAcc                                                              <  2:788235/148‑1 (MQ=255)
                                                                                   tGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCa                                                             >  2:1059336/1‑148 (MQ=255)
                                                                                    ggCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAg                                                            <  1:1595646/148‑1 (MQ=255)
                                                                                    ggCAGACAAAATCAGGGTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAg                                                            <  1:2545304/148‑1 (MQ=255)
                                                                                      cAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCg                                                          >  1:2592187/1‑148 (MQ=255)
                                                                                        gACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCggg                                                        >  2:1302724/1‑148 (MQ=255)
                                                                                        gACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGATCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCggg                                                        <  2:1899911/148‑1 (MQ=255)
                                                                                           aaaaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTc                                                     <  1:3373532/148‑1 (MQ=255)
                                                                                            aaaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCt                                                    <  2:3107413/148‑1 (MQ=255)
                                                                                             aaTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTg                                                   <  2:53934/148‑1 (MQ=255)
                                                                                              aTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTgg                                                  <  1:2233505/148‑1 (MQ=255)
                                                                                              aTCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTgg                                                  <  1:245323/148‑1 (MQ=255)
                                                                                               tCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgc                                                 >  2:838352/1‑148 (MQ=255)
                                                                                                cAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgca                                                >  2:230091/1‑148 (MQ=255)
                                                                                                  ggTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgcagc                                              >  2:1703352/1‑148 (MQ=255)
                                                                                                   gTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgcagca                                             <  2:2248851/148‑1 (MQ=255)
                                                                                                     tCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAg                                           <  1:3321778/148‑1 (MQ=255)
                                                                                                     tCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAg                                           <  1:3321744/148‑1 (MQ=255)
                                                                                                     tCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGAGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAg                                           >  2:3443125/1‑148 (MQ=255)
                                                                                                         cgACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTg                                       <  1:3274707/148‑1 (MQ=255)
                                                                                                         cgACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTg                                       >  2:2988268/1‑148 (MQ=255)
                                                                                                           aCTCGGCACCATAGCGTTTGAGGTTGAGCGACTGCCAGATATTTTTACGATGAATTACGCTCAATGTAGGTCTTGATCGGTTTCGAATTTTCGACACTTTGCCCGTAGTTTCGTCCTTTAACCAGCGGGTTCTGGCAGCAAGTTTGTa                                     <  2:3188477/148‑1 (MQ=255)
                                                                                                            cTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTAt                                    <  2:2924369/148‑1 (MQ=255)
                                                                                                                 cACCATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGt                               >  1:3395250/1‑148 (MQ=255)
                                                                                                                    cATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCa                            >  2:371756/1‑148 (MQ=255)
                                                                                                                    cATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCa                            <  2:1345320/148‑1 (MQ=255)
                                                                                                                    cATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCa                            <  2:2860864/148‑1 (MQ=255)
                                                                                                                    cATAGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCa                            >  2:2172651/1‑148 (MQ=255)
                                                                                                                      tAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATg                          >  2:1620445/1‑148 (MQ=255)
                                                                                                                       aGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCTTCAATGc                         >  2:2263756/1‑148 (MQ=255)
                                                                                                                       aGCGTTTGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGc                         >  2:468700/1‑148 (MQ=255)
                                                                                                                           tttGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCcgcg                     <  1:2815920/148‑1 (MQ=255)
                                                                                                                            ttGAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGt                    <  2:2026351/148‑1 (MQ=255)
                                                                                                                              gAGGTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTgg                  <  2:2592187/148‑1 (MQ=255)
                                                                                                                                ggTTGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGgt                >  2:197104/1‑148 (MQ=255)
                                                                                                                                  ttGAGCGCCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGgtgt              >  1:1744964/1‑148 (MQ=255)
                                                                                                                                   tGAGCGCCTGCCAGATATTTTTCCGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGtt             <  2:2666733/148‑1 (MQ=255)
                                                                                                                                      gcgcCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtg          <  2:798121/148‑1 (MQ=255)
                                                                                                                                      gcgcCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtg          <  2:798334/148‑1 (MQ=255)
                                                                                                                                       cgcCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtgt         <  2:3361734/148‑1 (MQ=255)
                                                                                                                                       cgcCTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtgt         <  2:1593542/148‑1 (MQ=255)
                                                                                                                                         ccTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtgtgt       <  2:3492124/148‑1 (MQ=255)
                                                                                                                                         ccTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtgtgt       <  2:1766518/148‑1 (MQ=255)
                                                                                                                                         ccTGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTtgtgt       >  1:267520/1‑148 (MQ=255)
                                                                                                                                           tGCCAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTgtgtgt     >  2:323226/1‑148 (MQ=255)
                                                                                                                                             ccAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCg   >  2:1373043/1‑148 (MQ=255)
                                                                                                                                              cAGATATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGc  <  2:1076160/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAGTTTTGCTGCGGTTACCCCCGCGGGACCAAAATGAGCGATGAACACATCGGCGCGAAACGGTGTTGCGACCTGGCCGCAAATGGCAGACAAAATCAGGTTCCGCGACTCGGCACCATAGCGTTTGAGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCATGCCGCGTGGGTGTTTTGTGTGTCGC  >  CP009273/2110201‑2110490

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: