Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,818,010 C→A 100% H46N (CAT→AAT)  trmH → tRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,818,0100CA99.6% 912.8 / ‑4.4 272H46N (CAT→AAT) trmHtRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (142/128);  total (143/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 272 total aligned reads displayed.

ATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCT  >  CP009273/3817864‑3818155
                                                                                                                                                  |                                                                                                                                                 
aTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAAt                                                                                                                                                 <  1:766093/149‑1 (MQ=255)
 ttAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATg                                                                                                                                                >  1:128816/1‑149 (MQ=255)
 ttAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATg                                                                                                                                                >  1:3990160/1‑149 (MQ=255)
 ttAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATg                                                                                                                                                >  1:1895806/1‑149 (MQ=255)
   aaTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                              <  2:2187928/149‑1 (MQ=255)
     tGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGt                                                                                                                                            >  2:3312324/1‑149 (MQ=255)
      gTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGtt                                                                                                                                           >  1:1666570/1‑149 (MQ=255)
        tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGTGTAAATGAAGTTCa                                                                                                                                         <  1:892402/149‑1 (MQ=255)
        tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCa                                                                                                                                         <  2:258022/149‑1 (MQ=255)
        tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAAACTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGTGTAAATGAAGTTCa                                                                                                                                         <  1:3282126/149‑1 (MQ=255)
         ttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTGCCGCAGATGCCGTTGGCGTAAATGAAGTTCAc                                                                                                                                        <  2:4034496/149‑1 (MQ=255)
          tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACg                                                                                                                                       >  1:855414/1‑149 (MQ=255)
          tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACg                                                                                                                                       >  1:856014/1‑149 (MQ=255)
           aTGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGc                                                                                                                                      <  1:2380047/149‑1 (MQ=255)
           aTGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGc                                                                                                                                      <  1:2598144/149‑1 (MQ=255)
              aaCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGt                                                                                                                                   <  2:2917085/149‑1 (MQ=255)
               aCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                  >  1:643843/1‑149 (MQ=255)
                 ccAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTg                                                                                                                                >  2:592730/1‑149 (MQ=255)
                  cAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTgg                                                                                                                               <  1:3538829/149‑1 (MQ=255)
                   aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGc                                                                                                                              <  1:4322779/149‑1 (MQ=255)
                     cacGTTAGGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCt                                                                                                                            <  2:275293/149‑1 (MQ=255)
                      acGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTg                                                                                                                           >  1:2633651/1‑149 (MQ=255)
                       cGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTgg                                                                                                                          <  1:4356979/149‑1 (MQ=255)
                            tGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGcc                                                                                                                     <  2:1423070/149‑1 (MQ=255)
                            tGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGcc                                                                                                                     >  2:4500027/1‑149 (MQ=255)
                            tGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGcc                                                                                                                     >  1:1775658/1‑149 (MQ=255)
                               ncGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCATGCCGTCTGGCCTGGTAGCCGCa                                                                                                                  >  1:4603505/2‑149 (MQ=255)
                                cGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCAt                                                                                                                 <  2:4778034/149‑1 (MQ=255)
                                cGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCAt                                                                                                                 >  1:1142422/1‑149 (MQ=255)
                                 gCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATg                                                                                                                >  1:581813/1‑149 (MQ=255)
                                  cATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgc                                                                                                               >  1:2111639/1‑149 (MQ=255)
                                   aTCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgcg                                                                                                              >  1:1567441/1‑149 (MQ=255)
                                       gCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCAcc                                                                                                          <  2:4329764/149‑1 (MQ=255)
                                         gAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCAt                                                                                                        <  2:4611831/149‑1 (MQ=255)
                                         gAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCAt                                                                                                        <  1:3956357/149‑1 (MQ=255)
                                          aaaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATg                                                                                                       >  2:3048454/1‑149 (MQ=255)
                                          aaaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATg                                                                                                       >  1:1978683/1‑149 (MQ=255)
                                          aaaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATg                                                                                                       >  2:2881485/1‑149 (MQ=255)
                                           aaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATgg                                                                                                      <  1:1094426/149‑1 (MQ=255)
                                           aaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATgg                                                                                                      <  2:495505/149‑1 (MQ=255)
                                               cTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTc                                                                                                  >  2:1526018/1‑149 (MQ=255)
                                                tCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCg                                                                                                 >  2:1591511/1‑149 (MQ=255)
                                                 cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCgg                                                                                                <  2:456193/149‑1 (MQ=255)
                                                   ccAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCa                                                                                              <  1:3798655/149‑1 (MQ=255)
                                                      ggcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcg                                                                                           >  2:1971833/1‑149 (MQ=255)
                                                            aGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgt                                                                                     <  2:4450127/149‑1 (MQ=255)
                                                             gCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                    >  2:1321762/1‑149 (MQ=255)
                                                              ccTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagt                                                                                  <  1:2506824/149‑1 (MQ=255)
                                                               cTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagt                                                                                  <  1:94176/149‑1 (MQ=255)
                                                                 gATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTaa                                                                                >  1:3366089/1‑149 (MQ=255)
                                                                  aTCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAAc                                                                               <  2:3618556/149‑1 (MQ=255)
                                                                        ccGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTgg                                                                         <  2:33221/149‑1 (MQ=255)
                                                                         cGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTggg                                                                        <  2:1705537/149‑1 (MQ=255)
                                                                            cTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTAc                                                                     <  2:3517678/149‑1 (MQ=255)
                                                                             tGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACa                                                                    <  2:652984/149‑1 (MQ=255)
                                                                                 tGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTCCAGGTg                                                                >  2:2503294/1‑149 (MQ=255)
                                                                                 tGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTg                                                                <  1:221363/149‑1 (MQ=255)
                                                                                 tGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTg                                                                <  2:3758835/149‑1 (MQ=255)
                                                                                    aGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGaaa                                                             <  1:3036473/149‑1 (MQ=255)
                                                                                      cAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAac                                                           <  2:1124681/149‑1 (MQ=255)
                                                                                       aGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAaca                                                          >  2:272593/1‑149 (MQ=255)
                                                                                         gTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacaca                                                        <  2:3260270/149‑1 (MQ=255)
                                                                                              caAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCa                                                   <  2:2742555/149‑1 (MQ=255)
                                                                                              caAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCa                                                   >  2:1952958/1‑149 (MQ=255)
                                                                                              caAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCa                                                   >  2:434327/1‑149 (MQ=255)
                                                                                                aaCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAcc                                                 >  2:868805/1‑149 (MQ=255)
                                                                                                aaCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAcc                                                 >  1:682829/1‑149 (MQ=255)
                                                                                                aaCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAcc                                                 >  1:3600030/1‑149 (MQ=255)
                                                                                                  ccTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAt                                               >  2:2585765/1‑149 (MQ=255)
                                                                                                  ccTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAt                                               <  1:1527784/149‑1 (MQ=255)
                                                                                                     cATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTgg                                            <  1:3451276/149‑1 (MQ=255)
                                                                                                         aCGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGAt                                        >  1:1763735/1‑149 (MQ=255)
                                                                                                          cGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATg                                       <  2:3607470/149‑1 (MQ=255)
                                                                                                           gTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGc                                      <  1:3434135/149‑1 (MQ=255)
                                                                                                             ttCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCg                                    >  1:4092825/1‑149 (MQ=255)
                                                                                                              tCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGt                                   >  2:2335202/1‑149 (MQ=255)
                                                                                                               cTGCGAGTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTc                                  >  1:3697562/1‑149 (MQ=255)
                                                                                                                 gCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGc                                >  2:4715242/1‑149 (MQ=255)
                                                                                                                    attattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCa                             >  2:792546/1‑149 (MQ=255)
                                                                                                                      tattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtc                           <  2:1740710/149‑1 (MQ=255)
                                                                                                                       attCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtct                          <  1:1024754/149‑1 (MQ=255)
                                                                                                                        ttCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtctc                         >  2:3368225/1‑149 (MQ=255)
                                                                                                                             aCCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAgg                    <  2:2050954/149‑1 (MQ=255)
                                                                                                                              ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGc                   <  2:3688953/149‑1 (MQ=255)
                                                                                                                              ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGc                   <  2:4187915/149‑1 (MQ=255)
                                                                                                                                 cAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAg                <  1:4455730/149‑1 (MQ=255)
                                                                                                                                 cAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAg                <  2:4722007/149‑1 (MQ=255)
                                                                                                                                  aGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAgg               >  1:4284984/1‑149 (MQ=255)
                                                                                                                                  aGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAgg               >  1:2725957/1‑149 (MQ=255)
                                                                                                                                  aGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAgg               <  2:1367097/149‑1 (MQ=255)
                                                                                                                                   gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAggg              >  1:3262970/1‑149 (MQ=255)
                                                                                                                                     tGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTCCAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCGCATCTCAAAGGCCAGGGCa            >  2:3985512/1‑149 (MQ=255)
                                                                                                                                      gCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCAt           <  1:4745148/149‑1 (MQ=255)
                                                                                                                                      gCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGAGATGCCGTCGCTCATCTCAAAGGCCAGGGCAt           >  2:4645190/1‑149 (MQ=255)
                                                                                                                                          ttGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAg       >  1:3334768/1‑149 (MQ=255)
                                                                                                                                            ggCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt     >  1:137814/1‑149 (MQ=255)
                                                                                                                                            ggCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt     <  1:4500357/149‑1 (MQ=255)
                                                                                                                                              cGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTc   <  1:989525/149‑1 (MQ=255)
                                                                                                                                               gTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCt  <  1:3175004/149‑1 (MQ=255)
                                                                                                                                               gTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCt  <  1:2906102/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCT  >  CP009273/3817864‑3818155

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: