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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,931,7170CA99.0% 640.6 / ‑1.9 198N307K (AAC→AAAfucKL‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (125/71);  total (125/73)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 200 total aligned reads displayed.

GGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTA  >  CP009273/2931571‑2931861
                                                                                                                                                  |                                                                                                                                                
ggCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAcc                                                                                                                                                 <  2:2197574/148‑1 (MQ=255)
 gcgcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAacc                                                                                                                                                >  2:2166710/1‑148 (MQ=255)
  cgcTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCa                                                                                                                                               <  2:2866960/148‑1 (MQ=255)
   gcTTGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAg                                                                                                                                              <  1:328481/148‑1 (MQ=255)
            gAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAt                                                                                                                                     >  1:3444356/1‑148 (MQ=255)
              aCAAAATGAACCCGTGCTCTCTTCCGGTACTTGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGTATAGCAAGGCAGGGTTGTATTAACCAGGTATTCAAtgg                                                                                                                                   <  2:2353346/148‑1 (MQ=255)
              aCAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAAtgg                                                                                                                                   <  1:2467211/148‑1 (MQ=255)
                 aaaTGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCACtggctg                                                                                                                                >  2:3031203/1‑148 (MQ=255)
                    tGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCa                                                                                                                             >  2:2067883/1‑148 (MQ=255)
                      aaCCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATc                                                                                                                           <  2:2230865/148‑1 (MQ=255)
                         ccGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCgg                                                                                                                        >  1:2997873/1‑148 (MQ=255)
                             gCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTg                                                                                                                    <  2:2298446/148‑1 (MQ=255)
                             gCTCTCTTCCGGTACATGGGAAATTTTAATGGGTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTg                                                                                                                    <  2:571982/148‑1 (MQ=255)
                              ctctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGc                                                                                                                   <  1:1894170/148‑1 (MQ=255)
                              ctctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGc                                                                                                                   >  2:3124901/1‑148 (MQ=255)
                              ctctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGc                                                                                                                   <  1:2669673/148‑1 (MQ=255)
                                ctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTg                                                                                                                 >  2:2370023/1‑148 (MQ=255)
                                ctctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTg                                                                                                                 >  2:618285/1‑148 (MQ=255)
                                 tctTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTgg                                                                                                                >  1:1755141/1‑148 (MQ=255)
                                  cctCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                               <  1:132078/146‑1 (MQ=255)
                                    tCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGact                                                                                                             >  2:3247408/1‑146 (MQ=255)
                                    tCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAAt                                                                                                             >  2:3030207/1‑148 (MQ=255)
                                     ccGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCATTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATACGGCGTGCTGGAATg                                                                                                            <  2:1560385/148‑1 (MQ=255)
                                     ccGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATg                                                                                                            <  1:1609495/148‑1 (MQ=255)
                                     ccGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATg                                                                                                            <  2:2219280/148‑1 (MQ=255)
                                     ccGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATg                                                                                                            <  2:3421272/148‑1 (MQ=255)
                                       ggTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATggg                                                                                                          <  1:3251711/148‑1 (MQ=255)
                                           cATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTgag                                                                                                      >  2:737732/1‑148 (MQ=255)
                                           cATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTgag                                                                                                      >  1:960373/1‑148 (MQ=255)
                                              gggAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGaaa                                                                                                   >  1:78503/1‑148 (MQ=255)
                                              gggAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGaaa                                                                                                   >  2:2963044/1‑148 (MQ=255)
                                               gtaaaTTTTAATGGTTCGCAGCGCCCAGGTTTATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTTGAATGGGTGAGaaaa                                                                                                  <  2:1837221/146‑1 (MQ=255)
                                                gAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                               <  2:1259321/134‑1 (MQ=255)
                                                gAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGa                                                                                                               >  1:1259321/1‑134 (MQ=255)
                                                gAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAAc                                                                                                 >  1:3508469/1‑148 (MQ=255)
                                                  aaTTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTg                                                                                               <  2:3127496/148‑1 (MQ=255)
                                                  aaTTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTg                                                                                               <  2:2826732/148‑1 (MQ=255)
                                                  aaTTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTg                                                                                               <  1:1539259/148‑1 (MQ=255)
                                                    ttttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGtt                                                                                             <  1:1484515/148‑1 (MQ=255)
                                                    ttttAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGtt                                                                                             >  2:888240/1‑148 (MQ=255)
                                                        aaTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGTAATGGGTGAGAAAACTGTTCTgg                                                                                         >  2:96523/1‑148 (MQ=255)
                                                        aaTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTgg                                                                                         <  2:457101/148‑1 (MQ=255)
                                                         aTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGa                                                                                        >  1:2375718/1‑148 (MQ=255)
                                                         aTGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGa                                                                                        >  2:1214171/1‑148 (MQ=255)
                                                              tCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCt                                                                                   >  2:705244/1‑148 (MQ=255)
                                                                   gcgcCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAac                                                                              <  1:2228010/148‑1 (MQ=255)
                                                                    cgcCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAaca                                                                             >  1:862153/1‑148 (MQ=255)
                                                                    cgcCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAaca                                                                             >  2:1774059/1‑148 (MQ=255)
                                                                     gcCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAacac                                                                            <  1:1146067/148‑1 (MQ=255)
                                                                         aGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                        >  2:2557300/1‑148 (MQ=255)
                                                                         aGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                        >  2:1502284/1‑148 (MQ=255)
                                                                         aGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTg                                                                        >  2:1501857/1‑148 (MQ=255)
                                                                          ggTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTgg                                                                       >  1:1742533/1‑148 (MQ=255)
                                                                          ggTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTgg                                                                       >  2:652033/1‑148 (MQ=255)
                                                                           gTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGc                                                                      >  1:3412413/1‑148 (MQ=255)
                                                                           gTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGc                                                                      >  1:3488936/1‑148 (MQ=255)
                                                                             tGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCaa                                                                    >  1:2407652/1‑148 (MQ=255)
                                                                             tGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCaa                                                                    <  2:3488936/148‑1 (MQ=255)
                                                                              gATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCaaa                                                                   <  2:1865040/148‑1 (MQ=255)
                                                                               aTACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAAt                                                                  >  2:895224/1‑148 (MQ=255)
                                                                                tACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATg                                                                 >  2:1319063/1‑148 (MQ=255)
                                                                                 aCTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGt                                                                >  2:1348137/1‑148 (MQ=255)
                                                                                  cTTCTCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGtt                                                               <  2:849795/148‑1 (MQ=255)
                                                                                  cTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGtt                                                               >  1:396042/1‑148 (MQ=255)
                                                                                    tCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGttga                                                             >  2:785871/1‑148 (MQ=255)
                                                                                         gTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaa                                                        >  1:2695605/1‑148 (MQ=255)
                                                                                         gTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaa                                                        <  1:245509/148‑1 (MQ=255)
                                                                                            aaGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaagaa                                                     >  1:1307136/1‑148 (MQ=255)
                                                                                            aaGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTgaagaa                                                     >  1:576831/1‑148 (MQ=255)
                                                                                                  gTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGt                                               >  1:2064098/1‑148 (MQ=255)
                                                                                                  gTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGt                                               <  1:1816628/148‑1 (MQ=255)
                                                                                                     cGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTg                                            >  2:3159076/1‑148 (MQ=255)
                                                                                                       ccGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGAt                                          <  1:2979035/148‑1 (MQ=255)
                                                                                                       ccGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGAt                                          <  1:2979032/148‑1 (MQ=255)
                                                                                                         ggTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATcg                                        >  1:806782/1‑148 (MQ=255)
                                                                                                           ttCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATcgcg                                      >  2:2797959/1‑148 (MQ=255)
                                                                                                            tCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCgcgc                                     >  1:2904071/1‑148 (MQ=255)
                                                                                                              cACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCt                                   >  2:1607847/1‑148 (MQ=255)
                                                                                                               aCCTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTg                                  >  2:1329285/1‑148 (MQ=255)
                                                                                                                 cTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgc                                >  1:3137212/1‑148 (MQ=255)
                                                                                                                 cTGCGAACTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGgc                                >  1:1362144/1‑148 (MQ=255)
                                                                                                                      aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGa                           >  1:3539063/1‑148 (MQ=255)
                                                                                                                      aaCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGa                           <  2:1118619/148‑1 (MQ=255)
                                                                                                                       aCTGGATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGAt                          >  2:2074449/1‑148 (MQ=255)
                                                                                                                           gATAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCg                      <  1:2067699/148‑1 (MQ=255)
                                                                                                                            aTAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGt                     <  2:3372861/148‑1 (MQ=255)
                                                                                                                            aTAGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGt                     >  2:406864/1‑148 (MQ=255)
                                                                                                                              aGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaa                   >  1:2130058/1‑148 (MQ=255)
                                                                                                                              aGCCAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaa                   >  1:2133870/1‑148 (MQ=255)
                                                                                                                                ccAGGCAGGGTTGTATAAACCAGGTATGCACTGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaaaa                 >  2:818740/1‑148 (MQ=255)
                                                                                                                                ccAGGCAGGGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTaaaa                 >  1:521666/1‑148 (MQ=255)
                                                                                                                                  aggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGNCGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAAt               <  1:1809647/148‑1 (MQ=255)
                                                                                                                                  aggcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAAt               >  2:1188067/1‑148 (MQ=255)
                                                                                                                                    gcaggGTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGACTGATCGCGCCTGGCGCGGATGGCGTAAAAATGc             >  2:506704/1‑148 (MQ=255)
                                                                                                                                       gggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgt          <  1:2208569/148‑1 (MQ=255)
                                                                                                                                       gggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgt          >  1:3465436/1‑148 (MQ=255)
                                                                                                                                        ggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtg         >  2:1931934/1‑148 (MQ=255)
                                                                                                                                        ggTTGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAgtg         >  1:1750349/1‑148 (MQ=255)
                                                                                                                                           tGTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGa      <  2:20648/148‑1 (MQ=255)
                                                                                                                                            gTATAAACCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGAt     <  1:1596150/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGCGCTGGTGCTGAACAAAATGAACCCGTGCTCTCTTCCGGTACATGGGAAATTTTAATGGTTCGCAGCGCCCAGGTTGATACTTCGCTGTTAAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTA  >  CP009273/2931571‑2931861

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: