Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,669,978 C→A 42.8% 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,9780CA42.8% 325.9 / inf 345A69A (GCG→GCTyphDinner membrane putative ABC superfamily sugar transporter permease
Reads supporting (aligned to +/- strand):  ref base C (102/95);  new base A (80/67);  total (182/163)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

CCATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGAC  >  CP009273/2669830‑2670127
                                                                                                                                                    |                                                                                                                                                     
    cAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGACCa                                                                                                                                                  >  2:1766157/1‑150 (MQ=255)
     aGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGTGTCATCGCCCAg                                                                                                                                                 <  2:3238898/150‑1 (MQ=255)
     aGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAg                                                                                                                                                 >  1:1201807/1‑150 (MQ=255)
     aGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAg                                                                                                                                                 <  2:4582992/150‑1 (MQ=255)
        gcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAggcg                                                                                                                                              <  1:3323230/150‑1 (MQ=255)
             tAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAAt                                                                                                                                         <  2:1424772/150‑1 (MQ=255)
             tAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAAATTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAAt                                                                                                                                         >  2:1662697/1‑150 (MQ=255)
              aaCAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATc                                                                                                                                        >  2:4231541/1‑150 (MQ=255)
                 aGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCa                                                                                                                                     <  2:4745347/150‑1 (MQ=255)
                    aaCACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtg                                                                                                                                  >  1:4592387/1‑150 (MQ=255)
                      cacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtggt                                                                                                                                <  2:1508260/150‑1 (MQ=255)
                            cagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTggcggc                                                                                                                          >  2:3974451/1‑150 (MQ=255)
                              ggCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCAt                                                                                                                        <  2:3219469/150‑1 (MQ=255)
                              ggCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGAACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCACTCCCAATGGTGGCGGCAt                                                                                                                        >  2:830806/1‑150 (MQ=255)
                                cACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATcg                                                                                                                      <  1:4115201/150‑1 (MQ=255)
                                   gCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCgcgc                                                                                                                   >  2:1225056/1‑150 (MQ=255)
                                     cACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgcag                                                                                                                 >  2:18703/1‑150 (MQ=255)
                                       ccGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgcagca                                                                                                               <  1:1538566/150‑1 (MQ=255)
                                        cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                              >  1:828205/1‑150 (MQ=255)
                                            aGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTc                                                                                                          <  1:2111639/150‑1 (MQ=255)
                                             gCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGAGCAGCACGTTCa                                                                                                         >  2:2564091/1‑150 (MQ=255)
                                                gAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCt                                                                                                      <  1:3442117/150‑1 (MQ=255)
                                                 aaCTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTg                                                                                                     <  1:4527820/150‑1 (MQ=255)
                                                  aCTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTgg                                                                                                    >  1:3030135/1‑150 (MQ=255)
                                                    ttCAAAGTGCAGCACACATGCCAGGCACACCGAGACAAAAGCCACCAGCGGCCCAACGCTGACATCAATTTCACCGGAGATAAGAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGtt                                                                                                  <  2:545868/150‑1 (MQ=255)
                                                      cAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTAt                                                                                                <  2:2631399/150‑1 (MQ=255)
                                                      cAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTAt                                                                                                <  2:261658/150‑1 (MQ=255)
                                                         aTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCa                                                                                             >  2:1651678/1‑150 (MQ=255)
                                                         aTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCa                                                                                             <  1:635336/150‑1 (MQ=255)
                                                         aTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCa                                                                                             >  2:2194226/1‑150 (MQ=255)
                                                         aTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCa                                                                                             <  2:1994706/150‑1 (MQ=255)
                                                           tGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAAt                                                                                           <  2:4601156/150‑1 (MQ=255)
                                                           tGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAAt                                                                                           <  1:4034885/150‑1 (MQ=255)
                                                               gcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATgaga                                                                                       <  1:1971567/150‑1 (MQ=255)
                                                                 aaaaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATa                                                                                     >  1:591522/1‑150 (MQ=255)
                                                                   aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                   >  1:4295730/1‑150 (MQ=255)
                                                                   aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                   >  2:3862951/1‑150 (MQ=255)
                                                                   aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                   <  2:3498764/150‑1 (MQ=255)
                                                                    aaTGCAAGGCACACCGAGACAAAAGACACCATCGGCCCAACTCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATTGTGTCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAg                                                                                  <  2:4345919/150‑1 (MQ=255)
                                                                     aTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGc                                                                                 >  1:4515031/1‑150 (MQ=255)
                                                                     aTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGc                                                                                 >  1:4680856/1‑150 (MQ=255)
                                                                       gCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCa                                                                               <  1:2662388/150‑1 (MQ=255)
                                                                        ccAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAg                                                                              >  2:3099802/1‑150 (MQ=255)
                                                                        ccAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAg                                                                              >  2:4144784/1‑150 (MQ=255)
                                                                         cAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAgg                                                                             >  2:3714227/1‑150 (MQ=255)
                                                                         cAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAgg                                                                             <  1:4155157/150‑1 (MQ=255)
                                                                          aGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgc                                                                            >  1:4580975/1‑150 (MQ=255)
                                                                          aGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgc                                                                            >  2:4694094/1‑150 (MQ=255)
                                                                           ggCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgcg                                                                           <  1:722009/150‑1 (MQ=255)
                                                                                  cGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAg                                                                    <  2:2777392/150‑1 (MQ=255)
                                                                                   gagaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAgg                                                                   <  1:2117750/150‑1 (MQ=255)
                                                                                   gagaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAgg                                                                   >  2:3789277/1‑150 (MQ=255)
                                                                                     gaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGga                                                                 <  2:1971005/150‑1 (MQ=255)
                                                                                      aCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGgag                                                                >  2:1172815/1‑150 (MQ=255)
                                                                                      aCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGgag                                                                >  2:2975670/1‑150 (MQ=255)
                                                                                           aGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGAc                                                           <  2:69884/150‑1 (MQ=255)
                                                                                           aGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGAc                                                           >  2:72453/1‑150 (MQ=255)
                                                                                            gCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGAcc                                                          >  1:2210309/1‑150 (MQ=255)
                                                                                             ccaccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCa                                                         >  1:4789681/1‑150 (MQ=255)
                                                                                             ccaccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCa                                                         <  1:4376242/150‑1 (MQ=255)
                                                                                                ccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGAt                                                      >  2:3291052/1‑150 (MQ=255)
                                                                                                ccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGAt                                                      >  2:3340916/1‑150 (MQ=255)
                                                                                                 caTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATa                                                     >  1:2358315/1‑150 (MQ=255)
                                                                                                 caTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATa                                                     >  2:861861/1‑150 (MQ=255)
                                                                                                     ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                                 <  2:412083/150‑1 (MQ=255)
                                                                                                     ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                                 <  2:3378290/150‑1 (MQ=255)
                                                                                                      gCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGaa                                                >  1:2842522/1‑150 (MQ=255)
                                                                                                        ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                              <  1:3889425/150‑1 (MQ=255)
                                                                                                         cAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCACTGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATgg                                             <  2:4064033/150‑1 (MQ=255)
                                                                                                         cAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATgg                                             <  2:2434648/150‑1 (MQ=255)
                                                                                                         cAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATgg                                             <  1:109308/150‑1 (MQ=255)
                                                                                                          aaCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGc                                            >  2:625488/1‑150 (MQ=255)
                                                                                                          aaCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGc                                            <  1:3655925/150‑1 (MQ=255)
                                                                                                            cGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCaa                                          >  1:4211302/1‑150 (MQ=255)
                                                                                                             gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                         <  1:333372/150‑1 (MQ=255)
                                                                                                             gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                         >  2:3482498/1‑150 (MQ=255)
                                                                                                                gACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcac                                      >  1:2334199/1‑150 (MQ=255)
                                                                                                                   aTCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcaccac                                   >  1:3949661/1‑150 (MQ=255)
                                                                                                                        tttCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCa                              >  1:248708/1‑150 (MQ=255)
                                                                                                                        tttCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCa                              >  2:4199224/1‑150 (MQ=255)
                                                                                                                          tCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCaaa                            >  2:1347304/1‑150 (MQ=255)
                                                                                                                          tCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCaaa                            >  2:3532013/1‑150 (MQ=255)
                                                                                                                           cACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAAc                           <  1:3769225/150‑1 (MQ=255)
                                                                                                                                   ataataATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCa                   >  1:1419562/1‑150 (MQ=255)
                                                                                                                                   ataataATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCa                   >  2:2030566/1‑150 (MQ=255)
                                                                                                                                   ataataATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCa                   <  2:3120978/150‑1 (MQ=255)
                                                                                                                                     aataatCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAtt                 >  2:1295903/1‑150 (MQ=255)
                                                                                                                                      ataatCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTa                >  2:1066307/1‑150 (MQ=255)
                                                                                                                                       taatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTaa               <  2:2260209/150‑1 (MQ=255)
                                                                                                                                          tCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAatat            <  1:552939/150‑1 (MQ=255)
                                                                                                                                          tCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAatat            >  2:1523852/1‑150 (MQ=255)
                                                                                                                                           cAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATg           >  2:2574860/1‑150 (MQ=255)
                                                                                                                                            aGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg          >  1:172071/1‑150 (MQ=255)
                                                                                                                                            aGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATgg          >  2:240001/1‑150 (MQ=255)
                                                                                                                                              cGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCg        >  1:1048995/1‑150 (MQ=255)
                                                                                                                                              cGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCg        >  1:3149139/1‑150 (MQ=255)
                                                                                                                                                tCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGAc      >  2:2884170/1‑150 (MQ=255)
                                                                                                                                                 cATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATTGCGACt     <  2:170764/150‑1 (MQ=255)
                                                                                                                                                   tCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGa   <  2:4446160/150‑1 (MQ=255)
                                                                                                                                                    cGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGAc  >  2:74701/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACTGAC  >  CP009273/2669830‑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: