Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,926,831 G→T 100% D41Y (GAT→TAT)  rbsD → putative cytoplasmic sugar‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,926,8310GT99.4% 994.9 / ‑2.4 308D41Y (GAT→TAT) rbsDputative cytoplasmic sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (158/148);  total (158/150)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-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 308 total aligned reads displayed.

CGAAACGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAG  >  CP009273/3926684‑3926977
                                                                                                                                                   |                                                                                                                                                  
cGAAACGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATATGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTat                                                                                                                                                  <  1:5062803/150‑1 (MQ=255)
  aaaCGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTatat                                                                                                                                                <  2:3903113/150‑1 (MQ=255)
     cGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGc                                                                                                                                             <  2:1392211/150‑1 (MQ=255)
      gTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCa                                                                                                                                            >  1:1634434/1‑150 (MQ=255)
           gCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAAc                                                                                                                                       <  1:3532451/150‑1 (MQ=255)
             tGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAAccc                                                                                                                                     <  1:4567222/150‑1 (MQ=255)
              gATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCa                                                                                                                                    >  1:4580102/1‑150 (MQ=255)
              gATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCa                                                                                                                                    <  1:435921/150‑1 (MQ=255)
              gATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCa                                                                                                                                    <  1:3873797/150‑1 (MQ=255)
              gATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCa                                                                                                                                    <  2:2044229/150‑1 (MQ=255)
               aTGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAg                                                                                                                                   >  1:2137971/1‑150 (MQ=255)
                 ggAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAggg                                                                                                                                 <  2:1792804/150‑1 (MQ=255)
                   agaAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGgtg                                                                                                                               >  1:1905156/1‑150 (MQ=255)
                     aaaaaaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTa                                                                                                                             <  1:4789534/150‑1 (MQ=255)
                      aaaaaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTAc                                                                                                                            >  1:2417014/1‑150 (MQ=255)
                       aaaaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGGTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTAcc                                                                                                                           >  1:1218115/1‑150 (MQ=255)
                        aaaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACct                                                                                                                          <  1:368946/150‑1 (MQ=255)
                         aaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACctt                                                                                                                         >  1:1504097/1‑150 (MQ=255)
                         aaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACctt                                                                                                                         >  2:4638890/1‑150 (MQ=255)
                           aTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACcttct                                                                                                                       <  2:5318122/150‑1 (MQ=255)
                           aTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACcttct                                                                                                                       <  1:1020602/150‑1 (MQ=255)
                              aaaaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCtttt                                                                                                                    <  1:4020700/150‑1 (MQ=255)
                              aaaaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCtttt                                                                                                                    >  2:2536373/1‑150 (MQ=255)
                                aanaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTAt                                                                                                                  >  1:4294685/1‑150 (MQ=255)
                                aaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTAt                                                                                                                  >  1:1982553/1‑150 (MQ=255)
                                 aaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATg                                                                                                                 >  2:2786141/1‑150 (MQ=255)
                                 aaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATg                                                                                                                 >  1:2555131/1‑150 (MQ=255)
                                  aaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGc                                                                                                                >  1:740301/1‑150 (MQ=255)
                                    ggCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAg                                                                                                              >  2:4549640/1‑150 (MQ=255)
                                         cGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCt                                                                                                         >  2:2245326/1‑150 (MQ=255)
                                         cGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCt                                                                                                         >  2:2104215/1‑150 (MQ=255)
                                          gTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTg                                                                                                        >  2:4887328/1‑150 (MQ=255)
                                           ttCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTgg                                                                                                       <  1:4613717/150‑1 (MQ=255)
                                               tAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGcgt                                                                                                   >  2:4759056/1‑150 (MQ=255)
                                               tAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGcgt                                                                                                   <  2:4954086/150‑1 (MQ=255)
                                                aaTTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGcgtc                                                                                                  >  1:248806/1‑150 (MQ=255)
                                                      gATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCCTTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTcaca                                                                                            >  2:3147908/1‑150 (MQ=255)
                                                       atatTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaa                                                                                           >  1:1939915/1‑150 (MQ=255)
                                                             cATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaaatgaaa                                                                                     <  1:3181170/150‑1 (MQ=255)
                                                              aTCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaaatgaaat                                                                                    >  2:929918/1‑150 (MQ=255)
                                                              aTCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaaatgaaat                                                                                    <  1:3853726/150‑1 (MQ=255)
                                                              aTCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaaatgaaat                                                                                    <  1:1324390/150‑1 (MQ=255)
                                                              aTCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaaatgaaat                                                                                    >  1:234456/1‑150 (MQ=255)
                                                                cGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGc                                                                                  <  1:3607767/150‑1 (MQ=255)
                                                                cGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGG‑‑‑CGTCACAAATGAAATGCAgg                                                                               <  2:370072/150‑1 (MQ=255)
                                                                 ggTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCa                                                                                 >  2:164019/1‑150 (MQ=255)
                                                                  gTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAg                                                                                >  2:5224021/1‑150 (MQ=255)
                                                                   tGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAgg                                                                               >  1:5280470/1‑150 (MQ=255)
                                                                   tGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAgg                                                                               >  2:753984/1‑150 (MQ=255)
                                                                    gATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGt                                                                              >  2:261917/1‑150 (MQ=255)
                                                                     aTCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTc                                                                             <  1:4935548/150‑1 (MQ=255)
                                                                     aTCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTc                                                                             <  2:1214674/150‑1 (MQ=255)
                                                                      tctcCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACACATGAAATGCAGGTCg                                                                            >  1:3331094/1‑150 (MQ=255)
                                                                        tcCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAg                                                                          >  1:3063612/1‑150 (MQ=255)
                                                                         cccGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAgg                                                                         <  2:2505342/150‑1 (MQ=255)
                                                                         cccGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAgg                                                                         <  2:916410/150‑1 (MQ=255)
                                                                         accGTCTGGGACATACCGATACGCTTGTGGTGTGTGATGCTGGTTTACCCATACCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAgg                                                                         <  1:1037554/149‑1 (MQ=255)
                                                                             tCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAggcggc                                                                     >  2:943426/1‑150 (MQ=255)
                                                                               tGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCa                                                                   >  2:4351002/1‑150 (MQ=255)
                                                                               tGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCa                                                                   >  2:1071497/1‑150 (MQ=255)
                                                                               tGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCa                                                                   >  2:5164972/1‑150 (MQ=255)
                                                                                gggACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCAt                                                                  >  2:1669157/1‑150 (MQ=255)
                                                                                gggACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCAt                                                                  <  1:2664365/150‑1 (MQ=255)
                                                                                 ggACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCAtt                                                                 <  1:3217250/150‑1 (MQ=255)
                                                                                 ggACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCAtt                                                                 <  1:449394/150‑1 (MQ=255)
                                                                                   aCATTCCGTTTCGCGGGTGGGGTGTGTTTGTGGTTTACCCATCCCCAAAAGTACCACGCGTATCTTTTTGGGATTAACCCAGGGGGTTCCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTAt                                                               <  2:789400/150‑1 (MQ=255)
                                                                                    cATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGTTTTTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATc                                                              <  1:653271/150‑1 (MQ=255)
                                                                                    cATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATc                                                              <  1:4849194/150‑1 (MQ=255)
                                                                                      tACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATcgc                                                            <  2:2391339/150‑1 (MQ=255)
                                                                                       aCCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTACTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATAATctcg                                                           >  1:2708981/1‑150 (MQ=255)
                                                                                        ccGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCgg                                                          >  2:1737454/1‑150 (MQ=255)
                                                                                           aTACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAag                                                       >  1:3758705/1‑150 (MQ=255)
                                                                                             aCGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAagag                                                     <  2:1923736/150‑1 (MQ=255)
                                                                                             aCGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAagag                                                     <  2:3558360/150‑1 (MQ=255)
                                                                                               gCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGAt                                                   <  1:3659824/150‑1 (MQ=255)
                                                                                                cTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATc                                                  <  1:536729/150‑1 (MQ=255)
                                                                                                cTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATc                                                  <  1:4435233/150‑1 (MQ=255)
                                                                                                 tggtggTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCa                                                 <  1:2523134/150‑1 (MQ=255)
                                                                                                   gtggtgTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCaaa                                               >  1:290927/1‑150 (MQ=255)
                                                                                                    tggtgTGTGATTCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAac                                              >  1:2756422/1‑150 (MQ=255)
                                                                                                        gtgtgATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCa                                          <  2:3775904/150‑1 (MQ=255)
                                                                                                          gtgATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATa                                        >  2:3250968/1‑150 (MQ=255)
                                                                                                            gATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAAt                                      >  2:1422561/1‑150 (MQ=255)
                                                                                                             aTGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATc                                     <  1:1872173/150‑1 (MQ=255)
                                                                                                                         ccATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACg                         >  1:2772612/1‑150 (MQ=255)
                                                                                                                         ccATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACg                         <  1:259928/150‑1 (MQ=255)
                                                                                                                          cATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGa                        >  2:4511813/1‑150 (MQ=255)
                                                                                                                          cATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGa                        <  2:397957/150‑1 (MQ=255)
                                                                                                                           aTCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGaa                       <  1:4248077/150‑1 (MQ=255)
                                                                                                                                 aaaaGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCACAACTCCACGAAACGTTg                 >  2:4922952/1‑150 (MQ=255)
                                                                                                                                    aGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctc              >  2:4220744/1‑150 (MQ=255)
                                                                                                                                    aGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctc              >  2:4100949/1‑150 (MQ=255)
                                                                                                                                    aGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctc              >  1:2080734/1‑150 (MQ=255)
                                                                                                                                     gTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctca             >  2:2873224/1‑150 (MQ=255)
                                                                                                                                       aCAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcact           <  1:1812860/150‑1 (MQ=255)
                                                                                                                                         aaCGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcactca         <  2:3426174/150‑1 (MQ=255)
                                                                                                                                         aaCGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcactca         <  2:2959857/150‑1 (MQ=255)
                                                                                                                                             cgTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCtt     <  2:1566543/150‑1 (MQ=255)
                                                                                                                                             cgTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCtt     <  2:44745/150‑1 (MQ=255)
                                                                                                                                              gTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTg    <  1:3917465/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGAAACGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAG  >  CP009273/3926684‑3926977

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: