Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,926,831 G→T 92.1% 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,8310GT92.1% 748.6 / 54.9 254D41Y (GAT→TAT) rbsDputative cytoplasmic sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (8/12);  new base T (138/96);  total (146/108)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 254 total aligned reads displayed.

CGAAACGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAG  >  CP009273/3926684‑3926977
                                                                                                                                                   |                                                                                                                                                  
cGAAACGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTat                                                                                                                                                  >  2:1693261/1‑150 (MQ=255)
   aaCGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCTATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCGATATg                                                                                                                                               >  1:629785/1‑150 (MQ=255)
   aaCGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCTATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCGATATg                                                                                                                                               >  1:2023891/1‑150 (MQ=255)
   aaCGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATg                                                                                                                                               >  2:494819/1‑150 (MQ=255)
   aaCGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATg                                                                                                                                               <  2:3859155/150‑1 (MQ=255)
    aCGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATgg                                                                                                                                              <  1:3169593/150‑1 (MQ=255)
     cGTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGc                                                                                                                                             <  1:3442671/150‑1 (MQ=255)
      gTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCa                                                                                                                                            >  1:948788/1‑150 (MQ=255)
      gTTTCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCa                                                                                                                                            >  1:1278423/1‑150 (MQ=255)
         tCGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTa                                                                                                                                         <  2:3686047/150‑1 (MQ=255)
          cGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTaa                                                                                                                                        >  1:2218578/1‑150 (MQ=255)
          cGCTGATGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTaa                                                                                                                                        >  2:129629/1‑150 (MQ=255)
               aTGGAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAg                                                                                                                                   >  2:3921130/1‑150 (MQ=255)
                 ggAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAggg                                                                                                                                 >  2:1453342/1‑150 (MQ=255)
                 ggAGAAAAAAATGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAggg                                                                                                                                 >  2:3365118/1‑150 (MQ=255)
                      aaaaaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTAc                                                                                                                            <  1:2818021/150‑1 (MQ=255)
                        aaaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACct                                                                                                                          >  1:3046226/1‑150 (MQ=255)
                        aaaaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTAGATGGCATTAACCCAGGGTGTACct                                                                                                                          >  1:4495812/1‑150 (MQ=255)
                          aaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACcttc                                                                                                                        <  2:2985709/150‑1 (MQ=255)
                          aaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACcttc                                                                                                                        >  2:2557604/1‑150 (MQ=255)
                          aaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACcttc                                                                                                                        <  1:2140895/150‑1 (MQ=255)
                          aaTGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACcttc                                                                                                                        >  1:969340/1‑150 (MQ=255)
                            tGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACcttctt                                                                                                                      >  1:2857739/1‑150 (MQ=255)
                            tGAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACcttctt                                                                                                                      >  1:2241386/1‑150 (MQ=255)
                             gAAAAAAGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCttt                                                                                                                     <  1:3011010/150‑1 (MQ=255)
                              aaaaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCtttt                                                                                                                    >  1:1148453/1‑150 (MQ=255)
                              aaaaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCtttt                                                                                                                    <  1:2380715/150‑1 (MQ=255)
                              aaaaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCtttt                                                                                                                    <  1:2346323/150‑1 (MQ=255)
                              aaaaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCtttt                                                                                                                    >  1:1875642/1‑150 (MQ=255)
                               aaaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTa                                                                                                                   <  1:2767435/150‑1 (MQ=255)
                                aaaaGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTAt                                                                                                                  <  2:1704745/150‑1 (MQ=255)
                                   aGGCACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCa                                                                                                               >  1:3620765/1‑150 (MQ=255)
                                      cACCGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGt                                                                                                            >  2:2032375/1‑150 (MQ=255)
                                         cGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCt                                                                                                         >  2:1800592/1‑150 (MQ=255)
                                         cGTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCt                                                                                                         >  2:3510898/1‑150 (MQ=255)
                                          gTTCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTg                                                                                                        >  2:1317525/1‑150 (MQ=255)
                                           ttCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTgg                                                                                                       <  1:2433919/150‑1 (MQ=255)
                                            tCTTAATTCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTggg                                                                                                      >  2:2370719/1‑150 (MQ=255)
                                                   tCTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCgtcgtc                                                                                               <  2:618506/150‑1 (MQ=255)
                                                    cTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTca                                                                                              >  1:526121/1‑150 (MQ=255)
                                                    cTGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTca                                                                                              >  1:1020834/1‑150 (MQ=255)
                                                     tGATATTTCATCGGTGATCTCCCGTCTGGGNCATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTcac                                                                                             >  1:4276574/1‑150 (MQ=255)
                                                     tGATATTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTcac                                                                                             >  1:140416/1‑150 (MQ=255)
                                                       atatTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaa                                                                                           >  1:3105544/1‑150 (MQ=255)
                                                       atatTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaa                                                                                           >  2:244316/1‑150 (MQ=255)
                                                         atTTCATCGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACaaat                                                                                         <  2:2243350/150‑1 (MQ=255)
                                                           ttCATCGGTGATCTCACGTCTGGGACAAACCGATAAGCTGGAGGAGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGAGTATCTAAATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCAGGGCGACGTcacacaaga                                                                                       >  2:2265707/1‑145 (MQ=255)
                                                                cGGTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGc                                                                                  <  1:3234907/150‑1 (MQ=255)
                                                                 ggTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCa                                                                                 >  1:4590963/1‑150 (MQ=255)
                                                                 ggTGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCa                                                                                 <  2:2376788/150‑1 (MQ=255)
                                                                   tGATCTCCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAgg                                                                               >  1:1839049/1‑150 (MQ=255)
                                                                      tctcCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCg                                                                            <  2:2781717/150‑1 (MQ=255)
                                                                       ctcCCGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGa                                                                           >  2:3109118/1‑150 (MQ=255)
                                                                          ccGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAggc                                                                        <  2:103511/150‑1 (MQ=255)
                                                                           cGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAggcg                                                                       <  2:4311675/150‑1 (MQ=255)
                                                                           cGTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAggcg                                                                       <  2:1985545/150‑1 (MQ=255)
                                                                            gTCTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAggcgg                                                                      >  2:3816065/1‑150 (MQ=255)
                                                                              cTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGcc                                                                    >  1:721332/1‑150 (MQ=255)
                                                                              cTGGGACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGcc                                                                    >  1:4192744/1‑150 (MQ=255)
                                                                                  gACATACCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTa                                                                <  1:855995/150‑1 (MQ=255)
                                                                                       aCCGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATcgcg                                                           >  1:573214/1‑150 (MQ=255)
                                                                                         cGATACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGa                                                         >  1:3404002/1‑150 (MQ=255)
                                                                                           aTACGCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAag                                                       >  2:693256/1‑150 (MQ=255)
                                                                                               gCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGAt                                                   >  1:853885/1‑150 (MQ=255)
                                                                                               gCTGGTGGTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGAt                                                   >  1:1703195/1‑150 (MQ=255)
                                                                                                 tggtggTGTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCa                                                 >  2:407765/1‑150 (MQ=255)
                                                                                                   gtggtgTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCaaa                                               <  1:4523976/150‑1 (MQ=255)
                                                                                                    tggtgTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAac                                              >  2:1055054/1‑150 (MQ=255)
                                                                                                    tggtgTGTGATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAac                                              <  1:2539546/150‑1 (MQ=255)
                                                                                                        gtgtgATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCa                                          <  2:2745796/150‑1 (MQ=255)
                                                                                                            gATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAAt                                      <  2:4070946/150‑1 (MQ=255)
                                                                                                            gATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAAt                                      <  2:2213083/150‑1 (MQ=255)
                                                                                                            gATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAAt                                      <  1:744180/150‑1 (MQ=255)
                                                                                                            gATGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAAt                                      <  1:4364152/150‑1 (MQ=255)
                                                                                                             aTGCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATc                                     <  1:3903901/150‑1 (MQ=255)
                                                                                                               gCTGGTTTACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCg                                   >  1:4235112/1‑150 (MQ=255)
                                                                                                                   gTTTACCCATCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAAc                               <  2:1953407/150‑1 (MQ=255)
                                                                                                                   gTTTACCCATCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAAc                               <  2:786180/150‑1 (MQ=255)
                                                                                                                      tACCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTcc                            >  1:3553220/1‑150 (MQ=255)
                                                                                                                       aCCCATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCa                           <  1:3860714/150‑1 (MQ=255)
                                                                                                                         ccATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACg                         >  1:4299136/1‑150 (MQ=255)
                                                                                                                         ccATCCCCAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACg                         >  1:983258/1‑150 (MQ=255)
                                                                                                                           aTCCCCAAAAGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGaa                       <  1:3645352/150‑1 (MQ=255)
                                                                                                                             ccccAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAAc                     <  1:1743741/150‑1 (MQ=255)
                                                                                                                             ccccAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAAc                     <  2:2861338/150‑1 (MQ=255)
                                                                                                                             ccccAAAAGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAAc                     >  2:2224610/1‑150 (MQ=255)
                                                                                                                                 aaaaGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTg                 >  1:1656744/1‑150 (MQ=255)
                                                                                                                                 aaaaGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTg                 >  2:3958611/1‑150 (MQ=255)
                                                                                                                                    aGTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctc              >  2:3655776/1‑150 (MQ=255)
                                                                                                                                    aGTACAACGCGTATCGATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAACTGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCCCGAAACGTTGctc              <  2:1901174/150‑1 (MQ=255)
                                                                                                                                     gTACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctca             >  1:1761675/1‑150 (MQ=255)
                                                                                                                                      tACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcac            <  2:1511688/150‑1 (MQ=255)
                                                                                                                                      tACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcac            >  2:800301/1‑150 (MQ=255)
                                                                                                                                      tACAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcac            <  2:891924/150‑1 (MQ=255)
                                                                                                                                       aCAACGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcact           >  2:580075/1‑150 (MQ=255)
                                                                                                                                         aaCGCGTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGctcactca         >  2:4126835/1‑150 (MQ=255)
                                                                                                                                           cgcgTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCAcc       <  2:4618390/150‑1 (MQ=255)
                                                                                                                                             cgTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGTAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCtt     >  2:4051395/1‑150 (MQ=255)
                                                                                                                                             cgTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCtt     >  2:3184960/1‑150 (MQ=255)
                                                                                                                                             cgTATCTATATGGCATTAACCCAGGGTGTACCTTCTTTTATGCAGGTGCTGGGCGTCGTCACAAATGAAATGCAGGTCGAGGCGGCCATTATCGCGGAAGAGATCAAACACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCtt     <  1:1969586/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: