Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,952,842 A→C 10.9% intergenic (+37/+50) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,952,8420AC10.9% 309.8 / 4.5 135intergenic (+37/+50)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base A (36/79);  new base C (17/3);  total (53/82)
Fisher's exact test for biased strand distribution p-value = 8.16e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

ACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAC  >  CP009273/3952693‑3952991
                                                                                                                                                     |                                                                                                                                                     
aCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                       <  2:3782666/150‑1 (MQ=255)
 cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                                      <  2:5158276/150‑1 (MQ=255)
   ggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATACGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACaggg                                                                                                                                                    >  2:5108191/1‑150 (MQ=255)
     cAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggag                                                                                                                                                  <  1:1965132/150‑1 (MQ=255)
      aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCCCCCCGACCCACACCCCCcgggggg                                                                                                                                                 >  1:1244153/1‑147 (MQ=255)
        cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGTCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                               <  1:389999/150‑1 (MQ=255)
        cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                               <  1:4542065/150‑1 (MQ=255)
        cTGCGTGAGGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                               <  1:4542087/150‑1 (MQ=255)
         tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGAGGCCCCCACCCCCACCCACAACCACacggggcggg                                                                                                                                              >  1:2040461/1‑150 (MQ=255)
            gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAgggagagggaga                                                                                                                                           <  1:4216667/150‑1 (MQ=255)
            gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAgggagagggaga                                                                                                                                           <  1:3738810/150‑1 (MQ=255)
             tGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATCTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACTCCCACAGGGGGAGGGAGaa                                                                                                                                          >  1:2195544/1‑150 (MQ=255)
                tgtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGTGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAac                                                                                                                                       <  1:4263155/150‑1 (MQ=255)
                    aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCTCCCACAGGGAGAGGGCGAAAACACTc                                                                                                                                   >  2:4703362/1‑150 (MQ=255)
                      cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGCGAGGGAGAAAACACTaaa                                                                                                                                 >  2:4469371/1‑150 (MQ=255)
                       gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                <  2:4117500/150‑1 (MQ=255)
                       gAAGCGAGTCGCAGCCCTGCGATTGAGCAGGTAGGTAAGACACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                <  2:4117425/150‑1 (MQ=255)
                         aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                              >  2:5099904/1‑150 (MQ=255)
                           cGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCGCTCCCACAGGGGGAGGGAGAAAACACACAAGGcca                                                                                                                            >  2:2437842/1‑149 (MQ=255)
                             aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGCGAGGGAGAAAAAACACAAGGCCAtc                                                                                                                          >  2:758595/1‑150 (MQ=39)
                              ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                                                                         <  1:3283156/150‑1 (MQ=255)
                              ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                                                                         <  2:3981585/150‑1 (MQ=255)
                               tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                                                        <  2:3537406/150‑1 (MQ=255)
                                 gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                      <  1:3618943/150‑1 (MQ=255)
                                 gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                      <  2:4396126/150‑1 (MQ=255)
                                 gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                      <  1:1649967/150‑1 (MQ=255)
                                    gcgaTGCGATTGAGCAGGTAGGTAAGAAACTACGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                   <  2:755939/147‑1 (MQ=255)
                                     ccAAGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgag                                                                                                                  <  1:2853314/150‑1 (MQ=255)
                                         gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCCCCCACACCCACCCGGAGCGGGAGAAAACACACAAAGCCATCTCCTGGCGACgg                                                                                                              >  2:1869605/1‑150 (MQ=255)
                                          cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCCCCCCGACCCTCAGCCACAGGGAGTGGGGGAAAAAACTCAAGGCCGTCTCCGGGGGAAggg                                                                                                             >  2:3678501/1‑149 (MQ=255)
                                            aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCAACCCGACCCACTCCCACAGGGATAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGcca                                                                                                           >  2:1718027/1‑149 (MQ=255)
                                            aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACTCCCACAGGGAGCGGGAGAAAACACACAACGCCGTCTCCGGGAGAAGGccg                                                                                                           >  2:4587721/1‑149 (MQ=255)
                                              tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTACCCCAGGGGGAGGGAGTAAACACTCAAGGCCCTCTCCTGGTGAAGGCCCTg                                                                                                         >  1:4730193/1‑150 (MQ=255)
                                                aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                       <  1:1148094/150‑1 (MQ=255)
                                                      taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                                                                                 <  1:2614557/150‑1 (MQ=255)
                                                        ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACATGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                               <  2:1947621/150‑1 (MQ=255)
                                                        ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCAACCACACCCCCCGGGTGGGGGAGAAAAAACACAAGGCCGTATACTGGAGAAGGCCATGCAATTAtttt                                                                                               >  2:2059118/1‑149 (MQ=255)
                                                         gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                              <  2:2223706/150‑1 (MQ=255)
                                                          taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCTCACCCCCAGGGCGAGGGCGAAAACACGCAAGGCCATGTGCAGGCGAAGGCCAGGCTAATAGTTTCg                                                                                             >  2:712333/1‑150 (MQ=255)
                                                              aaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCAACCCCACCCACGACCACAGGGAGAGGGGGTAAAAAAACAAGGCCGTCTCCGGGGGAAGGCCGGGCGATTAGTTTGGGTaa                                                                                         >  1:3875807/1‑149 (MQ=255)
                                                               aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                                                        <  1:4960844/150‑1 (MQ=255)
                                                                    cgcgGCTATATGACAGAGATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                   <  2:1952845/150‑1 (MQ=255)
                                                                      cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                 <  2:1276133/150‑1 (MQ=255)
                                                                      cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                 <  2:3332243/150‑1 (MQ=255)
                                                                      cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                 <  1:1642312/150‑1 (MQ=255)
                                                                      cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                 <  1:1681727/150‑1 (MQ=255)
                                                                      cgGCGATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCGCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCGGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                 <  2:3332263/150‑1 (MQ=255)
                                                                        gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                               <  1:4898774/150‑1 (MQ=255)
                                                                            tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACCGGGAGAGGGAGAAAACAAACAAGGCCATCACCGGGCGAAGGCCCTGCTATTAGTTTCGGTACAGCCCCTTAatgat                                                                           >  2:4590068/1‑150 (MQ=37)
                                                                            tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACACCCACCGGGACAGGGCGAAAACACACAACGCCGTCTCCGGGGGAAGGCCGTGCAATTAGTTGCGGTACAGCACCATAatgtt                                                                           >  2:4589706/1‑148 (MQ=255)
                                                                             atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                          <  2:4819196/150‑1 (MQ=255)
                                                                             atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAATGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                          >  2:3223624/1‑150 (MQ=255)
                                                                              tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCACACCCCCACCCACAACCACACGGAGACGGAGAAAACACACAAGGCCATAGCCGGGAGACCGCCATGCAAATAGTTTCGGTTCAGCACCATAatggttt                                                                         >  2:2818542/1‑146 (MQ=255)
                                                                               gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGTGAAAACACTCAAGGCCGTCAACGGGAGAAGGCCGTGCTATTAGTTGCGGGGCCGCACCATAATGAggtg                                                                        >  2:4940887/1‑150 (MQ=255)
                                                                                 cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                      <  2:4024014/150‑1 (MQ=255)
                                                                                 cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                      <  1:3187155/150‑1 (MQ=255)
                                                                                    atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTAACCCCGACCCTAGCCCACAGGGGGGGGGTGAAAAAACTCAAGGCCGTAACCGGGCGCAGGGCGTGCTATTAGGTTCGGTAACGCCCCATAATGTTGTGGTATc                                                                   >  2:3650504/1‑150 (MQ=255)
                                                                                      atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCAAAGGGAGAGGGAGAAAAAACTCAAGGCCTACTCCTGGAGAAGGCCTTTCTATTATATGCGGTACAGCACCAAAATTATGTGatatacc                                                                 >  2:554653/1‑147 (MQ=37)
                                                                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACACAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                              >  2:3854221/1‑150 (MQ=255)
                                                                                           aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTg                                                            <  1:3864667/150‑1 (MQ=255)
                                                                                             gtcttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACACCCACCGGGAGGGGGAGAAAAAACTAAATGCCGTATCCTGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCAACTTAATGATGTGTTATCCGAACTggg                                                          >  2:2581515/4‑150 (MQ=37)
                                                                                                 tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCACCCCCACCCACACCCCCACGGAGAGGGCTAAAACCCACAAAGCCATCTACAGGCGAAGGCCAGGCGATTAGTTGCGGTACAGCACCATAATGATGTGAATTCCCAACTGGGggtg                                                      >  2:3477162/1‑150 (MQ=255)
                                                                                                     gTTGCGGGTTAAGTGCGCGCTGATTCCCACACCCCGCACCTCACCCACAGGGCGAGGGAGAAAAAACACAAGGCCTTCACCAGGACAAAGCCATGCTAATAGTTTCGGTACAGCACCTTAAAGATGTGATATACGCACAGGGTGTGCAGc                                                  >  2:4044060/1‑150 (MQ=255)
                                                                                                     gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGc                                                  <  2:83336/150‑1 (MQ=255)
                                                                                                         cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACCGGGAGAGGGAGAAAACACGCAACGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTTCGGTACAGCACCTTAATGATTTGATATCCGAACAGGGGGTGCAGCGGGc                                              >  2:5080402/1‑150 (MQ=37)
                                                                                                          gggTTAAGTGCGCGCTGATGCCCACACCCCGACCCACACCCACCGGGAGAGGGAGAAAACAATCAATGCCTTCACCATGAGAAGGCCATGCTAATAGTTGCGGTACAGCACCTTAATGATGTGATATACCAACTGGGTGTGCAGCGGGcc                                             >  2:4728442/1‑150 (MQ=255)
                                                                                                                aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAACCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                       <  1:1771994/150‑1 (MQ=255)
                                                                                                                    cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                   <  2:3098496/150‑1 (MQ=255)
                                                                                                                    cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                   <  2:3098444/150‑1 (MQ=255)
                                                                                                                      cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAg                                 <  1:2452290/150‑1 (MQ=255)
                                                                                                                       gcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGt                                <  2:3245002/150‑1 (MQ=255)
                                                                                                                       gcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGt                                <  2:3244894/150‑1 (MQ=255)
                                                                                                                       gcTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACGGGGTGTGCAGCGGGCCGGTCGGCTCCAGt                                >  2:2169034/1‑150 (MQ=255)
                                                                                                                        cTGAGGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                               <  2:241456/150‑1 (MQ=255)
                                                                                                                         tGATGCCCTCACCCCGACCCTCTCCCACAGGGCGAGGGAGAAAACACTCAAGGCCTTCTCCTGGCGAAGGCCTTGCGTTTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                              >  2:4861894/1‑150 (MQ=255)
                                                                                                                         tGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                              >  2:618182/1‑150 (MQ=255)
                                                                                                                          gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                             >  2:1238574/1‑150 (MQ=255)
                                                                                                                           aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                            <  2:1015516/150‑1 (MQ=255)
                                                                                                                           aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                            <  2:1015883/150‑1 (MQ=255)
                                                                                                                             gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                          <  1:2129593/150‑1 (MQ=255)
                                                                                                                              cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCGATTCGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAc                         <  2:1621159/150‑1 (MQ=255)
                                                                                                                               ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                        <  1:2647574/150‑1 (MQ=255)
                                                                                                                                cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                       <  2:943505/150‑1 (MQ=255)
                                                                                                                                cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                       <  2:2509848/150‑1 (MQ=255)
                                                                                                                                  cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                     <  2:1543739/150‑1 (MQ=255)
                                                                                                                                    ccccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTTCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                   >  1:2150762/1‑150 (MQ=255)
                                                                                                                                     cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                  <  1:2960374/150‑1 (MQ=255)
                                                                                                                                      ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCCAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaa                 <  1:4353651/150‑1 (MQ=255)
                                                                                                                                        gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc               <  1:1930386/150‑1 (MQ=255)
                                                                                                                                        gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc               <  1:5134378/150‑1 (MQ=255)
                                                                                                                                         aCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAcc              <  1:3900133/150‑1 (MQ=255)
                                                                                                                                          cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa             <  1:294454/150‑1 (MQ=255)
                                                                                                                                           ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAAAACCAc            >  1:316651/1‑150 (MQ=255)
                                                                                                                                             tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt          <  1:4958918/150‑1 (MQ=255)
                                                                                                                                              ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACttt         >  1:828375/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc      <  1:3375687/150‑1 (MQ=255)
                                                                                                                                                 ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc      <  1:770950/150‑1 (MQ=255)
                                                                                                                                                   acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGa    <  1:268185/150‑1 (MQ=255)
                                                                                                                                                     agggagagggagAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAc  <  2:2345344/150‑1 (MQ=255)
                                                                                                                                                     agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGACCTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGACAACCACTTGATCGAAc  >  2:682361/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAAC  >  CP009273/3952693‑3952991

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑