Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,916,830 C→A 89.9% D71Y (GAT→TAT)  rsmG ← 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,916,8300CA89.9% 784.5 / 78.3 286D71Y (GAT→TAT) rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base C (12/17);  new base A (117/139);  total (129/157)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 286 total aligned reads displayed.

CTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGT  >  CP009273/3916683‑3916976
                                                                                                                                                   |                                                                                                                                                  
cTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGa                                                                                                                                                  <  2:1975583/150‑1 (MQ=255)
    ttCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATaaa                                                                                                                                              >  2:458381/1‑150 (MQ=255)
      cAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAcc                                                                                                                                            <  2:2502479/150‑1 (MQ=255)
       aaTATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCCGTCCTGGTCCGGTGCCGACATAGATAAACCg                                                                                                                                           >  2:3919062/1‑150 (MQ=255)
       aaTATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCg                                                                                                                                           >  1:2645562/1‑150 (MQ=255)
         tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGtt                                                                                                                                         >  2:3505947/1‑150 (MQ=255)
            tccccAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAc                                                                                                                                      >  1:608459/4‑150 (MQ=255)
               ccAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCtt                                                                                                                                   <  1:3593609/150‑1 (MQ=255)
               ccAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCtt                                                                                                                                   <  1:3595056/150‑1 (MQ=255)
                cAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTg                                                                                                                                  <  1:2536193/150‑1 (MQ=255)
                 aGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGc                                                                                                                                 >  1:3334005/1‑150 (MQ=255)
                      aaGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATa                                                                                                                            >  1:2654982/1‑150 (MQ=255)
                       aGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATAc                                                                                                                           >  1:2065538/1‑150 (MQ=255)
                           cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTg                                                                                                                       >  2:205266/1‑150 (MQ=255)
                            aTGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGc                                                                                                                      <  2:3058662/150‑1 (MQ=255)
                             ttttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGcc                                                                                                                     >  1:382762/3‑150 (MQ=255)
                             ttttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGcc                                                                                                                     >  1:383712/3‑150 (MQ=255)
                                   aCCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACcaca                                                                                                               <  2:2586697/150‑1 (MQ=255)
                                    ccTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACaa                                                                                                              <  1:873303/150‑1 (MQ=255)
                                        aCGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCt                                                                                                          <  2:187780/150‑1 (MQ=255)
                                         cGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGcca                                                                                                                    >  1:2051023/1‑139 (MQ=255)
                                         cGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTa                                                                                                         <  1:979517/150‑1 (MQ=255)
                                          gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                        >  1:1011614/1‑150 (MQ=255)
                                           aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                                                                                       <  1:1050303/150‑1 (MQ=255)
                                           aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                                                                                       <  2:1521781/150‑1 (MQ=255)
                                           aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                                                                                       <  1:4012319/150‑1 (MQ=255)
                                           aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                                                                                       <  2:1699129/150‑1 (MQ=255)
                                              gAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCgag                                                                                                    <  2:2766260/150‑1 (MQ=255)
                                              gAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCgag                                                                                                    <  1:914445/150‑1 (MQ=255)
                                                aaCGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGaa                                                                                                  >  1:3917435/1‑150 (MQ=255)
                                                  cgcacgcGTTTACCAAGGCTATCCTACAGTGTGTAATGGGCTTCATGACGCACGATAGAGAGTGGAATGCCTGTCATTCATGGTACTGTGCCGACATTGATAAACCGTTATCCTTTCAGATACGTTGCCACCACAATGCTATCGAGAata                                                                                                <  1:708880/150‑1 (MQ=255)
                                                   gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAatat                                                                                               <  1:2360580/150‑1 (MQ=255)
                                                   gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAatat                                                                                               <  2:4147063/150‑1 (MQ=255)
                                                   gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAatat                                                                                               <  1:915053/150‑1 (MQ=255)
                                                      cgcgTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGc                                                                                            <  1:2921968/150‑1 (MQ=255)
                                                       gcgTTTACCAAGGCTATCCAACAGAGTTAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCg                                                                                           >  1:4264852/1‑150 (MQ=255)
                                                         gTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTa                                                                                         <  2:1173066/150‑1 (MQ=255)
                                                         gTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTa                                                                                         >  2:4512117/1‑150 (MQ=255)
                                                          tttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAc                                                                                        <  1:1395109/150‑1 (MQ=255)
                                                           ttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAcc                                                                                       <  2:24575/150‑1 (MQ=255)
                                                             aCCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAg                                                                                     >  2:2977913/1‑150 (MQ=255)
                                                               cAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCa                                                                                   <  1:700528/150‑1 (MQ=255)
                                                               cAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCa                                                                                   <  1:1296964/150‑1 (MQ=255)
                                                                aaGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAt                                                                                  >  2:109121/1‑150 (MQ=255)
                                                                  ggCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAtct                                                                                >  1:2767754/1‑150 (MQ=255)
                                                                     tATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAt                                                                             >  2:3352690/1‑150 (MQ=255)
                                                                             agagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcg                                                                     >  1:637788/1‑150 (MQ=255)
                                                                               agTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTTCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATAAGGATcgcg                                                                   >  2:340427/1‑150 (MQ=255)
                                                                                 tGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGa                                                                 <  1:3962089/150‑1 (MQ=255)
                                                                                  gAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAc                                                                >  1:2207677/1‑150 (MQ=255)
                                                                                  gAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAc                                                                >  1:4445447/1‑150 (MQ=255)
                                                                                   aaaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAcc                                                               <  1:1134248/150‑1 (MQ=255)
                                                                                    aaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTCGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGAGCCACCACAATGCTATCGAGAATAAGTCGTACCAGAATCTCATTAGGATCGCGGACCg                                                              >  2:1963592/1‑150 (MQ=255)
                                                                                      tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                            >  1:2139343/1‑150 (MQ=255)
                                                                                       gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGAGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAg                                                           >  1:1817039/1‑150 (MQ=255)
                                                                                          cTTCAGGACGCACGATAGAGAGTGGTATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCa                                                        <  2:744408/150‑1 (MQ=255)
                                                                                           ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAg                                                       >  2:674549/1‑150 (MQ=255)
                                                                                           ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAg                                                       <  2:842028/150‑1 (MQ=255)
                                                                                             cAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGt                                                     <  2:216273/150‑1 (MQ=255)
                                                                                              aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                    >  2:3611908/1‑150 (MQ=255)
                                                                                              aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                    >  2:2904778/1‑150 (MQ=255)
                                                                                              aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                    >  2:2001047/1‑150 (MQ=255)
                                                                                              aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                    <  2:3806707/150‑1 (MQ=255)
                                                                                               ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTg                                                   >  1:767200/1‑150 (MQ=255)
                                                                                               ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTg                                                   >  2:4251488/1‑150 (MQ=255)
                                                                                               ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTg                                                   <  1:3112754/150‑1 (MQ=255)
                                                                                                  cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTAc                                                >  1:1351721/1‑150 (MQ=255)
                                                                                                   gCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACg                                               >  2:4334420/1‑150 (MQ=255)
                                                                                                     aCGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTTCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGATTATGGCGTACCAGATTCTCATTAGGATCGCGGACCGATGTCAGGTTGTACGCt                                             <  2:2683523/150‑1 (MQ=255)
                                                                                                     aCGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCt                                             <  2:1859570/150‑1 (MQ=255)
                                                                                                     aCGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCt                                             <  1:458486/150‑1 (MQ=255)
                                                                                                      cGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCtt                                            >  2:1779640/1‑150 (MQ=255)
                                                                                                       gATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCCCCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCttt                                           <  1:3695252/150‑1 (MQ=255)
                                                                                                        aTAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTg                                          >  2:2219658/1‑150 (MQ=255)
                                                                                                          agagagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGtt                                        >  2:2949279/1‑150 (MQ=255)
                                                                                                               gTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAtt                                   <  2:3804868/150‑1 (MQ=255)
                                                                                                                tGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAttt                                  >  2:708062/1‑150 (MQ=255)
                                                                                                                tGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAttt                                  <  2:2653681/150‑1 (MQ=255)
                                                                                                                tGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGAATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAttt                                  >  2:1662249/1‑150 (MQ=255)
                                                                                                                tGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAttt                                  <  2:1223090/150‑1 (MQ=255)
                                                                                                                 ggAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTa                                 >  1:275314/1‑150 (MQ=255)
                                                                                                                 ggAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGGGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTa                                 <  2:3791758/150‑1 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                <  1:1714742/150‑1 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                <  1:844102/150‑1 (MQ=255)
                                                                                                                   aaTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATg                               <  2:3050717/150‑1 (MQ=255)
                                                                                                                   aaTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATg                               <  2:3049455/150‑1 (MQ=255)
                                                                                                                   aaTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATg                               <  1:287993/150‑1 (MQ=255)
                                                                                                                       ccTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcagc                           <  1:4381813/150‑1 (MQ=255)
                                                                                                                         tGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCat                         <  1:2899307/150‑1 (MQ=255)
                                                                                                                            cAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATAtt                      <  1:4334420/150‑1 (MQ=255)
                                                                                                                             aGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTc                     <  2:1940221/150‑1 (MQ=255)
                                                                                                                             aGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTc                     >  2:950830/1‑150 (MQ=255)
                                                                                                                               tCCTTGTCCGGTGCAGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCAc                   <  2:2789656/150‑1 (MQ=255)
                                                                                                                                   ggTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAg               <  2:4458055/150‑1 (MQ=255)
                                                                                                                                       cGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTGTGTTCCTGTTATGCAGCATAGTCACGTAGGCac           >  2:1139076/1‑149 (MQ=255)
                                                                                                                                        ggTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATATCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAAt          >  2:3176538/1‑150 (MQ=255)
                                                                                                                                              gACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTg    <  2:1302407/150‑1 (MQ=255)
                                                                                                                                              gACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTg    >  1:2570362/1‑150 (MQ=255)
                                                                                                                                              gACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTg    >  1:1565892/1‑150 (MQ=255)
                                                                                                                                               aCATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTgg   >  2:3524612/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGT  >  CP009273/3916683‑3916976

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: