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.2% 855.4 / ‑1.6 248H46N (CAT→AAT) trmHtRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (148/98);  total (148/100)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-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 248 total aligned reads displayed.

AATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTG  >  CP009273/3817863‑3818156
                                                                                                                                                   |                                                                                                                                                  
aaTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAAt                                                                                                                                                  >  1:656084/1‑150 (MQ=255)
   tAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                               <  1:3779902/150‑1 (MQ=255)
    aaTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAg                                                                                                                                              <  1:2297807/150‑1 (MQ=255)
      tGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGtt                                                                                                                                            <  1:632697/150‑1 (MQ=255)
        ttttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCa                                                                                                                                          >  2:2698385/1‑150 (MQ=255)
        ttttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCa                                                                                                                                          <  1:1950949/150‑1 (MQ=255)
           tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGc                                                                                                                                       >  1:2422304/1‑150 (MQ=255)
           tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGc                                                                                                                                       >  2:1239523/1‑150 (MQ=255)
           tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGc                                                                                                                                       >  2:3737283/1‑150 (MQ=255)
             tGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCg                                                                                                                                     >  2:3116497/1‑150 (MQ=255)
             tGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCg                                                                                                                                     >  1:2928338/1‑150 (MQ=255)
               aaCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                   >  2:1295300/1‑150 (MQ=255)
               aaCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                   >  2:3852373/1‑150 (MQ=255)
               aaCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTATCTGCGATTATTCGAACCGCAGATGCCGTAGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                   >  1:2051877/1‑150 (MQ=255)
                  ccAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTgg                                                                                                                                >  2:4085306/1‑150 (MQ=255)
                    aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGcc                                                                                                                              >  2:245927/1‑150 (MQ=255)
                    aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGAGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGcc                                                                                                                              >  1:3124505/1‑150 (MQ=255)
                     acacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCt                                                                                                                             >  2:3918745/1‑150 (MQ=255)
                      cacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTg                                                                                                                            <  2:627686/150‑1 (MQ=255)
                         gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                         >  2:3892384/1‑150 (MQ=255)
                         gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                         <  2:1566396/150‑1 (MQ=255)
                            aTGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGcc                                                                                                                      <  1:910305/150‑1 (MQ=255)
                            aTGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGcc                                                                                                                      >  1:2694645/1‑150 (MQ=255)
                               cACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCa                                                                                                                   >  2:824769/1‑150 (MQ=255)
                               cACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCa                                                                                                                   >  1:258234/1‑150 (MQ=255)
                                aCGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCAt                                                                                                                  >  2:3375440/1‑150 (MQ=255)
                                          gAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTATCCGCATGCGCACCATg                                                                                                        <  2:3764736/150‑1 (MQ=255)
                                          gAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATg                                                                                                        >  2:2340784/1‑150 (MQ=255)
                                           aaaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATgg                                                                                                       >  2:1662341/1‑150 (MQ=255)
                                            aaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGc                                                                                                      <  2:2782491/150‑1 (MQ=255)
                                            aaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGc                                                                                                      >  1:862618/1‑150 (MQ=255)
                                             aTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCt                                                                                                     >  2:447806/1‑150 (MQ=255)
                                              tGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCtt                                                                                                    >  2:1185319/1‑150 (MQ=255)
                                                  cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGc                                                                                                >  1:3333522/1‑150 (MQ=255)
                                                      aGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcg                                                                                            >  1:2766424/1‑150 (MQ=255)
                                                      aGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcg                                                                                            >  1:1506404/1‑150 (MQ=255)
                                                        gcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcggc                                                                                          <  2:1824628/150‑1 (MQ=255)
                                                             aGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                     <  2:1364020/150‑1 (MQ=255)
                                                             aGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                     >  2:3051768/1‑150 (MQ=255)
                                                             aGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                     >  1:1863103/1‑150 (MQ=255)
                                                              gCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtag                                                                                    <  1:2113676/150‑1 (MQ=255)
                                                                cTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagta                                                                                  <  2:2488950/150‑1 (MQ=255)
                                                                cTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagta                                                                                  >  2:2514414/1‑150 (MQ=255)
                                                                   aTCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACa                                                                               >  1:2872326/1‑150 (MQ=255)
                                                                   aTCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACa                                                                               <  2:1058620/150‑1 (MQ=255)
                                                                   aTCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACa                                                                               >  1:2492619/1‑150 (MQ=255)
                                                                    tCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATCACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAg                                                                              <  2:36219/150‑1 (MQ=255)
                                                                     cTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGc                                                                             >  2:478643/1‑150 (MQ=255)
                                                                     cTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGc                                                                             >  2:1538866/1‑150 (MQ=255)
                                                                           gTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTa                                                                       <  1:2006473/150‑1 (MQ=255)
                                                                              tGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAg                                                                    >  1:1315047/1‑150 (MQ=255)
                                                                               gCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAgg                                                                   <  1:4103387/150‑1 (MQ=255)
                                                                                cATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGt                                                                  <  2:629735/150‑1 (MQ=255)
                                                                                cATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGt                                                                  >  1:3799960/1‑150 (MQ=255)
                                                                                     aGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGaaaa                                                             <  2:789554/150‑1 (MQ=255)
                                                                                        aGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacac                                                          >  1:1999415/1‑150 (MQ=255)
                                                                                          gTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacacac                                                        <  2:2981609/150‑1 (MQ=255)
                                                                                          gTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacacac                                                        >  1:1242101/1‑150 (MQ=255)
                                                                                           tCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACAcc                                                       <  2:427805/150‑1 (MQ=255)
                                                                                                aaaCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAcc                                                  >  2:1701260/1‑150 (MQ=255)
                                                                                                aaaCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAcc                                                  >  1:769119/1‑150 (MQ=255)
                                                                                                aaaCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAcc                                                  >  1:2635350/1‑150 (MQ=255)
                                                                                                 aaCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCa                                                 <  1:1008997/150‑1 (MQ=255)
                                                                                                  aCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAt                                                <  1:1403974/150‑1 (MQ=255)
                                                                                                   ccTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAtt                                               <  2:2123698/150‑1 (MQ=255)
                                                                                                    cTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTg                                              <  2:3182115/150‑1 (MQ=255)
                                                                                                     tCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTgg                                             <  1:1485720/150‑1 (MQ=255)
                                                                                                      cATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGc                                            >  2:2222437/1‑150 (MQ=255)
                                                                                                       aTAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCg                                           >  1:34044/1‑150 (MQ=255)
                                                                                                         aaCGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGAt                                         >  2:1007830/1‑150 (MQ=255)
                                                                                                           cGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGc                                       <  1:172970/150‑1 (MQ=255)
                                                                                                              ttCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGt                                    <  1:1652132/150‑1 (MQ=255)
                                                                                                               tCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTc                                   >  1:3887159/1‑150 (MQ=255)
                                                                                                               tCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTc                                   <  2:3688361/150‑1 (MQ=255)
                                                                                                                  gCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCt                                <  2:2003112/150‑1 (MQ=255)
                                                                                                                    gATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCa                              >  2:226293/1‑150 (MQ=255)
                                                                                                                    gATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCa                              <  2:1302846/150‑1 (MQ=255)
                                                                                                                    gATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCa                              >  1:523010/1‑150 (MQ=255)
                                                                                                                     attattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAt                             <  1:2741249/150‑1 (MQ=255)
                                                                                                                           cGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCaaa                       >  1:454847/1‑150 (MQ=255)
                                                                                                                            gTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAg                      <  2:2905988/150‑1 (MQ=255)
                                                                                                                             tACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAgg                     >  2:3071685/1‑150 (MQ=255)
                                                                                                                             tACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAgg                     <  2:2602548/150‑1 (MQ=255)
                                                                                                                              aCCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGc                    <  2:2936324/150‑1 (MQ=255)
                                                                                                                              aCCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGc                    >  1:3550705/1‑150 (MQ=255)
                                                                                                                               ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                   <  2:3874153/150‑1 (MQ=255)
                                                                                                                                cGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCa                  >  2:1097377/1‑150 (MQ=255)
                                                                                                                                cGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCa                  >  2:1097367/1‑150 (MQ=255)
                                                                                                                                 gCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAg                 >  1:750372/1‑150 (MQ=255)
                                                                                                                                    gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGc              >  1:2961921/1‑150 (MQ=255)
                                                                                                                                    gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGc              >  1:3919352/1‑150 (MQ=255)
                                                                                                                                          gTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAGTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAg        >  2:2354567/1‑150 (MQ=255)
                                                                                                                                           ttGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGa       >  1:1620890/1‑150 (MQ=255)
                                                                                                                                           ttGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGa       >  1:2990753/1‑150 (MQ=255)
                                                                                                                                            tGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt      <  1:2029702/150‑1 (MQ=255)
                                                                                                                                            tGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt      >  1:4096862/1‑150 (MQ=255)
                                                                                                                                             ggCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAtt     <  2:2282356/150‑1 (MQ=255)
                                                                                                                                              gCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTc    <  2:2559413/150‑1 (MQ=255)
                                                                                                                                               cGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCt   <  1:3354683/150‑1 (MQ=255)
                                                                                                                                                gTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTg  >  1:31986/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTG  >  CP009273/3817863‑3818156

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: