Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,916,830 C→A 93.8% 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,8300CA93.8% 819.1 / 40.7 260D71Y (GAT→TAT) rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base A (109/135);  total (117/143)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 260 total aligned reads displayed.

CTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTT  >  CP009273/3916683‑3916977
                                                                                                                                                   |                                                                                                                                                   
cTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGa                                                                                                                                                   >  2:257893/1‑150 (MQ=255)
 tGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGAt                                                                                                                                                  >  1:3986344/1‑150 (MQ=255)
  ggTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATa                                                                                                                                                 >  2:2856465/1‑150 (MQ=255)
    ttCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATaaa                                                                                                                                               <  2:582985/150‑1 (MQ=255)
       aaTATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCg                                                                                                                                            >  2:134796/1‑150 (MQ=255)
        atatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGt                                                                                                                                           >  2:242975/1‑150 (MQ=255)
         tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGtt                                                                                                                                          >  2:3985987/1‑150 (MQ=255)
            tctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAc                                                                                                                                       >  2:2984091/1‑150 (MQ=255)
            tctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAc                                                                                                                                       >  2:5063205/1‑150 (MQ=255)
               ccAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCtt                                                                                                                                    <  2:1992964/150‑1 (MQ=255)
                cAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTg                                                                                                                                   <  1:1928264/150‑1 (MQ=255)
                  gTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCa                                                                                                                                 <  2:1714948/150‑1 (MQ=255)
                  gTTTAAGCTCATGGTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCa                                                                                                                                 <  2:283/150‑1 (MQ=255)
                   tttAAGCTCATGGTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAg                                                                                                                                <  2:4908074/150‑1 (MQ=255)
                      aaGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAGATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATa                                                                                                                             <  2:824588/150‑1 (MQ=255)
                      aaGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATa                                                                                                                             >  2:4323273/1‑150 (MQ=255)
                      aaGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATa                                                                                                                             <  2:1894407/150‑1 (MQ=255)
                       aGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATAc                                                                                                                            <  2:4302506/150‑1 (MQ=255)
                        gCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACg                                                                                                                           >  1:4442808/1‑150 (MQ=255)
                        gCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACg                                                                                                                           <  1:1207442/150‑1 (MQ=255)
                          tCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGt                                                                                                                         >  2:4986444/1‑150 (MQ=255)
                          tCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGt                                                                                                                         >  1:2373240/1‑150 (MQ=255)
                              gTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGcca                                                                                                                     <  1:2535473/150‑1 (MQ=255)
                                tGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGccacc                                                                                                                   <  2:4990824/150‑1 (MQ=255)
                                tGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGccacc                                                                                                                   <  2:2023795/150‑1 (MQ=255)
                                     cTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAAt                                                                                                              >  1:4574284/1‑150 (MQ=255)
                                     cTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAAt                                                                                                              >  1:5005012/1‑150 (MQ=255)
                                      tGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATg                                                                                                             >  1:3393348/1‑150 (MQ=255)
                                          gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                         >  1:2581343/1‑150 (MQ=255)
                                          gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                         >  1:1822196/1‑150 (MQ=255)
                                             ggAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCga                                                                                                      <  2:3259152/150‑1 (MQ=255)
                                             ggAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCga                                                                                                      <  1:944658/150‑1 (MQ=255)
                                             ggAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCga                                                                                                      >  1:712099/1‑150 (MQ=255)
                                                    cacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATg                                                                                               >  1:4396374/1‑150 (MQ=255)
                                                      cgcgTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGc                                                                                             >  1:3266678/1‑150 (MQ=255)
                                                        cgTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGt                                                                                           >  1:2351038/1‑150 (MQ=255)
                                                          tttACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTAc                                                                                         <  1:2550599/150‑1 (MQ=255)
                                                            tACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCa                                                                                       >  2:3735800/1‑150 (MQ=255)
                                                            tACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCa                                                                                       <  1:154778/150‑1 (MQ=255)
                                                             aCCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAg                                                                                      <  1:1976181/150‑1 (MQ=255)
                                                              ccAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGc                                                                                     <  1:1091772/150‑1 (MQ=255)
                                                                aaGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAt                                                                                   >  2:3161310/1‑150 (MQ=255)
                                                                 aGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAtc                                                                                  >  1:4526964/1‑150 (MQ=255)
                                                                     tATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAt                                                                              <  1:2793468/150‑1 (MQ=255)
                                                                      aTCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAtt                                                                             <  1:2072062/150‑1 (MQ=255)
                                                                        ccAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAg                                                                           <  2:896111/150‑1 (MQ=255)
                                                                         cAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAgg                                                                          >  2:685153/1‑150 (MQ=255)
                                                                          aaCAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGa                                                                         >  2:3089294/1‑150 (MQ=255)
                                                                              gagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcgc                                                                     <  2:1193259/150‑1 (MQ=255)
                                                                                gTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCgg                                                                   <  1:2506843/150‑1 (MQ=255)
                                                                                 tGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGa                                                                  >  2:3683995/1‑150 (MQ=255)
                                                                                 tGAAATGGGCGTCAGGACGCACGAGAGAGAGTGGAAGGCCTGGCAGTCCGGGTCCGGTGCCGACATAGATACACCGTTCACCGTGCAGATACGGTGCCACCACCATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGa                                                                  <  2:4604982/150‑1 (MQ=255)
                                                                                   aaaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAcc                                                                <  1:4446481/150‑1 (MQ=255)
                                                                                      tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                             <  1:2939145/150‑1 (MQ=255)
                                                                                      tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGaa                                                             >  1:4428223/1‑150 (MQ=255)
                                                                                       gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAg                                                            >  2:3483663/1‑150 (MQ=255)
                                                                                          cTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCa                                                         >  2:796705/1‑150 (MQ=255)
                                                                                          cTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCa                                                         <  2:4724276/150‑1 (MQ=255)
                                                                                          cTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCa                                                         <  2:1818735/150‑1 (MQ=255)
                                                                                           ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAg                                                        >  1:1216486/1‑150 (MQ=255)
                                                                                           ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAg                                                        >  2:2101406/1‑150 (MQ=255)
                                                                                             cAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGt                                                      <  1:2183455/150‑1 (MQ=255)
                                                                                             cAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGt                                                      <  2:4122733/150‑1 (MQ=255)
                                                                                             cAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGt                                                      >  1:4019705/1‑150 (MQ=255)
                                                                                              aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                     >  2:4803965/1‑150 (MQ=255)
                                                                                               ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTg                                                    >  2:4559410/1‑150 (MQ=255)
                                                                                               ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTg                                                    <  1:4842473/150‑1 (MQ=255)
                                                                                                gACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGt                                                   <  1:3463305/150‑1 (MQ=255)
                                                                                                 aCGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTa                                                  <  1:2103253/150‑1 (MQ=255)
                                                                                                  cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTAc                                                 <  1:3013363/150‑1 (MQ=255)
                                                                                                  cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTAc                                                 <  2:4376568/150‑1 (MQ=255)
                                                                                                   gCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACg                                                >  2:3123349/1‑150 (MQ=255)
                                                                                                   gCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACg                                                <  1:3403843/150‑1 (MQ=255)
                                                                                                       gATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCttt                                            <  1:1454910/150‑1 (MQ=255)
                                                                                                       gATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCttt                                            <  1:1253569/150‑1 (MQ=255)
                                                                                                        aTAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTg                                           <  1:353079/150‑1 (MQ=255)
                                                                                                         tAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGt                                          <  2:1341968/150‑1 (MQ=255)
                                                                                                          agagagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACAAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGtt                                         <  1:1381665/150‑1 (MQ=255)
                                                                                                           gagagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTc                                        <  2:3357797/150‑1 (MQ=255)
                                                                                                             gagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCa                                      <  2:4035057/150‑1 (MQ=255)
                                                                                                               gTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAtt                                    <  1:4284629/150‑1 (MQ=255)
                                                                                                                 ggAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTa                                  <  1:1877299/150‑1 (MQ=255)
                                                                                                                   aaTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATg                                <  1:3339853/150‑1 (MQ=255)
                                                                                                                    aTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgc                               <  1:1109550/150‑1 (MQ=255)
                                                                                                                        cTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcagca                           <  1:3355005/150‑1 (MQ=255)
                                                                                                                           gCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCatat                        <  1:2856465/150‑1 (MQ=255)
                                                                                                                            cAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATAtt                       >  2:2213217/1‑150 (MQ=255)
                                                                                                                              gTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCa                     >  2:3643990/1‑150 (MQ=255)
                                                                                                                                ccTGGTCCGGTGCCGACATAGATAAACCGGTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACg                   >  1:4464467/1‑150 (MQ=255)
                                                                                                                                  tGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTa                 <  1:3089294/150‑1 (MQ=255)
                                                                                                                                   ggTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAg                >  1:819891/1‑150 (MQ=255)
                                                                                                                                    gTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAgg               <  2:1271422/150‑1 (MQ=255)
                                                                                                                                    gTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAgg               <  2:2593814/150‑1 (MQ=255)
                                                                                                                                       cGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCaa            <  1:2385109/150‑1 (MQ=255)
                                                                                                                                        ggTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAAt           >  1:3792082/1‑150 (MQ=255)
                                                                                                                                        ggTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAAt           <  1:5144115/150‑1 (MQ=255)
                                                                                                                                          tGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATaa         >  2:4251835/1‑150 (MQ=255)
                                                                                                                                          tGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATaa         >  1:1367036/1‑150 (MQ=255)
                                                                                                                                             ctaCATAGATAAACCTATCACCTTGCATATTCGTTGCCACCAAAATGCTATAGATTATATGGAGTACCATCATCTCATTATGATAGAGTAACGAATTCTGTTTGTACGCTTTGTTCCATTTATGAAGCATATTCACGTACGCAATAAGCt      <  2:3897826/148‑1 (MQ=255)
                                                                                                                                              gACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTg     >  1:1726212/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTT  >  CP009273/3916683‑3916977

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: