Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,242,574 C→A 9.2% G87* (GGA→TGA)  treA ← periplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,242,5740CA 9.2% 1056.0 / 51.0 281G87* (GGA→TGA) treAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (132/123);  new base A (13/13);  total (145/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 287 total aligned reads displayed.

CATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTG  >  CP009273/1242426‑1242723
                                                                                                                                                    |                                                                                                                                                     
cATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCg                                                                                                                                                      >  1:33739/1‑150 (MQ=255)
cATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCg                                                                                                                                                      >  1:30550/1‑150 (MQ=255)
 aTTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGc                                                                                                                                                     >  2:3678857/1‑150 (MQ=255)
   ttttCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGctc                                                                                                                                                   >  2:2881056/1‑150 (MQ=255)
    tttCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGctct                                                                                                                                                  <  1:3204127/150‑1 (MQ=255)
    tttCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGctct                                                                                                                                                  <  2:3591547/150‑1 (MQ=255)
     ttCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTg                                                                                                                                                 >  1:423079/1‑150 (MQ=255)
       cGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGt                                                                                                                                               <  2:1222901/150‑1 (MQ=255)
       cGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGt                                                                                                                                               <  1:2588428/150‑1 (MQ=255)
          tgtTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGtttt                                                                                                                                            <  2:3303857/150‑1 (MQ=255)
             tttCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTtgct                                                                                                                                         >  1:1999606/1‑150 (MQ=255)
              ttCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTtgctg                                                                                                                                        >  2:2552625/1‑150 (MQ=255)
                cGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCa                                                                                                                                      <  2:4131373/150‑1 (MQ=255)
                cGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCa                                                                                                                                      <  1:3688854/150‑1 (MQ=255)
                cGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCa                                                                                                                                      <  2:1684847/150‑1 (MQ=255)
                 ggTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCAt                                                                                                                                     >  2:2920930/1‑150 (MQ=255)
                 ggTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCAt                                                                                                                                     >  1:962502/1‑150 (MQ=255)
                 ggTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCAt                                                                                                                                     <  1:408448/150‑1 (MQ=255)
                   tAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATcc                                                                                                                                   <  2:3583714/150‑1 (MQ=255)
                     gAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGa                                                                                                                                 <  2:3154333/150‑1 (MQ=255)
                      aaCGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGAt                                                                                                                                >  2:815575/1‑150 (MQ=255)
                      aaCGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGAt                                                                                                                                <  2:2296526/150‑1 (MQ=255)
                       aCGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATa                                                                                                                               <  2:3706683/150‑1 (MQ=255)
                           gTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCa                                                                                                                           >  1:3858384/1‑150 (MQ=255)
                           gTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCa                                                                                                                           >  1:4573977/1‑150 (MQ=255)
                            ttAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAg                                                                                                                          <  2:688170/150‑1 (MQ=255)
                              aaTACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCa                                                                                                                        <  2:3891604/150‑1 (MQ=255)
                                tACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAg                                                                                                                      <  2:3528769/150‑1 (MQ=255)
                                    ggCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatc                                                                                                                  >  2:1783557/1‑150 (MQ=255)
                                      ccAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatcat                                                                                                                >  2:3548939/1‑150 (MQ=255)
                                        aaaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGAtcatca                                                                                                              <  2:4061175/150‑1 (MQ=255)
                                        aaaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGAtcatca                                                                                                              <  2:520948/150‑1 (MQ=255)
                                        aaaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGAtcatca                                                                                                              <  1:2528720/150‑1 (MQ=255)
                                        aaaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGAtcatca                                                                                                              >  1:1788587/1‑150 (MQ=255)
                                         aaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAg                                                                                                             >  1:2840083/1‑150 (MQ=255)
                                         aaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAg                                                                                                             >  2:3984241/1‑150 (MQ=255)
                                          aGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGc                                                                                                            >  2:1480622/1‑150 (MQ=255)
                                           gTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCg                                                                                                           >  1:4277654/1‑150 (MQ=255)
                                           gTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCg                                                                                                           >  1:630428/1‑150 (MQ=255)
                                            tCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCgg                                                                                                          <  1:3306429/150‑1 (MQ=255)
                                               gTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATc                                                                                                       >  1:3227431/1‑150 (MQ=255)
                                                  aaTATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCt                                                                                                    <  1:1179262/150‑1 (MQ=255)
                                                    tatGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGt                                                                                                  <  1:4393987/150‑1 (MQ=255)
                                                     atGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGtt                                                                                                 >  2:3210742/1‑150 (MQ=255)
                                                       gTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcg                                                                                               >  1:334201/1‑150 (MQ=255)
                                                       gTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcg                                                                                               <  2:1824735/150‑1 (MQ=255)
                                                        ttCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcgg                                                                                              <  2:4545292/150‑1 (MQ=255)
                                                          cgcgCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggca                                                                                            >  2:2737841/1‑150 (MQ=255)
                                                          cgcgCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggca                                                                                            <  2:974864/150‑1 (MQ=255)
                                                          cgcgCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggca                                                                                            >  2:794273/1‑150 (MQ=255)
                                                            cgcAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacg                                                                                          <  1:210897/150‑1 (MQ=255)
                                                            cgcAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacg                                                                                          <  1:4474686/150‑1 (MQ=255)
                                                             gcAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacgg                                                                                         <  1:1154552/150‑1 (MQ=255)
                                                               aGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacggca                                                                                       >  2:2960839/1‑150 (MQ=255)
                                                                gTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCAt                                                                                      <  2:1113373/150‑1 (MQ=255)
                                                                 tGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATc                                                                                     <  2:1480371/150‑1 (MQ=255)
                                                                  gACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCg                                                                                    <  2:1581286/150‑1 (MQ=255)
                                                                  gACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCg                                                                                    <  1:14705/150‑1 (MQ=255)
                                                                   aCTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCgg                                                                                   <  2:197539/150‑1 (MQ=255)
                                                                    cTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGc                                                                                  >  1:575646/1‑150 (MQ=255)
                                                                     tGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCa                                                                                 <  2:2109412/150‑1 (MQ=255)
                                                                        cccTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAg                                                                              <  2:3229187/150‑1 (MQ=255)
                                                                         ccTCTGGCGGAACATCTTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAgg                                                                             >  1:3775513/1‑150 (MQ=255)
                                                                         ccTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAgg                                                                             >  2:4441509/1‑150 (MQ=255)
                                                                         ccTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAgg                                                                             >  2:2686367/1‑150 (MQ=255)
                                                                          ctctGGCGGAACATATTTCTCGCCTTCTTTCTGCAGGGTGAAATTGACGTTAACGAATTGTCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATTATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAATGGt                                                                            <  1:1552036/150‑1 (MQ=255)
                                                                                cGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTg                                                                      <  2:3409570/150‑1 (MQ=255)
                                                                                   aaCATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTc                                                                   >  1:3910114/1‑150 (MQ=255)
                                                                                     cATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg                                                                 >  1:904693/1‑150 (MQ=255)
                                                                                          ttCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                            <  1:2335620/150‑1 (MQ=255)
                                                                                           tctcGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaaa                                                           >  1:1352644/1‑150 (MQ=255)
                                                                                           tctcGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaaa                                                           >  1:3168979/1‑150 (MQ=255)
                                                                                              cGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGtt                                                        >  1:1592006/1‑150 (MQ=255)
                                                                                                 cttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTgg                                                     >  1:3793088/1‑150 (MQ=255)
                                                                                                  ttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGc                                                    >  2:3910417/1‑150 (MQ=255)
                                                                                                       tCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGtttt                                               <  1:3951242/150‑1 (MQ=255)
                                                                                                         ggCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGc                                             <  1:37436/150‑1 (MQ=255)
                                                                                                           cAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGcac                                           >  2:4588807/1‑150 (MQ=255)
                                                                                                            aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGAAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGcaca                                          >  1:4157155/1‑150 (MQ=255)
                                                                                                             gggTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCAcat                                         >  1:1407104/1‑150 (MQ=255)
                                                                                                              ggTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCAcatc                                        <  2:2500209/150‑1 (MQ=255)
                                                                                                                 gAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCAtt                                     <  1:1428377/150‑1 (MQ=255)
                                                                                                                  aaaTTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTAGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTa                                    >  2:4007624/1‑150 (MQ=255)
                                                                                                                  aaaTTGACGTTAACGAAATGGCGCAGATCAAATCAGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTa                                    <  1:3864603/150‑1 (MQ=255)
                                                                                                                   aaTTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaa                                   <  1:3511289/150‑1 (MQ=255)
                                                                                                                    aTTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                  >  1:4443632/1‑150 (MQ=255)
                                                                                                                     ttGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAAc                                 <  1:3123828/150‑1 (MQ=255)
                                                                                                                         cGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCg                             >  1:2766415/1‑150 (MQ=255)
                                                                                                                         cGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCg                             <  2:430987/150‑1 (MQ=255)
                                                                                                                            tAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGcc                          >  2:3434426/1‑150 (MQ=255)
                                                                                                                               cGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCta                       <  2:614592/150‑1 (MQ=255)
                                                                                                                                 aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCtaat                     <  2:2925567/150‑1 (MQ=255)
                                                                                                                                 aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCtaat                     >  1:4562891/1‑150 (MQ=255)
                                                                                                                                       cgcAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAata               <  1:1494788/150‑1 (MQ=255)
                                                                                                                                          aGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCa            >  2:262893/1‑150 (MQ=255)
                                                                                                                                             tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAggc         >  1:2695214/1‑150 (MQ=255)
                                                                                                                                              cAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAggcg        <  1:1783557/150‑1 (MQ=255)
                                                                                                                                                 aTCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCt     >  2:2654603/1‑150 (MQ=255)
                                                                                                                                                 aTCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCt     <  2:1161758/150‑1 (MQ=255)
                                                                                                                                                    cGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCtgtg  >  1:1879183/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATTTTTCGGTGTTTTCGGTAGAACGCGTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTG  >  CP009273/1242426‑1242723

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: ‑