Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,471,930 T→A 100% P575P (CCA→CCT valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,471,9300TA96.2% 1004.6 / 17.4 308P575P (CCA→CCTvalSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (2/10);  new base A (152/144);  total (154/154)
Fisher's exact test for biased strand distribution p-value = 3.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 308 total aligned reads displayed.

GGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTTG  >  CP009273/4471783‑4472074
                                                                                                                                                   |                                                                                                                                                
ggTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAgg                                                                                                                                                >  2:3617976/1‑150 (MQ=255)
 gTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGc                                                                                                                                               <  1:4716042/150‑1 (MQ=255)
 gTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGc                                                                                                                                               <  1:549563/150‑1 (MQ=255)
 gTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGc                                                                                                                                               <  2:332533/150‑1 (MQ=255)
    gACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGc                                                                                                                                            >  1:4746708/1‑150 (MQ=255)
    gACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGc                                                                                                                                            >  1:4487066/1‑150 (MQ=255)
     aCGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCg                                                                                                                                           >  2:4586744/1‑150 (MQ=255)
     aCGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCg                                                                                                                                           <  1:3087539/150‑1 (MQ=255)
      cGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGa                                                                                                                                          >  1:488034/1‑150 (MQ=255)
       gCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGaa                                                                                                                                         >  1:3459496/1‑150 (MQ=255)
         cAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAAt                                                                                                                                       >  2:1025549/1‑150 (MQ=255)
          aGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATa                                                                                                                                      <  2:2968590/150‑1 (MQ=255)
          aGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATa                                                                                                                                      <  2:3278976/150‑1 (MQ=255)
           gcgcCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATAc                                                                                                                                     >  1:2728160/1‑150 (MQ=255)
                gccAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCa                                                                                                                                <  1:2540995/150‑1 (MQ=255)
                    gggTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCa                                                                                                                            >  2:3177363/1‑150 (MQ=255)
                     ggTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCat                                                                                                                           >  2:3033680/1‑150 (MQ=255)
                     ggTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCat                                                                                                                           <  2:4488926/150‑1 (MQ=255)
                     ggTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAAGACCGTTCGGGAACTGCGTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCat                                                                                                                           <  2:1821696/150‑1 (MQ=255)
                       tGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCatat                                                                                                                         >  1:1640077/1‑150 (MQ=255)
                            cgcagcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGt                                                                                                                    >  2:5053277/1‑150 (MQ=255)
                            cgcagcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGt                                                                                                                    >  2:2488470/1‑150 (MQ=255)
                              cagcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTgg                                                                                                                  >  1:1333714/1‑150 (MQ=255)
                              cagcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTgg                                                                                                                  >  1:4001460/1‑150 (MQ=255)
                                gcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGt                                                                                                                <  2:4995794/150‑1 (MQ=255)
                                gcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGt                                                                                                                <  2:3956529/150‑1 (MQ=255)
                                 cgcgTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTc                                                                                                               >  1:5060192/1‑150 (MQ=255)
                                   cgTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGa                                                                                                             <  2:825312/150‑1 (MQ=255)
                                    gTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGAt                                                                                                            >  2:1593671/1‑150 (MQ=255)
                                    gTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGAt                                                                                                            <  1:1909167/150‑1 (MQ=255)
                                     tCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATa                                                                                                           <  2:3412320/150‑1 (MQ=255)
                                       aGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAAc                                                                                                         >  1:403056/1‑150 (MQ=255)
                                        gTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACg                                                                                                        >  2:3128928/1‑150 (MQ=255)
                                           ccGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTa                                                                                                     >  2:3785980/1‑150 (MQ=255)
                                                  gCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTgg                                                                                              <  2:1171010/150‑1 (MQ=255)
                                                  gCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTgg                                                                                              >  1:1505004/1‑150 (MQ=255)
                                                  gCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTgg                                                                                              <  1:4808195/150‑1 (MQ=255)
                                                   cTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTggn                                                                                             <  1:5129332/150‑2 (MQ=255)
                                                   cTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGa                                                                                             >  1:3217269/1‑150 (MQ=255)
                                                    tCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAt                                                                                            >  1:2824788/1‑150 (MQ=255)
                                                        tACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTgg                                                                                        <  1:789751/150‑1 (MQ=255)
                                                           cGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACa                                                                                     <  1:4276455/150‑1 (MQ=255)
                                                           cGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACa                                                                                     <  1:3496861/150‑1 (MQ=255)
                                                             ttagggAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtc                                                                                   <  2:636465/147‑1 (MQ=255)
                                                             ttCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtc                                                                                   >  2:1652315/1‑150 (MQ=255)
                                                             ttCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtc                                                                                   <  1:2667383/150‑1 (MQ=255)
                                                               cGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtctt                                                                                 <  1:1411935/150‑1 (MQ=255)
                                                                gggAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtcttc                                                                                <  1:3818564/150‑1 (MQ=255)
                                                                   aaCTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTgg                                                                             >  2:68393/1‑150 (MQ=255)
                                                                    aCTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGc                                                                            <  1:2469813/150‑1 (MQ=255)
                                                                        cTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTc                                                                        <  2:4146153/150‑1 (MQ=255)
                                                                         ttCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCg                                                                       >  2:4893896/1‑150 (MQ=255)
                                                                         ttCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCg                                                                       >  2:4419889/1‑150 (MQ=255)
                                                                          tctcGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGt                                                                      <  2:4757535/150‑1 (MQ=255)
                                                                           ctcGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGtc                                                                     <  1:2084549/150‑1 (MQ=255)
                                                                            tcGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGtca                                                                    <  1:2606104/150‑1 (MQ=255)
                                                                                 gcgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCGAGCAGTTCAGGCAGCGAAATTCCGTCAACCATATCCAGTGGGTCGATATCGTTACCCTTGGATTTTGACATCTTCTGGCCTTCGTCTTCTCg                                                               <  2:2293263/150‑1 (MQ=255)
                                                                                 gcgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCAca                                                               <  1:3722647/150‑2 (MQ=255)
                                                                                 gcgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACg                                                               >  1:4534622/1‑150 (MQ=255)
                                                                                 gcgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACg                                                               <  2:1597658/150‑1 (MQ=255)
                                                                                 gcgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACg                                                               >  2:131055/1‑150 (MQ=255)
                                                                                 gcgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACg                                                               <  1:4773305/150‑1 (MQ=255)
                                                                                  cgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGa                                                              >  2:2202822/1‑150 (MQ=255)
                                                                                   gcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGaa                                                             <  1:354934/150‑1 (MQ=255)
                                                                                     ttACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATc                                                           >  2:2810845/1‑150 (MQ=255)
                                                                                         ggATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGc                                                       >  1:2143657/1‑150 (MQ=255)
                                                                                          gATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGcc                                                      >  1:3364233/1‑150 (MQ=255)
                                                                                               tGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCa                                                 <  2:1700102/150‑1 (MQ=255)
                                                                                                   cGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTa                                             <  1:1593671/150‑1 (MQ=255)
                                                                                                    gCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTaa                                            <  2:4746708/150‑1 (MQ=255)
                                                                                                    gCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTaa                                            >  1:3913796/1‑150 (MQ=255)
                                                                                                      cAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAAc                                          <  1:2747634/150‑1 (MQ=255)
                                                                                                        gCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACgg                                        <  2:2302555/150‑1 (MQ=255)
                                                                                                         cTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCATCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGACGGTCATGTAAACGgt                                       <  2:2439434/150‑1 (MQ=255)
                                                                                                         cTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGgt                                       <  1:922302/150‑1 (MQ=255)
                                                                                                          tGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGgtg                                      >  1:1941515/1‑150 (MQ=255)
                                                                                                            cGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGtgtg                                    >  1:1804864/1‑150 (MQ=255)
                                                                                                             ggCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTAGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTgg                                   >  1:1211612/1‑150 (MQ=255)
                                                                                                              gCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGa                                  >  1:177095/1‑150 (MQ=255)
                                                                                                                tGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAAc                                >  1:469100/1‑150 (MQ=255)
                                                                                                                tGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAAc                                <  2:1395995/150‑1 (MQ=255)
                                                                                                                 tcatcatATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACg                               >  2:4446104/2‑150 (MQ=255)
                                                                                                                 gCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACg                               <  1:2373140/150‑1 (MQ=255)
                                                                                                                     catATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCAc                           >  1:1300988/1‑150 (MQ=255)
                                                                                                                      atatTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCAcc                          >  2:1974069/1‑150 (MQ=255)
                                                                                                                        atTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTg                        <  2:3846844/150‑1 (MQ=255)
                                                                                                                           gCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCgg                     <  2:5067119/150‑1 (MQ=255)
                                                                                                                           gCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCgg                     >  1:1165085/1‑150 (MQ=255)
                                                                                                                             cGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTTGACATCTTCTGGCCTTAGTCAACACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGAGGtt                   >  1:1748075/1‑150 (MQ=255)
                                                                                                                             cGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGtt                   <  2:1752847/150‑1 (MQ=255)
                                                                                                                               gTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTg                 <  1:1067547/150‑1 (MQ=255)
                                                                                                                                tACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGc                >  2:4242817/1‑150 (MQ=255)
                                                                                                                                tACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGc                >  2:3821507/1‑150 (MQ=255)
                                                                                                                                     ttttCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCAttt           <  1:4343873/150‑1 (MQ=255)
                                                                                                                                      tttCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCAtttt          <  2:527031/150‑1 (MQ=255)
                                                                                                                                       ttCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTc         <  2:465110/150‑1 (MQ=255)
                                                                                                                                           agcagTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCt     <  2:2712973/150‑1 (MQ=255)
                                                                                                                                             cagTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCttt   <  2:1505004/150‑1 (MQ=255)
                                                                                                                                             cagTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCttt   >  1:4994520/1‑150 (MQ=255)
                                                                                                                                              agTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTtg  <  1:2972786/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTTG  >  CP009273/4471783‑4472074

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: