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,3630TA97.9% 624.8 / 2.2 192T138S (ACC→TCC) tsaCtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base A (103/85);  total (106/86)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-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 194 total aligned reads displayed.

AGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGG  >  CP009273/3424218‑3424507
                                                                                                                                                 |                                                                                                                                                
  ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                              >  1:343620/1‑148 (MQ=255)
  ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCAACTCAAGTTGGCACTGgaagaa                                                                                                                                              >  1:344788/1‑148 (MQ=255)
   tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGaaa                                                                                                                                             >  2:190337/1‑148 (MQ=255)
   tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                             <  2:1287088/148‑1 (MQ=255)
   tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                             >  1:3130730/1‑148 (MQ=255)
   tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                             <  2:3498523/148‑1 (MQ=255)
   tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                             >  1:1743385/1‑148 (MQ=255)
    cACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAAc                                                                                                                                            >  1:352589/1‑148 (MQ=255)
     aCCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAAcc                                                                                                                                           >  1:1345075/1‑148 (MQ=255)
           cAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAg                                                                                                                                         <  1:820517/144‑1 (MQ=255)
           cAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAg                                                                                                                                         >  2:820259/1‑144 (MQ=255)
            aGGGCATCGCGGATTTCTGAAGTATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGtt                                                                                                                                    <  1:2486386/148‑1 (MQ=255)
              ggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTa                                                                                                                                  >  1:671028/1‑148 (MQ=255)
              ggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTa                                                                                                                                  >  1:2229967/1‑148 (MQ=255)
                cATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTAcc                                                                                                                                <  1:1502633/148‑1 (MQ=255)
                  tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                              >  2:191555/1‑148 (MQ=255)
                  tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                              >  2:191697/1‑148 (MQ=255)
                  tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCAt                                                                                                                              >  2:1991528/1‑148 (MQ=255)
                    gcgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATaa                                                                                                                            <  1:1490026/148‑1 (MQ=255)
                       gATTTCTTATGGATTTAAACTCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGcc                                                                                                                         <  2:2993371/148‑1 (MQ=255)
                       gATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGcc                                                                                                                         <  2:1966983/148‑1 (MQ=255)
                        aTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCt                                                                                                                        <  1:462146/148‑1 (MQ=255)
                        aTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCt                                                                                                                        >  2:1187750/1‑148 (MQ=255)
                               aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                 <  1:219887/148‑1 (MQ=255)
                                 ggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                                                                               <  2:155403/148‑1 (MQ=255)
                                 ggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                                                                               >  2:140307/1‑148 (MQ=255)
                                   aTTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCa                                                                                                             >  2:2398068/1‑148 (MQ=255)
                                        aaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccac                                                                                                        <  1:3339588/148‑1 (MQ=255)
                                        aaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccac                                                                                                        <  2:2471681/148‑1 (MQ=255)
                                         aCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                                                                                       >  2:772719/1‑148 (MQ=255)
                                         aCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                                                                                       >  1:1636382/1‑148 (MQ=255)
                                          cGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAAccacca                                                                                                      >  1:2422973/1‑148 (MQ=255)
                                          cGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAAccacca                                                                                                      >  1:832400/1‑148 (MQ=255)
                                          cGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAAccacca                                                                                                      >  2:640387/1‑148 (MQ=255)
                                          cGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAAccacca                                                                                                      >  1:3328884/1‑148 (MQ=255)
                                               ccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGa                                                                                                 <  1:2460095/148‑1 (MQ=255)
                                               ccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGa                                                                                                 <  1:2742529/148‑1 (MQ=255)
                                                cccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGAt                                                                                                >  2:373257/1‑148 (MQ=255)
                                                cccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGAt                                                                                                >  2:3278605/1‑148 (MQ=255)
                                                cccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGAt                                                                                                <  1:2856652/148‑1 (MQ=255)
                                                 ccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATg                                                                                               >  2:345826/1‑148 (MQ=255)
                                                  cGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATgg                                                                                              >  2:661483/1‑148 (MQ=255)
                                                     ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCg                                                                                           >  1:1303864/1‑148 (MQ=255)
                                                        aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                                                                        >  2:3041866/1‑148 (MQ=255)
                                                        aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                                                                        >  1:760829/1‑148 (MQ=255)
                                                             gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                   >  2:893/1‑148 (MQ=255)
                                                             gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                   >  2:1873542/1‑148 (MQ=255)
                                                             gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                   >  1:3446589/1‑148 (MQ=255)
                                                             gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                   >  2:2606071/1‑148 (MQ=255)
                                                                caACCGGGAACGCCGCGCCAACTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTAc                                                                                <  2:106433/148‑1 (MQ=255)
                                                                caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTAc                                                                                >  2:1715682/1‑148 (MQ=255)
                                                                caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTAc                                                                                <  2:106411/148‑1 (MQ=255)
                                                                caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTAc                                                                                <  2:107006/148‑1 (MQ=255)
                                                                   ccGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGc                                                                             >  2:1993533/1‑148 (MQ=255)
                                                                    cGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCa                                                                            >  2:2529199/1‑148 (MQ=255)
                                                                               cgcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCaaa                                                                 <  2:1743900/148‑1 (MQ=255)
                                                                                   aaaTTGTGCGCGATTTTCGTCTACTTTTCTAAATGTTGTAAATCAACTAAAGTTGTCACTGGAATAAACCAGCGGTTTACCATAAGCCGGGCACAATGCAACAACCAACGGTTGTTCGGTGACTCGTACAGCAAGCTAATCAAAGCgg                                                             <  2:3426508/148‑1 (MQ=255)
                                                                                     aTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc                                                           >  2:287870/1‑148 (MQ=255)
                                                                                     aTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc                                                           >  1:2504597/1‑148 (MQ=255)
                                                                                       tgtgCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCg                                                         >  2:1301594/1‑148 (MQ=255)
                                                                                       tgtgCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCg                                                         >  2:1675598/1‑148 (MQ=255)
                                                                                        gtgCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGt                                                        <  2:1129080/148‑1 (MQ=255)
                                                                                          gcgcGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCa                                                      >  1:1748739/1‑148 (MQ=255)
                                                                                               aaCTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAg                                                 >  1:590107/1‑148 (MQ=255)
                                                                                               aaCTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAg                                                 >  1:3221416/1‑148 (MQ=255)
                                                                                                 attCGTCTACTGTTCTACATGGTGTCAATCCACTCAAGTTGGCACTGTAATAAACCAGGGTTTTACCATAAGCCTGGCACAACGCAACCACCAACGTTTGGTCGTTGACTCGTACCGCAAGCTAATCAAAGCGGCCCGTCAACCAgcg                                               <  1:2728999/147‑1 (MQ=255)
                                                                                                    cGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCgg                                            <  2:1345075/148‑1 (MQ=255)
                                                                                                     gTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCACGCGAATCAAAGCGGCCCGTCAACCAGCGCGgt                                           <  2:2718808/148‑1 (MQ=255)
                                                                                                       cTACTTTTCGACATGGTGGCAATCCACTCAAGTTGGCACTGGAACAAACCAGCGTTTTACAATAAGCCTGCCACAAACCAACCACCAACGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtgt                                         <  1:1851888/148‑1 (MQ=255)
                                                                                                       cTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtgt                                         >  1:2497277/1‑148 (MQ=255)
                                                                                                         aCTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAACGCGGCCCGTCAACCAGCGCGGtgttg                                       <  1:2007486/148‑1 (MQ=255)
                                                                                                         aCTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttg                                       <  2:3156658/148‑1 (MQ=255)
                                                                                                          cTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttgt                                      >  2:1005272/1‑148 (MQ=255)
                                                                                                          cTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttgt                                      <  2:787432/148‑1 (MQ=255)
                                                                                                           tGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTc                                     >  2:2335748/1‑148 (MQ=255)
                                                                                                              tCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCa                                  >  1:2421654/1‑148 (MQ=255)
                                                                                                                 aCAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGgc                               >  2:2495234/1‑148 (MQ=255)
                                                                                                                  cAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGgcg                              >  2:1437437/1‑148 (MQ=255)
                                                                                                                   atggTGGCAATCCACTCAATTTGGCACTGGAAGAAACCAGGGGTTTACCATACGCCGGGCACAAAGCAACCACCACCGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGgcgc                             <  1:1002281/146‑1 (MQ=255)
                                                                                                                    aGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGcgcg                            >  1:1047134/1‑148 (MQ=255)
                                                                                                                     ggTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCgg                           >  1:1136980/1‑148 (MQ=255)
                                                                                                                      gTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCggg                          >  1:3249075/1‑148 (MQ=255)
                                                                                                                         gCAATCAACTCAAGTGGGCACTGTAATAAACCAGCGTTTTACCATACGCCGGGCACAAACCAACCCCCACCGGTGGGTGGGTGACTCGTACAGCACGCTAATCAACGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGaaa                       <  1:1206607/148‑1 (MQ=255)
                                                                                                                          cAATCCACTCCAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCACCGGTTGGTCGGTGACTCGTACCGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAg                      <  1:1767439/148‑1 (MQ=255)
                                                                                                                          cAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAg                      >  1:1920558/1‑148 (MQ=255)
                                                                                                                                 cTCAAGTTGGCAGTGGAAGAAACCAGCGTTTTACCATACGCCTGGCACAACGCAACCACCAACGGTTGGGCGGTGACTCTTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAgg               <  1:2837830/148‑1 (MQ=255)
                                                                                                                                 cTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAgg               <  2:2647029/148‑1 (MQ=255)
                                                                                                                                 cTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAgg               >  1:3369038/1‑148 (MQ=255)
                                                                                                                                  tCAAGTTGCCAGTGTAAGAAACCAGGGGTTTACAATACGCCTGGCACAAACCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAAAAAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGt              <  1:84230/148‑1 (MQ=255)
                                                                                                                                   cAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTg             <  2:3393693/148‑1 (MQ=255)
                                                                                                                                     tgtgggCACTGGAAGAAACCAGCGGTTTACCATACGCCTGGCACAAAGCAACCACCACCGGTTGGTCGGTGACTCGTACCGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGAc           <  1:716109/144‑1 (MQ=255)
                                                                                                                                     aGTTGGCACTGGAAGAAACCTGCGTTTTACCATACGCCTGGCACAAAGCAACCACCAACGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGAc           <  1:2137155/148‑1 (MQ=255)
                                                                                                                                      gTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAACGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa          <  2:2101467/148‑1 (MQ=255)
                                                                                                                                      gTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa          <  2:2103565/148‑1 (MQ=255)
                                                                                                                                       ttGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAg         >  2:525452/1‑148 (MQ=255)
                                                                                                                                           cACTGTAAGAAACCAGCGTTTTACCATAAGCCTGGCACAAAGCAACCACAAACGGTTGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAcc     <  2:1262553/148‑1 (MQ=255)
                                                                                                                                           cACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAcc     <  1:2939985/148‑1 (MQ=255)
                                                                                                                                             cTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTg   <  2:394501/148‑1 (MQ=255)
                                                                                                                                             cTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTg   >  1:799717/1‑148 (MQ=255)
                                                                                                                                              tGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTgg  <  1:140307/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGG  >  CP009273/3424218‑3424507

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: