Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,424,363 T→A 100% T138S (ACC→TCC)  tsaC ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,424,3630TA98.0% 665.8 / 2.0 202T138S (ACC→TCC) tsaCtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base A (110/88);  total (111/91)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
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 203 total aligned reads displayed.

AGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTG  >  CP009273/3424218‑3424506
                                                                                                                                                 |                                                                                                                                               
  ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                             >  2:3171502/1‑148 (MQ=255)
   tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                            >  2:1078260/1‑148 (MQ=255)
   tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                            <  1:235522/148‑1 (MQ=255)
   tCACCAGTAAGGGCATGGGGTATTTCTGACGTATTGAAACGCCCCCCCGTTTCAACATGAACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCTTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGtaa                                                                                                                                            <  1:235559/148‑1 (MQ=255)
    cACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAAc                                                                                                                                           <  2:838700/148‑1 (MQ=255)
       ccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAg                                                                                                                                        >  2:778270/1‑148 (MQ=255)
           cAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGt                                                                                                                                    >  2:1360501/1‑148 (MQ=255)
           cAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGt                                                                                                                                    >  1:1527340/1‑148 (MQ=255)
            aGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGtt                                                                                                                                   <  2:3449288/148‑1 (MQ=255)
                  tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                             >  1:2636198/1‑148 (MQ=255)
                  tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                             >  2:1588331/1‑148 (MQ=255)
                  tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                             <  2:1643127/148‑1 (MQ=255)
                  tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                             <  2:2853413/148‑1 (MQ=255)
                   cgcgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATa                                                                                                                            <  2:860284/148‑1 (MQ=255)
                     cgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAAGTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAg                                                                                                                          <  2:255899/148‑1 (MQ=255)
                       gATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGcc                                                                                                                        >  1:3212693/1‑148 (MQ=255)
                       gATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGcc                                                                                                                        >  1:1404643/1‑148 (MQ=255)
                          ttCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTgg                                                                                                                     >  1:1882486/1‑148 (MQ=255)
                          ttCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTgg                                                                                                                     >  2:3038236/1‑148 (MQ=255)
                               acggATTTACACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                <  1:628923/146‑1 (MQ=255)
                               aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                <  2:1939788/148‑1 (MQ=255)
                               aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                >  1:1937641/1‑148 (MQ=255)
                               aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                <  1:1495135/148‑1 (MQ=255)
                                aGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                                                                               >  2:1980152/1‑148 (MQ=255)
                                 ggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACACAg                                                                                                              >  2:3208803/1‑148 (MQ=255)
                                 ggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                                                                              >  1:2315460/1‑148 (MQ=255)
                                 ggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                                                                              >  2:3104572/1‑148 (MQ=255)
                                 ggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                                                                              >  2:1335704/1‑148 (MQ=255)
                                 ggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                                                                              >  1:1075062/1‑148 (MQ=255)
                                   aTTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCa                                                                                                            <  1:481752/148‑1 (MQ=255)
                                     ttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAac                                                                                                          <  2:469358/148‑1 (MQ=255)
                                        aaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccac                                                                                                       <  1:2830648/148‑1 (MQ=255)
                                        aaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccac                                                                                                       <  2:2945500/148‑1 (MQ=255)
                                         cgccccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                                                                                      >  2:930631/4‑148 (MQ=255)
                                         aCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                                                                                      >  2:1295357/1‑148 (MQ=255)
                                          cGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAAccacca                                                                                                     >  1:2065297/1‑148 (MQ=255)
                                          cGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAAccacca                                                                                                     >  2:929875/1‑148 (MQ=255)
                                            cccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACCAAGCAACCACCAAc                                                                                                   >  2:12607/1‑148 (MQ=255)
                                            cccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAAc                                                                                                   >  2:356527/1‑148 (MQ=255)
                                               ccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGa                                                                                                <  1:1711588/148‑1 (MQ=255)
                                                cccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGAt                                                                                               >  1:3111065/1‑148 (MQ=255)
                                                cccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGAt                                                                                               >  1:1638703/1‑148 (MQ=255)
                                                cccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGAt                                                                                               >  2:1625092/1‑148 (MQ=255)
                                                  cGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATgg                                                                                             <  1:1294310/148‑1 (MQ=255)
                                                     ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCg                                                                                          <  2:3537016/148‑1 (MQ=255)
                                                      tCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCgg                                                                                         <  2:1568222/148‑1 (MQ=255)
                                                        aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                                                                       <  1:674180/148‑1 (MQ=255)
                                                              cacaACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                                                                                 >  1:50317/1‑148 (MQ=255)
                                                              cacaACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                                                                                 <  1:2499796/148‑1 (MQ=255)
                                                               acaACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTa                                                                                <  2:2411334/148‑1 (MQ=255)
                                                               acaACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTa                                                                                <  1:2074629/148‑1 (MQ=255)
                                                                caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTAc                                                                               <  1:2946252/148‑1 (MQ=255)
                                                                caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTAc                                                                               >  2:607018/1‑148 (MQ=255)
                                                                 aaCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACa                                                                              <  1:2379974/148‑1 (MQ=255)
                                                                   ccGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGc                                                                            >  1:1625835/1‑148 (MQ=255)
                                                                    cGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCa                                                                           <  2:1485059/148‑1 (MQ=255)
                                                                    cGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCa                                                                           <  2:3462102/148‑1 (MQ=255)
                                                                       gAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGc                                                                        >  2:806133/1‑148 (MQ=255)
                                                                        aaCGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCg                                                                       >  2:1003480/1‑148 (MQ=255)
                                                                         aCGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGa                                                                      >  1:829824/1‑148 (MQ=255)
                                                                                gcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAg                                                               <  2:584782/148‑1 (MQ=255)
                                                                                 ccAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGc                                                              <  2:1494540/148‑1 (MQ=255)
                                                                                     aTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc                                                          >  1:1814167/1‑148 (MQ=255)
                                                                                     aTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc                                                          >  1:1008724/1‑148 (MQ=255)
                                                                                      ttGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGccc                                                         >  2:576512/1‑148 (MQ=255)
                                                                                       tgtgCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCg                                                        >  2:49838/1‑148 (MQ=255)
                                                                                       tgtgCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCg                                                        >  2:2209388/1‑148 (MQ=255)
                                                                                                      tCTACTGTGCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtg                                         <  2:1853346/148‑1 (MQ=255)
                                                                                                       cTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtgt                                        >  2:2328363/1‑148 (MQ=255)
                                                                                                       cTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtgt                                        >  2:2640750/1‑148 (MQ=255)
                                                                                                         aCTGTTCGACATGGTGGCAATCCACTCAAGTGGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                      <  1:968819/148‑1 (MQ=255)
                                                                                                         aCTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                      >  1:2772493/1‑148 (MQ=255)
                                                                                                         aCTGTTCGACAAGGTGGCAATCCACTAAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCCCCAACGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                      <  1:1130193/148‑1 (MQ=255)
                                                                                                          cTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttgt                                     >  1:1659751/1‑148 (MQ=255)
                                                                                                                  cAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGgcg                             <  2:1814167/148‑1 (MQ=255)
                                                                                                                  cAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCATGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGgcg                             >  2:84057/1‑148 (MQ=255)
                                                                                                                    aGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGcgcg                           >  1:1603582/1‑148 (MQ=255)
                                                                                                                     ggTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCgg                          >  1:3371817/1‑148 (MQ=255)
                                                                                                                     ggTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGGTGTCGCAGGCGCgg                          >  2:3260725/1‑148 (MQ=255)
                                                                                                                     ggTGGCAATACACTCAAGTTGGCACTGTAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCgg                          >  2:2448527/1‑148 (MQ=255)
                                                                                                                      gggtcAATCCACTAATGTTGTCAGTGTAATAAACCAGCTGTTTACCATAACCCTGGCACAAACCAACCACCACCGGTTGGTCGGTTACTCGAACAGCACGCAAAACAACGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCggg                         <  1:3292361/144‑1 (MQ=255)
                                                                                                                      gTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCggg                         >  2:1703431/1‑148 (MQ=255)
                                                                                                                       tGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGa                        <  2:1638703/148‑1 (MQ=255)
                                                                                                                       tGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGa                        <  1:1475185/148‑1 (MQ=255)
                                                                                                                         gCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGaaa                      >  2:1639955/1‑148 (MQ=255)
                                                                                                                           aaTCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCTAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGa                    <  1:2314187/148‑1 (MQ=255)
                                                                                                                           aaTCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGa                    >  1:1774843/1‑148 (MQ=255)
                                                                                                                             tCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACa                  >  1:479667/1‑148 (MQ=255)
                                                                                                                                aCTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAg               >  2:3282067/1‑148 (MQ=255)
                                                                                                                                   cAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAACGACAAAGGTg            >  1:2224312/1‑148 (MQ=255)
                                                                                                                                    acgtgggCACTGTAAGAAACCTGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGTTGGTCGGTGAGTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGa           <  2:3269680/143‑1 (MQ=255)
                                                                                                                                    aaGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGa           <  1:3038236/148‑1 (MQ=255)
                                                                                                                                     aGTTGGCACTGGAAGAAACCAGCGGTTTACCATACGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGAc          <  2:2287132/148‑1 (MQ=255)
                                                                                                                                     aGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGACAGACAAAGGTGAc          >  2:2379683/1‑148 (MQ=255)
                                                                                                                                      gTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa         >  1:3393894/1‑148 (MQ=255)
                                                                                                                                      gTTGGCACTGGAAGAAAACAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCTAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa         <  1:2843910/148‑1 (MQ=255)
                                                                                                                                          gCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAc     <  1:2625104/148‑1 (MQ=255)
                                                                                                                                           cACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCACCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAcc    <  1:3555558/148‑1 (MQ=255)
                                                                                                                                            aCTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGACAGACAAAGGTGACAGGACCt   >  1:641837/1‑148 (MQ=255)
                                                                                                                                             cTGGAAGAAACCAGCGTTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTg  <  2:1527340/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTG  >  CP009273/3424218‑3424506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: