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,9780CA98.1% 867.7 / 3.9 278A69A (GCG→GCTyphDinner membrane putative ABC superfamily sugar transporter permease
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base A (147/125);  total (148/130)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 278 total aligned reads displayed.

CATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACT  >  CP009273/2669831‑2670124
                                                                                                                                                   |                                                                                                                                                  
cATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGc                                                                                                                                                  <  1:4402804/150‑1 (MQ=255)
 aTCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGcc                                                                                                                                                 >  2:3037068/1‑150 (MQ=255)
 aTCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGcc                                                                                                                                                 >  1:1156459/1‑150 (MQ=255)
 aTCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGcc                                                                                                                                                 <  1:3842022/150‑1 (MQ=255)
       gcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggcg                                                                                                                                           <  1:4957720/150‑1 (MQ=255)
          ccTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCa                                                                                                                                        >  1:1005768/1‑150 (MQ=255)
          ccTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCa                                                                                                                                        >  1:1634454/1‑150 (MQ=255)
           cTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCac                                                                                                                                       >  2:5322726/1‑149 (MQ=255)
                aGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCa                                                                                                                                  >  1:1085457/1‑150 (MQ=255)
                 gCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCaa                                                                                                                                 <  1:870237/150‑1 (MQ=255)
                  cAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAt                                                                                                                                <  1:4131423/150‑1 (MQ=255)
                   aaCACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAtg                                                                                                                               >  1:2483379/1‑150 (MQ=255)
                    acacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAtgg                                                                                                                              <  1:5221496/150‑1 (MQ=255)
                     cacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAtggt                                                                                                                             <  2:5109382/150‑1 (MQ=255)
                      acCAGCAGGCACGCCACCGCCAGCGGAACTTCACATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAtggtg                                                                                                                            <  2:1207160/150‑1 (MQ=255)
                        cagcagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTggc                                                                                                                          <  2:2934836/150‑1 (MQ=255)
                            agGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCa                                                                                                                      >  2:3363912/1‑150 (MQ=255)
                                aCGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATcgc                                                                                                                  <  2:1052318/150‑1 (MQ=255)
                                  gCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCgcgc                                                                                                                >  2:1587724/1‑150 (MQ=255)
                                       cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                           <  1:3808852/150‑1 (MQ=255)
                                        gCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACg                                                                                                          <  2:4900557/150‑1 (MQ=255)
                                        gCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACg                                                                                                          <  1:2870978/150‑1 (MQ=255)
                                         ccAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGt                                                                                                         >  2:4190367/1‑150 (MQ=255)
                                         ccAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGt                                                                                                         >  2:2493811/1‑150 (MQ=255)
                                            gCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCa                                                                                                      <  1:4781389/150‑1 (MQ=255)
                                                 aCTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTgg                                                                                                 <  2:2326033/150‑1 (MQ=255)
                                                   ttCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGtt                                                                                               <  1:3383478/150‑1 (MQ=255)
                                                   ttCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGtt                                                                                               >  1:3465275/1‑150 (MQ=255)
                                                    tCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTa                                                                                              >  1:2364672/1‑150 (MQ=255)
                                                      aaaTTGAAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTAtt                                                                                            >  1:4962499/1‑150 (MQ=255)
                                                       aaTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTc                                                                                           <  1:4847439/150‑1 (MQ=255)
                                                          tGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAAt                                                                                        <  1:3882665/150‑1 (MQ=255)
                                                             agcaAACATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCACCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATgag                                                                                     <  2:4452477/150‑1 (MQ=255)
                                                               caAAAATGCCAGGCACAACGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGAt                                                                                   <  2:3783621/150‑1 (MQ=255)
                                                                  aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                >  1:3603699/1‑150 (MQ=255)
                                                                  aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                <  1:2640286/150‑1 (MQ=255)
                                                                  aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                <  2:774420/150‑1 (MQ=255)
                                                                    aTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGc                                                                              <  2:390275/150‑1 (MQ=255)
                                                                     tGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGcc                                                                             >  1:3069378/1‑150 (MQ=255)
                                                                       ccAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAg                                                                           <  1:1852825/150‑1 (MQ=255)
                                                                        cAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAgg                                                                          >  1:4703963/1‑150 (MQ=255)
                                                                         aGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgc                                                                         <  1:1497942/150‑1 (MQ=255)
                                                                          ggCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgcg                                                                        <  2:203341/150‑1 (MQ=255)
                                                                           gCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgcgc                                                                       >  2:944074/1‑150 (MQ=255)
                                                                           gCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgcgc                                                                       >  2:1879186/1‑150 (MQ=255)
                                                                              cacCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGtt                                                                    <  1:485438/150‑1 (MQ=255)
                                                                              cacCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGtt                                                                    >  1:3688110/1‑150 (MQ=255)
                                                                              cacCGAGACAAAAGCCACCATCGGCCCAACGCTGACAACAATTTCACCGGAGATAATAAACAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCTCGCAGCACGTGAATCTGGTTATTCAATGAGATAAAGCCAGGCGCGtt                                                                    >  1:3886925/1‑150 (MQ=255)
                                                                               acCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTc                                                                   >  2:2241823/1‑150 (MQ=255)
                                                                                 cGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAg                                                                 <  2:4741860/150‑1 (MQ=255)
                                                                                  gagaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAgg                                                                <  1:1924925/150‑1 (MQ=255)
                                                                                   agaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGTATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAggg                                                               >  2:2162838/1‑150 (MQ=255)
                                                                                   agaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAggg                                                               <  1:2586176/150‑1 (MQ=255)
                                                                                     aCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGgag                                                             >  1:1275337/1‑150 (MQ=255)
                                                                                      cAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGgaga                                                            <  1:5308283/150‑1 (MQ=255)
                                                                                       aaaaGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGagaa                                                           >  1:3084491/1‑150 (MQ=255)
                                                                                       aaaaGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGagaa                                                           >  1:1306472/1‑150 (MQ=255)
                                                                                               ccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGAt                                                   >  1:3071906/1‑150 (MQ=255)
                                                                                                 aTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATaa                                                 >  2:522461/1‑150 (MQ=255)
                                                                                                  tCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATaaa                                                <  2:1340994/150‑1 (MQ=255)
                                                                                                   cGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAg                                               >  2:3639495/1‑150 (MQ=255)
                                                                                                     gaacAACGCTTTCATCTATTTCAACGGTGATTATTATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGTCGGCTTCGCGCAGCTCGTTCATCTGGTTTTTCAATGAGATAAAGCCAGGCGCGTTCATGGTGAAGTCCTGATAAAGaa                                             <  2:4634476/147‑1 (MQ=255)
                                                                                                     gCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGaa                                             <  2:2988455/150‑1 (MQ=255)
                                                                                                     gCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGaa                                             >  2:3562719/1‑150 (MQ=255)
                                                                                                     gCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGaa                                             >  2:2574113/1‑150 (MQ=255)
                                                                                                      cccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAAt                                            <  2:1304401/150‑1 (MQ=255)
                                                                                                       ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                           >  2:3223971/1‑150 (MQ=255)
                                                                                                        cAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATgg                                          >  2:1797153/1‑150 (MQ=255)
                                                                                                        cAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATgg                                          >  2:3718019/1‑150 (MQ=255)
                                                                                                            gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                      <  2:3080192/150‑1 (MQ=255)
                                                                                                            gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCAGAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                      <  2:2704023/150‑1 (MQ=255)
                                                                                                               gACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcac                                   >  1:1020606/1‑150 (MQ=255)
                                                                                                               gACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcac                                   <  1:899184/150‑1 (MQ=255)
                                                                                                                aCATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcacc                                  <  1:319130/150‑1 (MQ=255)
                                                                                                                aCATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcacc                                  <  2:4972803/150‑1 (MQ=255)
                                                                                                                aCATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcacc                                  >  1:4082942/1‑150 (MQ=255)
                                                                                                                   tCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCaccacc                               <  2:1661213/150‑1 (MQ=255)
                                                                                                                      aTTTCACCGGAGATAATATTGAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTAATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGATTGGCAATCAGCACCAGc                            <  2:101265/150‑1 (MQ=255)
                                                                                                                      aTTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGc                            >  1:2828583/1‑150 (MQ=255)
                                                                                                                         tCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCaaa                         >  1:443532/1‑150 (MQ=255)
                                                                                                                              ggAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAt                    <  1:3284657/150‑1 (MQ=255)
                                                                                                                                agaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                  >  1:2397385/1‑150 (MQ=255)
                                                                                                                                agaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                  <  2:3190635/150‑1 (MQ=255)
                                                                                                                                agaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                  >  1:767678/1‑150 (MQ=255)
                                                                                                                                agaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                  >  1:5230086/1‑150 (MQ=255)
                                                                                                                                agaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                  >  2:4397961/1‑150 (MQ=255)
                                                                                                                                 gaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtctc                 >  2:4019978/1‑150 (MQ=255)
                                                                                                                                   taataaTCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAt               >  1:4764376/1‑150 (MQ=255)
                                                                                                                                    aataatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAtt              >  1:883506/1‑150 (MQ=255)
                                                                                                                                     ataatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTa             >  1:643579/1‑150 (MQ=255)
                                                                                                                                     ataatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTa             >  1:2044677/1‑150 (MQ=255)
                                                                                                                                      taatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCACTCACCACCAGCAAACCGATCTCATTaa            >  2:2820800/1‑150 (MQ=255)
                                                                                                                                       aatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAat           <  1:3733126/150‑1 (MQ=255)
                                                                                                                                         tCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAatat         <  2:3465275/150‑1 (MQ=255)
                                                                                                                                         tCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAatat         <  2:1547053/150‑1 (MQ=255)
                                                                                                                                           aGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg       >  2:4247230/1‑150 (MQ=255)
                                                                                                                                           aGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg       >  1:853690/1‑150 (MQ=255)
                                                                                                                                           aGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg       >  1:968788/1‑150 (MQ=255)
                                                                                                                                                cATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACt  <  1:4220152/150‑1 (MQ=255)
                                                                                                                                                cATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACt  >  1:2346688/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACT  >  CP009273/2669831‑2670124

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: