Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,295,611 C→T 7.8% A208T (GCG→ACG)  napA ← nitrate reductase, periplasmic, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,295,6110CT 7.8% 1095.4 / 56.6 320A208T (GCG→ACG) napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base C (166/129);  new base T (8/17);  total (174/146)
Fisher's exact test for biased strand distribution p-value = 2.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 326 total aligned reads displayed.

GATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGC  >  CP009273/2295462‑2295758
                                                                                                                                                     |                                                                                                                                                   
   gCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAGACGCGTc                                                                                                                                                  <  2:842068/150‑1 (MQ=255)
   gCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTc                                                                                                                                                  <  1:3863403/150‑1 (MQ=255)
   gCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTc                                                                                                                                                  <  2:4491370/150‑1 (MQ=255)
      gTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGc                                                                                                                                               >  1:1912865/1‑150 (MQ=255)
      gTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGc                                                                                                                                               >  2:3904332/1‑150 (MQ=255)
       ttATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGTGTCAGcc                                                                                                                                              <  1:111145/150‑1 (MQ=255)
       ttATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGTGTCAGcc                                                                                                                                              <  1:4313426/150‑1 (MQ=255)
        tATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCt                                                                                                                                             <  1:4304306/150‑1 (MQ=255)
        tATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCt                                                                                                                                             >  2:1126825/1‑150 (MQ=255)
           ccGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCt                                                                                                                                          <  1:2681520/150‑1 (MQ=255)
            cGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTc                                                                                                                                         <  1:2410022/150‑1 (MQ=255)
            cGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTc                                                                                                                                         >  1:3599184/1‑150 (MQ=255)
               cAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGAt                                                                                                                                      <  2:1000595/150‑1 (MQ=255)
                aGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGTGTCAGCCTGCTCGATg                                                                                                                                     <  1:3443527/150‑1 (MQ=255)
                    cGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATgtcgt                                                                                                                                 <  1:5134347/150‑1 (MQ=255)
                     gAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATgtcgtc                                                                                                                                >  2:645037/1‑150 (MQ=255)
                       aGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCAt                                                                                                                              >  1:1695277/1‑150 (MQ=255)
                       aGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCAt                                                                                                                              >  1:1208071/1‑150 (MQ=255)
                       aGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCAt                                                                                                                              >  1:4049370/1‑150 (MQ=255)
                         cTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATag                                                                                                                            >  2:1687714/1‑150 (MQ=255)
                            cGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATagcag                                                                                                                         <  2:207443/150‑1 (MQ=255)
                                 cTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCAt                                                                                                                    >  2:23085/1‑150 (MQ=255)
                                  tGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATc                                                                                                                   <  1:527883/150‑1 (MQ=255)
                                       ggtagAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTc                                                                                                              <  1:3662962/150‑1 (MQ=255)
                                           gAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATcc                                                                                                          >  2:1065138/1‑150 (MQ=255)
                                           gAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATcc                                                                                                          <  1:1498909/150‑1 (MQ=255)
                                            aaaGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCa                                                                                                         >  1:652132/1‑150 (MQ=255)
                                              aGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATa                                                                                                       <  2:1027298/150‑1 (MQ=255)
                                               gCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATAc                                                                                                      >  1:3201518/1‑150 (MQ=255)
                                                     ccACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAgg                                                                                                <  1:2258737/150‑1 (MQ=255)
                                                      cACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGt                                                                                               <  2:2468945/150‑1 (MQ=255)
                                                      cACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGt                                                                                               >  2:94133/1‑150 (MQ=255)
                                                      cACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGt                                                                                               >  1:1309028/1‑150 (MQ=255)
                                                         ggTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAGACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACg                                                                                            <  1:3078861/150‑1 (MQ=255)
                                                         ggTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACg                                                                                            >  1:462690/1‑150 (MQ=255)
                                                          gTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGc                                                                                           >  1:3272658/1‑150 (MQ=255)
                                                           tGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGGCAGCCTGCGCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCa                                                                                          >  2:143785/1‑150 (MQ=255)
                                                               gTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATaaa                                                                                      >  1:5673197/1‑150 (MQ=255)
                                                                ttCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAg                                                                                     <  2:885222/150‑1 (MQ=255)
                                                                  cTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGcc                                                                                   <  1:3904332/150‑1 (MQ=255)
                                                                    ggTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCaa                                                                                 >  1:4494767/1‑150 (MQ=255)
                                                                    ggTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCaa                                                                                 <  2:2916347/150‑1 (MQ=255)
                                                                      ttAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAact                                                                               <  1:2519043/150‑1 (MQ=255)
                                                                      ttAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAact                                                                               <  1:4939457/150‑1 (MQ=255)
                                                                       tAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAGACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAacta                                                                              >  1:3727695/1‑150 (MQ=255)
                                                                       tAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCCTAATGCCAacta                                                                              >  2:3035501/1‑150 (MQ=255)
                                                                       tAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAacta                                                                              <  1:3986133/150‑1 (MQ=255)
                                                                        agagagACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAactac                                                                             <  2:683478/150‑1 (MQ=255)
                                                                         gagagaCGACGGTTAGTGATGCGCGACCAGATGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACTAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAactact                                                                            <  2:4795497/150‑1 (MQ=255)
                                                                            agaCGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGcc                                                                         <  1:4152709/150‑1 (MQ=255)
                                                                             gacgacGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAGACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCAGACCAAAGGTACGCATAAAGCCAACTACTGCCg                                                                        <  2:3696212/150‑1 (MQ=255)
                                                                             gacgacGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGTGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCg                                                                        <  1:5552917/150‑1 (MQ=255)
                                                                              acgacgGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGa                                                                       >  1:5095546/1‑150 (MQ=255)
                                                                                gacgGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACg                                                                     >  1:5059402/1‑150 (MQ=255)
                                                                                  cgGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGcc                                                                   <  1:3224280/150‑1 (MQ=255)
                                                                                  cgGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGcc                                                                   <  2:2852766/150‑1 (MQ=255)
                                                                                      tAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATAc                                                               >  2:3075003/1‑150 (MQ=255)
                                                                                       aGTGATGCGCGACCAGAGGATCGGGAGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCATCCTGCTCGATTTCGACATAGCAGCCCATCGGCTCATCAATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACa                                                              >  2:5189970/1‑150 (MQ=255)
                                                                                        gTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAg                                                             >  1:86836/1‑150 (MQ=255)
                                                                                         tGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGt                                                            >  2:2622370/1‑150 (MQ=255)
                                                                                         tGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGt                                                            <  2:1551568/150‑1 (MQ=255)
                                                                                           aTGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGTGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGa                                                          >  2:4339793/1‑150 (MQ=255)
                                                                                           aTGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGa                                                          <  2:2305077/150‑1 (MQ=255)
                                                                                           aTGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGa                                                          <  1:294253/150‑1 (MQ=255)
                                                                                               gcgACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGAcgcg                                                      <  1:5652446/150‑1 (MQ=255)
                                                                                                      gagGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCgg                                               <  2:124831/150‑1 (MQ=255)
                                                                                                      gagGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCgg                                               >  2:2766131/1‑150 (MQ=255)
                                                                                                        ggATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGt                                             <  2:5052067/150‑1 (MQ=255)
                                                                                                        ggATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGt                                             >  1:2753578/1‑150 (MQ=255)
                                                                                                         gATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTc                                            >  2:4094136/1‑150 (MQ=255)
                                                                                                          aTCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGTGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCg                                           <  2:3422066/150‑1 (MQ=255)
                                                                                                             gggTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGAtg                                        <  1:4743176/150‑1 (MQ=255)
                                                                                                                tGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGAtgttg                                     <  2:4665925/150‑1 (MQ=255)
                                                                                                                 gCATTTCCGCCATGTTTGCGCCCCACAGCACAGACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGAtgttgt                                    >  1:2992316/1‑150 (MQ=255)
                                                                                                                 gCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGAtgttgt                                    >  2:5354657/1‑150 (MQ=255)
                                                                                                                 gCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGAtgttgt                                    >  2:3367266/1‑150 (MQ=255)
                                                                                                                 gCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGAtgttgt                                    >  1:485532/1‑150 (MQ=255)
                                                                                                                   aTTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTc                                  <  2:4249894/150‑1 (MQ=255)
                                                                                                                      tCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGaa                               >  2:692266/1‑150 (MQ=255)
                                                                                                                      tCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGaa                               >  2:1854463/1‑150 (MQ=255)
                                                                                                                       ccGCCATGTTTGCGCCCCACCGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAAc                              <  2:4957548/150‑1 (MQ=255)
                                                                                                                       ccGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAAc                              <  1:768265/150‑1 (MQ=255)
                                                                                                                       ccGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAAc                              <  2:365595/150‑1 (MQ=255)
                                                                                                                        cGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACg                             <  2:242551/150‑1 (MQ=255)
                                                                                                                        cGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACg                             >  1:2347816/1‑150 (MQ=255)
                                                                                                                          ccATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGa                           <  2:4864671/150‑1 (MQ=255)
                                                                                                                             tGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGc                        <  1:3158313/150‑1 (MQ=255)
                                                                                                                              gTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGcc                       >  2:1773800/1‑150 (MQ=255)
                                                                                                                                 tGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGc                    >  1:5373155/1‑150 (MQ=255)
                                                                                                                                   cgcCCCACAGCACAGACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCCTAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCtt                  >  2:5658970/1‑150 (MQ=255)
                                                                                                                                    gcCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCttt                 <  2:2130227/150‑1 (MQ=255)
                                                                                                                                        cacaGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAAc             <  2:3192268/150‑1 (MQ=255)
                                                                                                                                           aGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGc          >  2:2455228/1‑150 (MQ=255)
                                                                                                                                             cacaAACGTGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCtt        >  1:747768/1‑150 (MQ=255)
                                                                                                                                              acaAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTg       >  2:192989/1‑150 (MQ=255)
                                                                                                                                               caAACGTGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTgg      <  1:2156281/150‑1 (MQ=255)
                                                                                                                                                aaaCGTGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGa     <  2:59796/150‑1 (MQ=255)
                                                                                                                                                aaaCGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGa     <  1:1065138/150‑1 (MQ=255)
                                                                                                                                                 aaCGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGAc    >  1:2051589/1‑150 (MQ=255)
                                                                                                                                                  aCGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGAcg   >  1:2664811/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGC  >  CP009273/2295462‑2295758

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: