Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,669,978 C→A 46.7% 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,9780CA46.7% 182.0 / inf 286A69A (GCG→GCTyphDinner membrane putative ABC superfamily sugar transporter permease
Reads supporting (aligned to +/- strand):  ref base C (55/97);  new base A (81/53);  total (136/150)
Fisher's exact test for biased strand distribution p-value = 5.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

ATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGAC  >  CP009273/2669832‑2670127
                                                                                                                                                   |                                                                                                                                                     
   aGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAg                                                                                                                                                 <  1:730127/150‑1 (MQ=255)
   aGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAg                                                                                                                                                 >  1:2048233/1‑150 (MQ=255)
    gcgcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAgg                                                                                                                                                >  1:673301/1‑150 (MQ=255)
         ccTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCa                                                                                                                                           <  1:1163384/150‑1 (MQ=255)
         ccTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCa                                                                                                                                           >  1:527203/1‑150 (MQ=255)
             aCAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATcc                                                                                                                                       <  1:3225779/150‑1 (MQ=255)
               aGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCa                                                                                                                                     >  1:3774247/1‑150 (MQ=255)
                   acacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtgg                                                                                                                                 >  1:1031569/1‑150 (MQ=255)
                      ccAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACTCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtggtgg                                                                                                                              <  2:2422537/150‑1 (MQ=255)
                       cagcagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTggc                                                                                                                             >  1:1886094/1‑150 (MQ=255)
                             gCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATc                                                                                                                       <  1:797942/150‑1 (MQ=255)
                             gCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATc                                                                                                                       <  1:959912/150‑1 (MQ=255)
                             gCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATc                                                                                                                       <  1:4158536/150‑1 (MQ=255)
                             gCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATc                                                                                                                       <  1:3011635/150‑1 (MQ=255)
                                   ccccgcccgcCGAACTTCAACTTGCCGCAAAAATTCCCCGCACCCCGACTCAAAACGCAACCTCCGCCCAACGGCTTAATCAATTTACACGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgcag                                                                                                                 <  2:373523/148‑1 (MQ=255)
                                    aCCGCCAGAGGAACTTCAAATTTCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCAAAGGCCCAAAGCTGACATCAAAATCACCGGAGATAAAAATCATCGTCATCGCCCAGGCGGCAATCACAAACGTGGCGGCATCGCgcagc                                                                                                                >  1:4235545/1‑150 (MQ=255)
                                      cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                              <  2:4261302/150‑1 (MQ=255)
                                       gCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACg                                                                                                             <  2:3053738/150‑1 (MQ=255)
                                       gCCAGCGGAACTTCAAAGTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACg                                                                                                             <  2:2778335/150‑1 (MQ=255)
                                         cAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGtt                                                                                                           <  2:2825416/150‑1 (MQ=255)
                                         cAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGtt                                                                                                           <  2:612665/150‑1 (MQ=255)
                                            cGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATACTAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCAt                                                                                                        <  2:626247/150‑1 (MQ=255)
                                            cGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCAt                                                                                                        <  1:276806/150‑1 (MQ=255)
                                            cGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCAt                                                                                                        >  2:3502536/1‑150 (MQ=255)
                                                 cTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGt                                                                                                   >  2:3029294/1‑150 (MQ=255)
                                                 cTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCACAGCACGTTCATCTGGt                                                                                                   >  1:3532465/1‑150 (MQ=255)
                                                  ttCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGtt                                                                                                  >  1:1802352/1‑150 (MQ=255)
                                                  ttCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGtt                                                                                                  >  1:3952384/1‑150 (MQ=255)
                                                    cAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTAt                                                                                                <  2:1700635/150‑1 (MQ=255)
                                                    cAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTAt                                                                                                >  1:3088949/1‑150 (MQ=255)
                                                       aTTGCAGCAAACATGCCAGGCACCCCGAGACAACAGCCCCCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCa                                                                                             <  2:3311029/150‑1 (MQ=255)
                                                              caAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGAt                                                                                      >  1:2320591/1‑150 (MQ=255)
                                                              caAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGAt                                                                                      <  1:4106392/150‑1 (MQ=255)
                                                               aaaaaTGCCATGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGAAATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATa                                                                                     <  2:4235545/150‑1 (MQ=255)
                                                                aaaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaa                                                                                    <  1:4222345/150‑1 (MQ=255)
                                                                 aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                   >  2:2055997/1‑150 (MQ=255)
                                                                 aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                   <  2:1014402/150‑1 (MQ=255)
                                                                 aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                   >  2:2162868/1‑150 (MQ=255)
                                                                  aaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAg                                                                                  >  2:922887/1‑150 (MQ=255)
                                                                   aTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGc                                                                                 <  1:2379205/150‑1 (MQ=255)
                                                                     gCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCa                                                                               <  1:4023409/150‑1 (MQ=255)
                                                                       cAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCGGGTTATTCAATGAGATAAAGCCAgg                                                                             >  2:850859/1‑150 (MQ=255)
                                                                        aGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgc                                                                            >  2:2628157/1‑150 (MQ=255)
                                                                            acacCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGt                                                                        >  2:2710536/1‑150 (MQ=255)
                                                                               ccGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCa                                                                     <  1:436388/150‑1 (MQ=255)
                                                                                cGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAg                                                                    <  2:1314737/150‑1 (MQ=255)
                                                                                cGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAg                                                                    <  1:4180169/150‑1 (MQ=255)
                                                                                cGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAg                                                                    <  1:3356652/150‑1 (MQ=255)
                                                                                 gagaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAgg                                                                   >  2:3178459/1‑150 (MQ=255)
                                                                                 gagaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAgg                                                                   >  2:2436186/1‑150 (MQ=255)
                                                                                   gaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGga                                                                 <  2:755250/150‑1 (MQ=255)
                                                                                    aCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGgag                                                                <  1:3837066/150‑1 (MQ=255)
                                                                                      aaaaGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCCATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGagaa                                                              <  1:2198940/150‑1 (MQ=255)
                                                                                      aaaaGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGagaa                                                              >  1:3548605/1‑150 (MQ=255)
                                                                                       aaaGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGagaag                                                             >  2:604177/1‑150 (MQ=255)
                                                                                           ccaccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCa                                                         >  1:1614267/1‑150 (MQ=255)
                                                                                            caccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAGGCCAGGCGCGTTCAGGGAGAAGACCAg                                                        <  2:1429804/150‑1 (MQ=255)
                                                                                                   ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                                 <  1:3291727/150‑1 (MQ=255)
                                                                                                   ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                                 >  1:345194/1‑150 (MQ=255)
                                                                                                    gCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGaa                                                <  1:4243858/150‑1 (MQ=255)
                                                                                                      ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                              >  2:351377/1‑150 (MQ=255)
                                                                                                      ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                              <  1:1246530/150‑1 (MQ=255)
                                                                                                      ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                              <  1:905752/150‑1 (MQ=255)
                                                                                                       cAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATgg                                             <  1:1699724/150‑1 (MQ=255)
                                                                                                       cAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCGTCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATgg                                             >  1:2011672/1‑150 (MQ=255)
                                                                                                        aaCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGc                                            <  2:436341/150‑1 (MQ=255)
                                                                                                         aCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCa                                           <  2:137715/150‑1 (MQ=255)
                                                                                                         aCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCa                                           <  2:4028094/150‑1 (MQ=255)
                                                                                                            gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCACTGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                         <  2:2473978/150‑1 (MQ=255)
                                                                                                            gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                         <  1:2566709/150‑1 (MQ=255)
                                                                                                             cTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATc                                        >  2:362871/1‑150 (MQ=255)
                                                                                                              tGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCACTGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATca                                       <  1:3198997/150‑1 (MQ=255)
                                                                                                              tGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATca                                       >  1:1324823/1‑150 (MQ=255)
                                                                                                              tGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATca                                       >  1:3518386/1‑150 (MQ=255)
                                                                                                               gACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcac                                      <  1:511151/150‑1 (MQ=255)
                                                                                                                 cATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcacca                                    <  2:2484965/150‑1 (MQ=255)
                                                                                                                  aTCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcaccac                                   >  2:1567175/1‑150 (MQ=255)
                                                                                                                    cAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCAccacca                                 <  2:1936737/150‑1 (MQ=255)
                                                                                                                     aaTTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAg                                >  2:935036/1‑150 (MQ=255)
                                                                                                                      aTTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGc                               >  2:1842465/1‑150 (MQ=255)
                                                                                                                      aTTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGc                               >  2:3446291/1‑150 (MQ=255)
                                                                                                                        ttCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCaa                             >  1:4206971/1‑150 (MQ=255)
                                                                                                                        ttCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCaa                             >  1:3701779/1‑150 (MQ=255)
                                                                                                                         tCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCaaa                            >  2:3822989/1‑150 (MQ=255)
                                                                                                                           aCCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCCCCACCAGCAAAcc                          >  2:4189435/1‑150 (MQ=255)
                                                                                                                           aCCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAAcc                          >  2:3264904/1‑150 (MQ=255)
                                                                                                                             cGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGa                        <  1:604177/150‑1 (MQ=255)
                                                                                                                                agaTAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGATAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                     <  2:370726/150‑1 (MQ=255)
                                                                                                                                agaTAATAATCAGCGTCAGAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                     <  2:3225273/150‑1 (MQ=255)
                                                                                                                                    aataatCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAtt                 <  2:1795561/150‑1 (MQ=255)
                                                                                                                                     ataatCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTa                <  2:1618397/150‑1 (MQ=255)
                                                                                                                                     ataatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTa                <  2:2320591/150‑1 (MQ=255)
                                                                                                                                          cAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATg           >  1:1462202/1‑150 (MQ=255)
                                                                                                                                           agatacATAGCCCAGGCGGCACTCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGTGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAAACGATCTCATTAATATgg          >  2:2594439/6‑150 (MQ=255)
                                                                                                                                           aGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg          <  2:2083288/150‑1 (MQ=255)
                                                                                                                                           aGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg          >  1:831297/1‑150 (MQ=255)
                                                                                                                                             cGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCg        <  1:1969249/150‑1 (MQ=255)
                                                                                                                                              gTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGa       >  2:4341683/1‑150 (MQ=255)
                                                                                                                                              gTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGa       >  1:3741997/1‑150 (MQ=255)
                                                                                                                                              gTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGa       >  2:3968164/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGAC  >  CP009273/2669832‑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: