Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,931,717 C→A 94.9% N307K (AAC→AAA)  fucK → L‑fuculokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,931,7170CA94.9% 862.1 / 35.6 276N307K (AAC→AAA) fucKL‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base C (5/9);  new base A (135/127);  total (140/136)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 278 total aligned reads displayed.

TGGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATT  >  CP009273/2931570‑2931863
                                                                                                                                                   |                                                                                                                                                  
tGGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAacc                                                                                                                                                  <  2:1018412/150‑1 (MQ=255)
tGGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAacc                                                                                                                                                  <  1:2352487/150‑1 (MQ=255)
   cgcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAgg                                                                                                                                               >  1:3147636/1‑150 (MQ=255)
    gcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGt                                                                                                                                              <  2:248612/150‑1 (MQ=255)
     cTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTa                                                                                                                                             >  1:2373183/1‑150 (MQ=255)
      tGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTAt                                                                                                                                            <  2:844500/150‑1 (MQ=255)
      tGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTAt                                                                                                                                            >  2:1599156/1‑150 (MQ=255)
      tGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTAt                                                                                                                                            <  1:1612975/150‑1 (MQ=255)
           cTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAt                                                                                                                                       <  1:1995527/150‑1 (MQ=255)
              aaCAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggc                                                                                                                                    <  1:4166548/150‑1 (MQ=255)
               aCAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTAAGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggct                                                                                                                                   <  2:1837672/150‑1 (MQ=255)
                cAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctg                                                                                                                                  >  1:518253/1‑150 (MQ=255)
                cAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctg                                                                                                                                  <  2:2286339/150‑1 (MQ=255)
                cAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtggctg                                                                                                                                  <  2:1440403/150‑1 (MQ=255)
                   aaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCa                                                                                                                               >  1:3518559/1‑150 (MQ=255)
                   aaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCa                                                                                                                               >  2:1704159/1‑150 (MQ=255)
                    aTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCAt                                                                                                                              <  2:586387/150‑1 (MQ=255)
                    aTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCAt                                                                                                                              <  2:586536/150‑1 (MQ=255)
                    aTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCAt                                                                                                                              >  2:4466354/1‑150 (MQ=255)
                     tGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATc                                                                                                                             <  2:351782/150‑1 (MQ=255)
                      gAACCCGTGCTCACTTCCGGTACATGGGAAATTTTAATGGGACGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCACTGGCTGGCATcc                                                                                                                            >  1:3208446/1‑150 (MQ=255)
                            gTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTg                                                                                                                      >  1:4017582/1‑150 (MQ=255)
                             tGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGc                                                                                                                     >  2:4569329/1‑150 (MQ=255)
                             tGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGc                                                                                                                     <  1:624361/150‑1 (MQ=255)
                                 ctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                                 <  2:4169738/150‑1 (MQ=255)
                                 ctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                                 <  1:3839879/150‑1 (MQ=255)
                                        ggTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTg                                                                                                          >  2:3389317/1‑150 (MQ=255)
                                          tACATGGGAAATTTTAATTGTACGCAGCGACCAGGTTCATACTTCGCTGTTAAGACAGTACGCCGGTTCCACCTGCGAACTGGATATCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCAGGAGTGCTGGAATGGGTgag                                                                                                        >  2:2600418/1‑150 (MQ=255)
                                            cATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGaa                                                                                                      >  2:3032646/1‑150 (MQ=255)
                                               gggAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAAc                                                                                                   >  2:1473878/1‑150 (MQ=255)
                                                       ttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTgg                                                                                           >  2:1281731/1‑150 (MQ=255)
                                                         aaTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGAc                                                                                         <  1:1358649/150‑1 (MQ=255)
                                                           tGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACgg                                                                                       >  2:4641513/1‑150 (MQ=255)
                                                             gTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCt                                                                                     <  1:3512210/150‑1 (MQ=255)
                                                               tCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGa                                                                                   <  2:1752906/150‑1 (MQ=255)
                                                               tCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGa                                                                                   >  2:2177635/1‑150 (MQ=255)
                                                                cGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGaa                                                                                  <  2:2729170/150‑1 (MQ=255)
                                                                 gCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGaaa                                                                                 >  2:1687709/1‑150 (MQ=255)
                                                                 gCAGCGCCCAGGTTGATACTTCGCGGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGaaa                                                                                 >  1:229304/1‑150 (MQ=255)
                                                                   aGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAaca                                                                               >  2:3717502/1‑150 (MQ=255)
                                                                      gcCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACAccc                                                                            <  1:1444307/150‑1 (MQ=255)
                                                                       cccAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCt                                                                           <  1:2804654/150‑1 (MQ=255)
                                                                       cccAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCt                                                                           <  1:4496137/150‑1 (MQ=255)
                                                                        ccAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                          >  2:1183378/1‑150 (MQ=255)
                                                                        ccAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                          >  1:1376783/1‑150 (MQ=255)
                                                                          aGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGc                                                                        >  2:3789965/1‑150 (MQ=255)
                                                                          aGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGc                                                                        >  1:2754177/1‑150 (MQ=255)
                                                                            gTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCACTGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCaa                                                                      >  2:4613011/1‑150 (MQ=255)
                                                                              tGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAAt                                                                    >  2:3557421/1‑150 (MQ=255)
                                                                                  aCTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttg                                                                >  1:4128817/1‑150 (MQ=255)
                                                                                   cTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttga                                                               >  1:3029944/1‑150 (MQ=255)
                                                                                     tCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttgatt                                                             <  2:964224/150‑1 (MQ=255)
                                                                                        cTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaa                                                          >  1:247874/1‑150 (MQ=255)
                                                                                        cTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaa                                                          >  1:860163/1‑150 (MQ=255)
                                                                                         tGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaag                                                         >  1:67275/1‑150 (MQ=255)
                                                                                          gTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaaga                                                        >  2:2031123/1‑150 (MQ=255)
                                                                                            tAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGACATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGaagaag                                                      <  2:3301574/150‑1 (MQ=255)
                                                                                              aGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTAGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCt                                                    <  2:1071923/150‑1 (MQ=255)
                                                                                                 cAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGt                                                 >  1:4346680/1‑150 (MQ=255)
                                                                                                   gTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCt                                               <  1:1755731/150‑1 (MQ=255)
                                                                                                       gCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATc                                           >  2:4093034/1‑150 (MQ=255)
                                                                                                          ggTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATcgcg                                        >  1:450948/1‑150 (MQ=255)
                                                                                                            ttCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGcc                                      <  1:3474555/150‑1 (MQ=255)
                                                                                                              ccACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTg                                    >  2:1887799/1‑150 (MQ=255)
                                                                                                              ccACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTg                                    <  2:3629592/150‑1 (MQ=255)
                                                                                                              ccACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTg                                    <  2:3401418/150‑1 (MQ=255)
                                                                                                               cACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTgg                                   <  1:3680945/150‑1 (MQ=255)
                                                                                                                aCCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgc                                  <  1:4569329/150‑1 (MQ=255)
                                                                                                                  cTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgcgc                                >  2:31525/1‑150 (MQ=255)
                                                                                                                   tGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGcgcg                               <  2:3972118/150‑1 (MQ=255)
                                                                                                                   tGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGcgcg                               <  2:3370709/150‑1 (MQ=255)
                                                                                                                       aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATg                           >  1:2630787/1‑150 (MQ=255)
                                                                                                                       aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATg                           >  1:2630957/1‑150 (MQ=255)
                                                                                                                       aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATg                           >  2:4410979/1‑150 (MQ=255)
                                                                                                                       aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATg                           >  1:2482402/1‑150 (MQ=255)
                                                                                                                        aCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATgg                          >  2:3113592/1‑150 (MQ=255)
                                                                                                                        aCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATgg                          >  2:3919021/1‑150 (MQ=255)
                                                                                                                         cTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGc                         >  1:4425354/1‑150 (MQ=255)
                                                                                                                         cTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGc                         <  1:3317430/150‑1 (MQ=255)
                                                                                                                         cTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGc                         >  2:1933772/1‑150 (MQ=255)
                                                                                                                          tGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACTCCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCg                        <  2:3208446/150‑1 (MQ=255)
                                                                                                                           ggATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGt                       >  2:798871/1‑150 (MQ=255)
                                                                                                                            gATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTa                      >  1:100710/1‑150 (MQ=255)
                                                                                                                            gATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTa                      >  1:1579673/1‑150 (MQ=255)
                                                                                                                             aTAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaa                     <  2:411143/150‑1 (MQ=255)
                                                                                                                             aTAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaa                     >  1:137441/1‑150 (MQ=255)
                                                                                                                               aGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaaaa                   >  2:4532968/1‑150 (MQ=255)
                                                                                                                                 ccAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAAt                 <  1:106346/150‑1 (MQ=255)
                                                                                                                                 ccAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAAt                 <  1:3663306/150‑1 (MQ=255)
                                                                                                                                 ccAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAAt                 <  2:2843670/150‑1 (MQ=255)
                                                                                                                                  caggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATg                >  1:1097563/1‑150 (MQ=255)
                                                                                                                                  caggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACAGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATg                >  1:4238850/1‑150 (MQ=255)
                                                                                                                                   aggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGc               <  1:3768367/150‑1 (MQ=255)
                                                                                                                                    ggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCa              <  1:80013/150‑1 (MQ=255)
                                                                                                                                       aggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtg           >  2:181914/1‑150 (MQ=255)
                                                                                                                                        gggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtgt          <  1:181799/150‑1 (MQ=255)
                                                                                                                                             gTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGAttt     <  1:947947/150‑1 (MQ=255)
                                                                                                                                               ataACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAt   >  1:2232137/1‑150 (MQ=255)
                                                                                                                                                taAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGACGAAGATCGTCTGAACGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAta  >  2:1996612/1‑149 (MQ=255)
                                                                                                                                                taAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTAtt  <  1:3923679/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATT  >  CP009273/2931570‑2931863

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: ‑