Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,916,830 C→A 100% 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,8300CA96.5% 750.4 / 16.9 229D71Y (GAT→TAT) rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base A (113/108);  total (119/110)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 229 total aligned reads displayed.

CTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGT  >  CP009273/3916683‑3916976
                                                                                                                                                   |                                                                                                                                                  
cTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGa                                                                                                                                                  <  2:2871800/150‑1 (MQ=255)
  ggTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATa                                                                                                                                                >  1:1672975/1‑150 (MQ=255)
     tCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAc                                                                                                                                             >  1:3078760/1‑150 (MQ=255)
      cAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAcc                                                                                                                                            >  2:154146/1‑150 (MQ=255)
      cAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAcc                                                                                                                                            >  1:2012800/1‑150 (MQ=255)
         tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGtt                                                                                                                                         <  1:3399585/150‑1 (MQ=255)
          atTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGANTGCCTGGCANTCCTGGTCNGGTGCCGACATAGATAANCCGTTc                                                                                                                                        >  1:716109/1‑150 (MQ=255)
          atTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTc                                                                                                                                        >  2:2759667/1‑150 (MQ=255)
              tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCt                                                                                                                                    <  1:3995702/150‑1 (MQ=255)
               ccAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCtt                                                                                                                                   <  1:3852715/150‑1 (MQ=255)
                   tttAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAg                                                                                                                               >  2:3809579/1‑150 (MQ=255)
                    ttAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGa                                                                                                                              >  2:1484658/1‑150 (MQ=255)
                    ttAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGa                                                                                                                              >  1:4084133/1‑150 (MQ=255)
                     tAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGAt                                                                                                                             >  1:724940/1‑150 (MQ=255)
                        gCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACg                                                                                                                          <  1:1667640/150‑1 (MQ=255)
                          tCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGt                                                                                                                        >  1:3910666/1‑150 (MQ=255)
                           cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTg                                                                                                                       >  1:4640967/1‑150 (MQ=255)
                            aTGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGc                                                                                                                      <  2:1711273/150‑1 (MQ=255)
                                tGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGccacc                                                                                                                  >  1:1180750/1‑150 (MQ=255)
                                  cACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCcaccac                                                                                                                >  1:4141483/1‑150 (MQ=255)
                                  cACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCcaccac                                                                                                                >  2:4419260/1‑150 (MQ=255)
                                     cTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAAt                                                                                                             <  1:976328/150‑1 (MQ=255)
                                      tGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTGGCAGATACGGTGCCACCACAATg                                                                                                            <  1:1013692/150‑1 (MQ=255)
                                       gACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGc                                                                                                           >  1:2687767/1‑150 (MQ=255)
                                          gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                        <  1:3509021/150‑1 (MQ=255)
                                           aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                                                                                       <  1:4425713/150‑1 (MQ=255)
                                           aaGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc                                                                                                       >  2:2773264/1‑150 (MQ=255)
                                             ggAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCga                                                                                                     <  1:1483669/150‑1 (MQ=255)
                                               aaaCGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCgaga                                                                                                   >  2:1899225/1‑150 (MQ=255)
                                                aaCGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGaa                                                                                                  >  2:568180/1‑150 (MQ=255)
                                                aaCGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGaa                                                                                                  >  2:3885710/1‑150 (MQ=255)
                                                   gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAatat                                                                                               >  2:2798701/1‑150 (MQ=255)
                                                   gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAatat                                                                                               >  2:3614075/1‑150 (MQ=255)
                                                     acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATgg                                                                                             >  1:3975584/1‑150 (MQ=255)
                                                         gTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTa                                                                                         <  2:2897471/150‑1 (MQ=255)
                                                         gTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACTGTGCCACCACAATGCTATCGAGAATATGGCGTa                                                                                         <  1:2226550/150‑1 (MQ=255)
                                                         gTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTa                                                                                         <  2:2459865/150‑1 (MQ=255)
                                                         gTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTa                                                                                         >  2:2545106/1‑150 (MQ=255)
                                                           ttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAcc                                                                                       <  2:121341/150‑1 (MQ=255)
                                                           ttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAcc                                                                                       >  2:4514245/1‑150 (MQ=255)
                                                           ttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGAGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAcc                                                                                       >  2:3599518/1‑150 (MQ=255)
                                                            tACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCa                                                                                      >  1:3666552/1‑150 (MQ=255)
                                                            tACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCa                                                                                      <  1:4322699/150‑1 (MQ=255)
                                                             aCCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCANCACAATGCTATCGAGAATATGGCGTACCAg                                                                                     <  1:4286759/150‑1 (MQ=255)
                                                             aCCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGAGAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAg                                                                                     >  2:2573998/1‑150 (MQ=255)
                                                               cAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCa                                                                                   >  2:3003425/1‑150 (MQ=255)
                                                                aaGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAt                                                                                  >  2:2411704/1‑150 (MQ=255)
                                                                       tCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTa                                                                           <  1:103670/150‑1 (MQ=255)
                                                                          aaCAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGa                                                                        <  1:4255844/150‑1 (MQ=255)
                                                                          aaCAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGa                                                                        <  2:3681504/150‑1 (MQ=255)
                                                                            cAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATc                                                                      >  2:3244780/1‑150 (MQ=255)
                                                                                  gAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAc                                                                >  1:3575428/1‑150 (MQ=255)
                                                                                    aaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCg                                                              >  1:548242/1‑150 (MQ=255)
                                                                                      tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                            <  1:3807621/150‑1 (MQ=255)
                                                                                      tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                            <  2:2404413/150‑1 (MQ=255)
                                                                                      tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                            <  1:1499419/150‑1 (MQ=255)
                                                                                       gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAg                                                           <  1:2924598/150‑1 (MQ=255)
                                                                                        ggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGt                                                          <  1:4533370/150‑1 (MQ=255)
                                                                                           ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAg                                                       >  1:3858062/1‑150 (MQ=255)
                                                                                              aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                    >  1:1715554/1‑150 (MQ=255)
                                                                                               ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTg                                                   <  1:4211239/150‑1 (MQ=255)
                                                                                                  cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTAc                                                >  2:3168026/1‑150 (MQ=255)
                                                                                                  cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTAc                                                >  1:186917/1‑150 (MQ=255)
                                                                                                  cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTAc                                                >  1:3712707/1‑150 (MQ=255)
                                                                                                      cGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCtt                                            <  1:1007869/150‑1 (MQ=255)
                                                                                                       gATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCttt                                           <  1:4583253/150‑1 (MQ=255)
                                                                                                       gATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCttt                                           <  1:4037282/150‑1 (MQ=255)
                                                                                                        aTAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTg                                          <  2:1441947/150‑1 (MQ=255)
                                                                                                        aTAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTg                                          <  2:4019719/150‑1 (MQ=255)
                                                                                                            agagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTcc                                      >  1:3309575/1‑150 (MQ=255)
                                                                                                            agagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTcc                                      >  1:601587/1‑150 (MQ=255)
                                                                                                              agTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAt                                    <  1:2207622/150‑1 (MQ=255)
                                                                                                              agTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAt                                    <  1:2638654/150‑1 (MQ=255)
                                                                                                                 ggAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCGATTTa                                 >  2:4511395/1‑150 (MQ=255)
                                                                                                                 ggAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTa                                 >  2:2723423/1‑150 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                >  1:2993509/1‑150 (MQ=255)
                                                                                                                      gCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcag                            <  1:499340/150‑1 (MQ=255)
                                                                                                                       ccTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcagc                           >  2:214532/1‑150 (MQ=255)
                                                                                                                        cTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcagca                          >  2:3076385/1‑150 (MQ=255)
                                                                                                                        cTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcagca                          >  2:686602/1‑150 (MQ=255)
                                                                                                                          ggCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCata                        <  1:2113198/150‑1 (MQ=255)
                                                                                                                           gCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCatat                       >  2:1863664/1‑150 (MQ=255)
                                                                                                                            cAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATAtt                      <  2:2996990/150‑1 (MQ=255)
                                                                                                                             aGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTc                     >  2:1976704/1‑150 (MQ=255)
                                                                                                                             aGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTc                     <  2:4614908/150‑1 (MQ=255)
                                                                                                                             aGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTc                     <  2:430699/150‑1 (MQ=255)
                                                                                                                             aGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTc                     >  2:2507902/1‑150 (MQ=255)
                                                                                                                               tCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCAc                   <  1:3259572/150‑1 (MQ=255)
                                                                                                                                ccTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACg                  <  2:1950909/150‑1 (MQ=255)
                                                                                                                                 cTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGt                 >  1:1569197/1‑150 (MQ=255)
                                                                                                                                 cTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCAAGt                 >  2:1425601/1‑150 (MQ=255)
                                                                                                                                    gTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAgg              >  2:3420361/1‑150 (MQ=255)
                                                                                                                                      ccGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCa            <  1:1899225/150‑1 (MQ=255)
                                                                                                                                            ccGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGc      >  1:262213/1‑150 (MQ=255)
                                                                                                                                             cGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCt     >  1:1765177/1‑150 (MQ=255)
                                                                                                                                               aCATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTgg   >  1:464689/1‑150 (MQ=255)
                                                                                                                                               aCATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTgg   >  1:75311/1‑150 (MQ=255)
                                                                                                                                                cATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGt  <  1:1926611/150‑1 (MQ=255)
                                                                                                                                                cATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGt  <  2:2691400/150‑1 (MQ=255)
                                                                                                                                                cATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGt  >  2:2732929/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: ‑