Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,818,010 C→A 92.0% 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,0100CA92.0% 1051.0 / 79.7 339H46N (CAT→AAT) trmHtRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base A (173/139);  total (186/153)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 339 total aligned reads displayed.

ATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCT  >  CP009273/3817864‑3818155
                                                                                                                                                  |                                                                                                                                                 
aTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACAt                                                                                                                                                 >  1:33119/1‑149 (MQ=255)
   aaTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAACCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAACTGaa                                                                                                                                              >  1:3419906/1‑149 (MQ=255)
   aaTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                              >  2:5136110/1‑149 (MQ=255)
   aaTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                              >  1:828527/1‑149 (MQ=255)
    aTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAg                                                                                                                                             >  1:4472002/1‑149 (MQ=255)
     tGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGt                                                                                                                                            <  1:2771740/149‑1 (MQ=255)
       ttttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTc                                                                                                                                          >  2:4971149/1‑149 (MQ=255)
       ttttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTc                                                                                                                                          <  1:3672518/149‑1 (MQ=255)
        tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCa                                                                                                                                         <  1:1131712/149‑1 (MQ=255)
            tGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGcc                                                                                                                                     >  2:2304999/1‑149 (MQ=255)
              aaCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGt                                                                                                                                   >  1:810393/1‑149 (MQ=255)
               aCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                  <  2:1098500/149‑1 (MQ=255)
                cccAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCt                                                                                                                                 <  1:200564/149‑1 (MQ=255)
                   aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGTATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGc                                                                                                                              >  1:1665708/1‑149 (MQ=255)
                   aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGc                                                                                                                              <  2:4409091/149‑1 (MQ=255)
                   aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGc                                                                                                                              <  2:4179002/149‑1 (MQ=255)
                   aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGc                                                                                                                              >  2:369126/1‑149 (MQ=255)
                      acGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTg                                                                                                                           >  1:2192034/1‑149 (MQ=255)
                         ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                        >  1:1059859/1‑149 (MQ=255)
                         ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                        >  1:4982130/1‑149 (MQ=255)
                         ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                        <  2:3629015/149‑1 (MQ=255)
                          tATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAg                                                                                                                       <  1:2279511/149‑1 (MQ=255)
                           aTGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGGTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGc                                                                                                                      <  2:2801865/149‑1 (MQ=255)
                               aCGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCa                                                                                                                  >  2:1306794/1‑149 (MQ=255)
                                    tCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgcgc                                                                                                             >  1:5056807/1‑149 (MQ=255)
                                    tCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgcgc                                                                                                             >  1:405197/1‑149 (MQ=255)
                                         gAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCAt                                                                                                        >  2:4527589/1‑149 (MQ=255)
                                          aaaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATg                                                                                                       <  2:56045/149‑1 (MQ=255)
                                             tGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCt                                                                                                    <  1:25012/149‑1 (MQ=255)
                                             tGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCt                                                                                                    >  2:3518961/1‑149 (MQ=255)
                                              gCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCtt                                                                                                   <  2:2438024/149‑1 (MQ=255)
                                              gCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCtt                                                                                                   >  2:4867138/1‑149 (MQ=255)
                                                tCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCg                                                                                                 <  1:3105945/149‑1 (MQ=255)
                                                tCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCg                                                                                                 >  2:3918083/1‑149 (MQ=255)
                                                 cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCgg                                                                                                >  1:1208655/1‑149 (MQ=255)
                                                  gCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGc                                                                                               >  2:1062677/1‑149 (MQ=255)
                                                    cAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAg                                                                                             >  1:4105747/1‑149 (MQ=255)
                                                      ggcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcg                                                                                           >  2:1543796/1‑149 (MQ=255)
                                                      ggcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAAGGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcg                                                                                           >  2:583166/1‑149 (MQ=255)
                                                         ggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcggcg                                                                                        <  2:3082599/149‑1 (MQ=255)
                                                         ggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcggcg                                                                                        >  2:1145403/1‑149 (MQ=255)
                                                           cAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCggg                                                                                      <  1:2080402/149‑1 (MQ=255)
                                                           cAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCggg                                                                                      >  2:3193122/1‑149 (MQ=255)
                                                             gCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                    >  2:4965768/1‑149 (MQ=255)
                                                               cTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagt                                                                                  <  1:2345844/149‑1 (MQ=255)
                                                               cTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCAGAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagt                                                                                  <  2:1939573/149‑1 (MQ=255)
                                                                  aTCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAAc                                                                               <  1:5179771/149‑1 (MQ=255)
                                                                   tCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACa                                                                              <  1:3955504/149‑1 (MQ=255)
                                                                    cTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAg                                                                             >  2:4086859/1‑149 (MQ=255)
                                                                      gACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCt                                                                           >  2:4011771/1‑149 (MQ=255)
                                                                       aCCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTg                                                                          <  2:2951819/149‑1 (MQ=255)
                                                                       aCCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGGACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTg                                                                          <  1:667631/149‑1 (MQ=255)
                                                                        ccGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTgg                                                                         >  1:2135657/1‑149 (MQ=255)
                                                                         cGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTggg                                                                        >  1:4533861/1‑149 (MQ=255)
                                                                            cTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTAc                                                                     <  2:2460819/149‑1 (MQ=255)
                                                                             tGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACa                                                                    <  2:2170118/149‑1 (MQ=255)
                                                                             tGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGAGCACCATGGCTTCGGCAGCGGCGGGTCGTAACAGCTGGGTACa                                                                    >  1:188178/1‑149 (MQ=255)
                                                                               cATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTCCAgg                                                                  >  1:1764793/1‑149 (MQ=255)
                                                                                aTGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGt                                                                 <  1:3239777/149‑1 (MQ=255)
                                                                                           ccACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGACTTTGGCGTAAATGAAGTTCACGCAGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACAcc                                                      <  2:2268411/149‑1 (MQ=255)
                                                                                                 aCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCa                                                >  1:3721831/1‑149 (MQ=255)
                                                                                                 aCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCa                                                <  2:3935018/149‑1 (MQ=255)
                                                                                                 aCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCa                                                >  1:2869282/1‑149 (MQ=255)
                                                                                                  ccTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAt                                               >  1:5173337/1‑149 (MQ=255)
                                                                                                    tCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTg                                             >  2:2677649/1‑149 (MQ=255)
                                                                                                     cATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTgg                                            <  1:5102009/149‑1 (MQ=255)
                                                                                                      aTAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGc                                           <  1:3944443/149‑1 (MQ=255)
                                                                                                       tAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCg                                          <  1:3153136/149‑1 (MQ=255)
                                                                                                        aaCGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGa                                         <  2:3027998/149‑1 (MQ=255)
                                                                                                          cGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGTGATg                                       >  2:2344106/1‑149 (MQ=255)
                                                                                                           gTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGc                                      >  1:84986/1‑149 (MQ=255)
                                                                                                           gTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGc                                      <  2:4606771/149‑1 (MQ=255)
                                                                                                               cTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTc                                  >  2:1665582/1‑149 (MQ=255)
                                                                                                               cTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTc                                  >  1:4311930/1‑149 (MQ=255)
                                                                                                                   gATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTc                              <  2:841325/149‑1 (MQ=255)
                                                                                                                    attattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCa                             >  1:4765138/1‑149 (MQ=255)
                                                                                                                     ttattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAt                            <  1:2977938/149‑1 (MQ=255)
                                                                                                                     ttattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAt                            >  1:2588798/1‑149 (MQ=255)
                                                                                                                      tattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtc                           <  2:1762620/149‑1 (MQ=255)
                                                                                                                        ttCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtctc                         <  2:330769/149‑1 (MQ=255)
                                                                                                                          cGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCaa                       <  2:4719474/149‑1 (MQ=255)
                                                                                                                           gTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCaaa                      <  2:1840296/149‑1 (MQ=255)
                                                                                                                             aCCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAgg                    <  1:5343/149‑1 (MQ=255)
                                                                                                                             aCCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAgg                    >  1:3546355/1‑149 (MQ=255)
                                                                                                                              ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGc                   <  2:954984/149‑1 (MQ=255)
                                                                                                                              ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGc                   <  1:1398619/149‑1 (MQ=255)
                                                                                                                                gCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCa                 <  2:4322710/149‑1 (MQ=255)
                                                                                                                                gCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCa                 >  1:4444487/1‑149 (MQ=255)
                                                                                                                                  aGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAgg               >  2:5212119/1‑149 (MQ=255)
                                                                                                                                   gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAggg              >  2:4917935/1‑149 (MQ=255)
                                                                                                                                   gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAggg              >  1:1463767/1‑149 (MQ=255)
                                                                                                                                   gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAggg              >  1:1609793/1‑149 (MQ=255)
                                                                                                                                     tGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTACCAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCa            >  2:4821412/1‑149 (MQ=255)
                                                                                                                                     tGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCa            <  2:1827065/149‑1 (MQ=255)
                                                                                                                                      gCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCAt           <  1:4645849/149‑1 (MQ=255)
                                                                                                                                        cGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGc         <  2:376773/149‑1 (MQ=255)
                                                                                                                                            ggCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt     <  1:1840138/149‑1 (MQ=255)
                                                                                                                                              cGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTc   <  2:341983/149‑1 (MQ=255)
                                                                                                                                               gTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCt  >  2:2842293/1‑149 (MQ=255)
                                                                                                                                               gTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCt  <  2:689476/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: