Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 54,392 G→T 100% A104E (GCG→GAG)  surA ← peptidyl‑prolyl cis‑trans isomerase (PPIase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927354,3920GT96.8% 784.8 / 14.0 249A104E (GCG→GAG) surApeptidyl‑prolyl cis‑trans isomerase (PPIase)
Reads supporting (aligned to +/- strand):  ref base G (1/7);  new base T (112/128);  total (114/135)
Fisher's exact test for biased strand distribution p-value = 7.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 254 total aligned reads displayed.

AGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCC  >  CP009273/54245‑54538
                                                                                                                                                   |                                                                                                                                                  
aGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTctcg                                                                                                                                                  >  2:1619525/1‑150 (MQ=255)
aGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTctcg                                                                                                                                                  <  1:3531643/150‑1 (MQ=255)
aGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTctcg                                                                                                                                                  <  2:396354/150‑1 (MQ=255)
 ggATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGa                                                                                                                                                 >  2:793466/1‑150 (MQ=255)
      gTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTa                                                                                                                                            <  2:1911464/150‑1 (MQ=255)
      gTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTa                                                                                                                                            <  2:3492310/150‑1 (MQ=255)
          tGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCaa                                                                                                                                        >  2:3029357/1‑150 (MQ=255)
           gCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAAt                                                                                                                                       <  2:2601991/150‑1 (MQ=255)
             gacgacGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCg                                                                                                                                     >  1:2662587/1‑150 (MQ=255)
                gactCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCt                                                                                                                                  <  2:1316002/150‑1 (MQ=255)
                gacgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCt                                                                                                                                  <  1:1290250/150‑1 (MQ=255)
                 acgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTg                                                                                                                                 <  2:3564028/150‑1 (MQ=255)
                   gCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGAt                                                                                                                               >  2:4994185/1‑150 (MQ=255)
                   gCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGAt                                                                                                                               <  2:309401/150‑1 (MQ=255)
                    cACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATc                                                                                                                              >  2:3221822/1‑150 (MQ=255)
                         cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgct                                                                                                                         <  1:711250/150‑1 (MQ=255)
                            tgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCCTGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgctgct                                                                                                                      >  2:2436886/1‑150 (MQ=255)
                            tgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgctgct                                                                                                                      <  1:20395/150‑1 (MQ=255)
                             gttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTc                                                                                                                     <  1:2599575/150‑1 (MQ=255)
                             gttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTc                                                                                                                     <  2:1572130/150‑1 (MQ=255)
                                aCGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATc                                                                                                                  <  2:2675355/150‑1 (MQ=255)
                                   cACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGga                                                                                                               >  2:1052092/1‑150 (MQ=255)
                                   cACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGga                                                                                                               >  2:2773426/1‑150 (MQ=255)
                                   cACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGga                                                                                                               <  1:1121774/150‑1 (MQ=255)
                                    aCTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgag                                                                                                              <  2:1775138/150‑1 (MQ=255)
                                    aCTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgag                                                                                                              <  2:1176594/150‑1 (MQ=255)
                                     cTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgaga                                                                                                             <  1:690031/150‑1 (MQ=255)
                                     cTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgaga                                                                                                             <  2:2887241/150‑1 (MQ=255)
                                         agagATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGAtttt                                                                                                         <  1:529493/150‑1 (MQ=255)
                                           agaTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCa                                                                                                       <  1:2468317/150‑1 (MQ=255)
                                             aTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACt                                                                                                     >  2:5034849/1‑150 (MQ=255)
                                               aaTCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTcc                                                                                                   >  1:311177/1‑150 (MQ=255)
                                                  catcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAt                                                                                                <  1:3865634/150‑1 (MQ=255)
                                                    tctcTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAttt                                                                                              <  1:147557/150‑1 (MQ=255)
                                                     ctctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAtttt                                                                                             <  2:1713968/150‑1 (MQ=255)
                                                     ctctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAtttt                                                                                             <  1:4625207/150‑1 (MQ=255)
                                                     ctctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAtttt                                                                                             <  1:2342991/150‑1 (MQ=255)
                                                      tctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTc                                                                                            <  1:4197941/150‑1 (MQ=255)
                                                      tctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTc                                                                                            <  1:157442/150‑1 (MQ=255)
                                                      tctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTc                                                                                            >  1:5130048/1‑150 (MQ=255)
                                                       ctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCt                                                                                           <  1:1234201/150‑1 (MQ=255)
                                                       ctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCt                                                                                           >  1:4514255/1‑150 (MQ=255)
                                                        tttGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTg                                                                                          <  1:4477316/150‑1 (MQ=255)
                                                           gCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGccc                                                                                       >  1:5079348/1‑150 (MQ=255)
                                                            cGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGcccc                                                                                      >  1:2547388/1‑150 (MQ=255)
                                                             ggATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCa                                                                                     >  2:3429550/1‑150 (MQ=255)
                                                               aTCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATc                                                                                   <  2:155781/150‑1 (MQ=255)
                                                                tCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCt                                                                                  <  2:2540465/150‑1 (MQ=255)
                                                                   ggTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCa                                                                               >  2:2068450/1‑150 (MQ=255)
                                                                    gTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAg                                                                              >  1:1655034/1‑150 (MQ=255)
                                                                          aTAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatgat                                                                        >  1:733394/1‑150 (MQ=255)
                                                                          aTAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatgat                                                                        <  2:3675029/150‑1 (MQ=255)
                                                                            aGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGAttt                                                                      >  2:2314346/1‑150 (MQ=255)
                                                                                gttgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATc                                                                  >  1:1888285/1‑150 (MQ=255)
                                                                                  tgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCa                                                                >  2:1432287/1‑150 (MQ=255)
                                                                                     aGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGa                                                             <  1:4736449/150‑1 (MQ=255)
                                                                                      gTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAt                                                            >  1:4363275/1‑150 (MQ=255)
                                                                                           gTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAAc                                                       <  1:2513955/150‑1 (MQ=255)
                                                                                             ccATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGt                                                     >  1:2772726/1‑150 (MQ=255)
                                                                                              cATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGtt                                                    <  2:4463600/150‑1 (MQ=255)
                                                                                              cATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGGTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGtt                                                    <  2:4722760/150‑1 (MQ=255)
                                                                                                  gTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCAt                                                <  2:3716230/150‑1 (MQ=255)
                                                                                                   tAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATg                                               <  1:4091041/150‑1 (MQ=255)
                                                                                                   tAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATg                                               >  2:4900688/1‑150 (MQ=255)
                                                                                                    aaGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGa                                              >  1:3186952/1‑150 (MQ=255)
                                                                                                              gCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGgcg                                    <  2:2254729/150‑1 (MQ=255)
                                                                                                               cTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGgcgc                                   >  1:109122/1‑150 (MQ=255)
                                                                                                                tGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCa                                  >  2:1819829/1‑150 (MQ=255)
                                                                                                                tGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCa                                  <  1:892477/150‑1 (MQ=255)
                                                                                                                tGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCa                                  <  1:3380732/150‑1 (MQ=255)
                                                                                                                 gcgcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAg                                 >  2:4101492/1‑150 (MQ=255)
                                                                                                                  cgcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgc                                >  2:735368/1‑150 (MQ=255)
                                                                                                                   gcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcg                               <  1:3029357/150‑1 (MQ=255)
                                                                                                                   gcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcg                               >  2:2269915/1‑150 (MQ=255)
                                                                                                                   gcATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcg                               <  1:3725924/150‑1 (MQ=255)
                                                                                                                    cATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgt                              >  1:938432/1‑150 (MQ=255)
                                                                                                                      tCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcg                            >  1:1150408/1‑150 (MQ=255)
                                                                                                                      tCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcg                            <  2:4412407/150‑1 (MQ=255)
                                                                                                                       cTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgc                           >  1:1921100/1‑150 (MQ=255)
                                                                                                                       cTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgc                           >  1:4432479/1‑150 (MQ=255)
                                                                                                                        tGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcg                          <  1:2910805/150‑1 (MQ=255)
                                                                                                                         gATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcgt                         >  2:3596885/1‑150 (MQ=255)
                                                                                                                           tCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtca                       >  1:3322710/1‑150 (MQ=255)
                                                                                                                           tCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtca                       <  1:2797956/150‑1 (MQ=255)
                                                                                                                           tCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtca                       <  2:455293/150‑1 (MQ=255)
                                                                                                                            ccAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcat                      <  1:4222009/150‑1 (MQ=255)
                                                                                                                              aGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcatca                    >  1:1344502/1‑150 (MQ=255)
                                                                                                                               gCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAg                   <  2:626172/150‑1 (MQ=255)
                                                                                                                                 gTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGa                 <  1:2212903/150‑1 (MQ=255)
                                                                                                                                    aTGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGt              <  1:4103865/150‑1 (MQ=255)
                                                                                                                                     tgttgtTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGtt             <  2:174495/150‑1 (MQ=255)
                                                                                                                                     tgttgtTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGtt             >  2:4690450/1‑150 (MQ=255)
                                                                                                                                     tgttgtTCTGTGTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGtt             <  1:4955541/150‑1 (MQ=255)
                                                                                                                                      gttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtg            >  1:4462795/1‑150 (MQ=255)
                                                                                                                                       ttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgc           >  2:4640430/1‑150 (MQ=255)
                                                                                                                                        tgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgct          <  1:1052092/150‑1 (MQ=255)
                                                                                                                                         gttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctg         <  2:2987169/150‑1 (MQ=255)
                                                                                                                                          ttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctgc        <  1:3919711/150‑1 (MQ=255)
                                                                                                                                            cTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcct      >  2:2563889/1‑150 (MQ=255)
                                                                                                                                              gTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgccttg    >  2:2721264/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGGATGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCC  >  CP009273/54245‑54538

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: