Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,669,978 C→A 100% A69A (GCG→GCT)  yphD ← inner membrane putative ABC superfamily sugar transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,669,9780CA97.8% 979.8 / 7.9 312A69A (GCG→GCT) yphDinner membrane putative ABC superfamily sugar transporter permease
Reads supporting (aligned to +/- strand):  ref base C (6/1);  new base A (166/138);  total (172/140)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-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 313 total aligned reads displayed.

CATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGAC  >  CP009273/2669831‑2670127
                                                                                                                                                   |                                                                                                                                                     
cATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGc                                                                                                                                                     >  2:4868748/1‑150 (MQ=255)
      cgcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggc                                                                                                                                               <  2:1486588/150‑1 (MQ=255)
      cgcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggc                                                                                                                                               >  1:3414067/1‑150 (MQ=255)
       gcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggcg                                                                                                                                              >  1:5162869/1‑150 (MQ=255)
        cgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggcgg                                                                                                                                             >  2:1894929/1‑150 (MQ=255)
         gcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggcggc                                                                                                                                            <  1:1574600/150‑1 (MQ=255)
           cTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCaa                                                                                                                                          >  1:886564/1‑150 (MQ=255)
             aaCAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCACTc                                                                                                                                        >  2:4567493/1‑150 (MQ=255)
                 gCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCaa                                                                                                                                    >  2:4628493/1‑150 (MQ=255)
                  cAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAt                                                                                                                                   <  1:1700230/150‑1 (MQ=255)
                     cacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtggt                                                                                                                                >  2:1989467/1‑150 (MQ=255)
                           cagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTggcggc                                                                                                                          >  1:4489469/1‑150 (MQ=255)
                            agGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCa                                                                                                                         >  2:697885/1‑150 (MQ=255)
                             ggCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCAt                                                                                                                        <  2:4147707/150‑1 (MQ=255)
                              gCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATc                                                                                                                       <  2:3125676/150‑1 (MQ=255)
                               cACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACAACAATTTCACCTCAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGTTGGCGGCAAcg                                                                                                                      >  1:1161825/1‑150 (MQ=255)
                                aCGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATcgc                                                                                                                     <  1:4226455/150‑1 (MQ=255)
                                  gCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCgcgc                                                                                                                   >  1:3624386/1‑150 (MQ=255)
                                  cccACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCAgcgcgc                                                                                                                   >  2:173887/2‑150 (MQ=255)
                                   ccACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgca                                                                                                                  <  2:2067775/150‑1 (MQ=255)
                                   ccACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgca                                                                                                                  <  1:3532003/150‑1 (MQ=255)
                                    cACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgcag                                                                                                                 >  1:5115287/1‑150 (MQ=255)
                                      ccGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgcagca                                                                                                               >  1:1911346/1‑150 (MQ=255)
                                       cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                              <  1:5126677/150‑1 (MQ=255)
                                       cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                              <  1:3096274/150‑1 (MQ=255)
                                       cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                              <  1:3222814/150‑1 (MQ=255)
                                          cAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGtt                                                                                                           <  2:2442552/150‑1 (MQ=255)
                                          cAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGtt                                                                                                           <  1:131826/150‑1 (MQ=255)
                                          cAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGtt                                                                                                           <  1:1213968/150‑1 (MQ=255)
                                           aGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTc                                                                                                          <  1:4977427/150‑1 (MQ=255)
                                            gCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCa                                                                                                         <  2:2530917/150‑1 (MQ=255)
                                              ggAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATc                                                                                                       >  2:2132771/1‑150 (MQ=255)
                                              ggAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATc                                                                                                       >  2:4688255/1‑150 (MQ=255)
                                               gAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCt                                                                                                      <  1:4628493/150‑1 (MQ=255)
                                                aaCTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTg                                                                                                     >  2:4470095/1‑150 (MQ=255)
                                                 aCTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTgg                                                                                                    >  1:5058373/1‑150 (MQ=255)
                                                     cAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTAt                                                                                                <  1:428216/150‑1 (MQ=255)
                                                      aaaTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTAtt                                                                                               <  1:3903291/150‑1 (MQ=255)
                                                        aTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCa                                                                                             >  1:1758112/1‑150 (MQ=255)
                                                           gcagcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATg                                                                                          >  2:1852980/1‑150 (MQ=255)
                                                             agcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATgag                                                                                        >  1:2715091/1‑150 (MQ=255)
                                                                aaaaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATa                                                                                     <  1:4858573/150‑1 (MQ=255)
                                                                aaaaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATa                                                                                     >  1:3270212/1‑150 (MQ=255)
                                                                 aaaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaa                                                                                    >  1:849618/1‑150 (MQ=255)
                                                                    aTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGc                                                                                 <  2:3971377/150‑1 (MQ=255)
                                                                      gCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCa                                                                               <  2:2699115/150‑1 (MQ=255)
                                                                       ccAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAg                                                                              <  1:2945171/150‑1 (MQ=255)
                                                                        cAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAgg                                                                             >  2:891999/1‑150 (MQ=255)
                                                                        cAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAgg                                                                             >  1:1442718/1‑150 (MQ=255)
                                                                        cAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAgg                                                                             >  2:4823549/1‑150 (MQ=255)
                                                                         aGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgc                                                                            >  1:3349775/1‑150 (MQ=255)
                                                                          ggCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgcg                                                                           <  2:4259219/150‑1 (MQ=255)
                                                                          ggCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgcg                                                                           >  2:2247101/1‑150 (MQ=255)
                                                                           gCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgcgc                                                                          >  2:3701154/1‑150 (MQ=255)
                                                                            cacaCCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGcgcg                                                                         <  1:478841/150‑1 (MQ=255)
                                                                            cacaCCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGcgcg                                                                         <  1:2447961/150‑1 (MQ=255)
                                                                              cacCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGtt                                                                       <  2:403958/150‑1 (MQ=255)
                                                                                ccGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCa                                                                     <  2:2724816/150‑1 (MQ=255)
                                                                                ccGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCa                                                                     <  2:4269181/150‑1 (MQ=255)
                                                                                  gagaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAgg                                                                   >  1:2001580/1‑150 (MQ=255)
                                                                                         aaGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGagaaga                                                            >  2:1466956/1‑150 (MQ=255)
                                                                                            ccaccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCa                                                         <  2:3880875/150‑1 (MQ=255)
                                                                                             caccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATACAGCCAGGCGCGTTCAGGGAGAAGACCAg                                                        >  1:953540/1‑150 (MQ=255)
                                                                                              accaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGa                                                       <  1:3030080/150‑1 (MQ=255)
                                                                                               ccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGAt                                                      >  1:2980838/1‑150 (MQ=255)
                                                                                                caTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATa                                                     <  2:849340/150‑1 (MQ=255)
                                                                                                caTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATa                                                     >  2:393921/1‑150 (MQ=255)
                                                                                                    ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                                 >  2:5091966/1‑150 (MQ=255)
                                                                                                    ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                                 >  1:3843679/1‑150 (MQ=255)
                                                                                                     gCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGaa                                                >  2:3003901/1‑150 (MQ=255)
                                                                                                      cccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAAt                                               >  2:1679273/1‑150 (MQ=255)
                                                                                                       ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                              <  2:2491098/150‑1 (MQ=255)
                                                                                                          aCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCa                                           >  2:2646814/1‑150 (MQ=255)
                                                                                                            gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                         >  1:2384949/1‑150 (MQ=255)
                                                                                                            gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                         <  2:3333251/150‑1 (MQ=255)
                                                                                                            gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                         >  1:4060203/1‑150 (MQ=255)
                                                                                                              tGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATca                                       >  2:1522576/1‑150 (MQ=255)
                                                                                                                 cATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcacca                                    >  1:2779381/1‑150 (MQ=255)
                                                                                                                  aTCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcaccac                                   >  2:3777977/1‑150 (MQ=255)
                                                                                                                  aTCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcaccac                                   >  2:1333226/1‑150 (MQ=255)
                                                                                                                   tCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCaccacc                                  >  1:801301/1‑150 (MQ=255)
                                                                                                                    cAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCAccacca                                 <  2:2907114/150‑1 (MQ=255)
                                                                                                                    cAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCAccacca                                 >  2:3912370/1‑150 (MQ=255)
                                                                                                                      aTTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGc                               <  2:56824/150‑1 (MQ=255)
                                                                                                                            ccGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCg                         <  1:892035/150‑1 (MQ=255)
                                                                                                                               gagaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtc                      >  1:4041921/1‑150 (MQ=255)
                                                                                                                                 gaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtctc                    <  2:1107767/150‑1 (MQ=255)
                                                                                                                                 gaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtctc                    >  1:1481571/1‑150 (MQ=255)
                                                                                                                                   taataaTCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAt                  >  1:544402/1‑150 (MQ=255)
                                                                                                                                    aataatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAtt                 >  2:1248417/1‑150 (MQ=255)
                                                                                                                                    aataatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAtt                 <  1:1794610/150‑1 (MQ=255)
                                                                                                                                       actCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAat              <  1:3674121/148‑1 (MQ=255)
                                                                                                                                       aatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAat              <  1:1063072/150‑1 (MQ=255)
                                                                                                                                       aatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAat              >  2:4409662/1‑150 (MQ=255)
                                                                                                                                          cAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATg           >  2:1187870/1‑150 (MQ=255)
                                                                                                                                           aGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg          >  1:1001961/1‑150 (MQ=255)
                                                                                                                                             cGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCg        <  1:4299567/150‑1 (MQ=255)
                                                                                                                                              gTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGa       >  1:114164/1‑150 (MQ=255)
                                                                                                                                               tCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGAc      <  2:503600/150‑1 (MQ=255)
                                                                                                                                                cATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACt     >  2:2022878/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGAC  >  CP009273/2669831‑2670127

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