Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,669,978 C→A 60.5% 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,9780CA60.5% 82.1 / inf 251A69A (GCG→GCTyphDinner membrane putative ABC superfamily sugar transporter permease
Reads supporting (aligned to +/- strand):  ref base C (50/49);  new base A (74/78);  total (124/127)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 252 total aligned reads displayed.

CCCATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACT  >  CP009273/2669829‑2670124
                                                                                                                                                     |                                                                                                                                                  
cccATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATc                                                                                                                                                    <  2:2595142/150‑1 (MQ=255)
    tCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGccc                                                                                                                                                <  1:509850/150‑1 (MQ=255)
        cgcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATTCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAggc                                                                                                                                            <  2:591486/150‑1 (MQ=255)
         gcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggcg                                                                                                                                           >  1:4212655/1‑150 (MQ=255)
          cgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAggcgg                                                                                                                                          >  2:1508764/1‑150 (MQ=255)
            ccTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCa                                                                                                                                        <  1:2370111/150‑1 (MQ=255)
               aaCAGCACCACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATc                                                                                                                                     <  2:3026811/150‑1 (MQ=255)
               aaCAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATc                                                                                                                                     >  1:2251136/1‑150 (MQ=255)
               aaCAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATc                                                                                                                                     <  2:4524948/150‑1 (MQ=255)
                aCAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATcc                                                                                                                                    >  2:215605/1‑150 (MQ=255)
                  cgCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAACATGCCAGGCACCCCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCa                                                                                                                                  <  2:691007/149‑1 (MQ=255)
                     aaCACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAtg                                                                                                                               <  2:1034845/150‑1 (MQ=255)
                       cacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAAtggt                                                                                                                             >  2:3771773/1‑150 (MQ=255)
                          cagcagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTggc                                                                                                                          <  2:4078247/150‑1 (MQ=255)
                            gcagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTggcgg                                                                                                                        <  1:23135/150‑1 (MQ=255)
                              agGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCa                                                                                                                      >  2:2107966/1‑150 (MQ=255)
                                 cACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGTCGGCATcg                                                                                                                   <  2:4024428/150‑1 (MQ=255)
                                      cACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgcag                                                                                                              >  1:2520700/1‑150 (MQ=255)
                                         cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                           >  2:3588029/1‑150 (MQ=255)
                                         cGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCAc                                                                                                           >  2:3588075/1‑150 (MQ=255)
                                          gCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACg                                                                                                          >  1:416979/1‑150 (MQ=255)
                                           ccAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGt                                                                                                         >  1:363082/1‑150 (MQ=255)
                                            ccgCGGAACTTCAAATAGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGtt                                                                                                        <  2:4134530/148‑1 (MQ=255)
                                              gCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCa                                                                                                      >  1:753846/1‑150 (MQ=255)
                                                 gAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCt                                                                                                   <  1:4122273/150‑1 (MQ=255)
                                                  aaCTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTg                                                                                                  >  1:2174509/1‑150 (MQ=255)
                                                      tCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTa                                                                                              >  2:340570/1‑150 (MQ=255)
                                                         aaTTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTc                                                                                           >  1:2593546/1‑150 (MQ=255)
                                                           ctGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCaa                                                                                         <  2:4299876/149‑1 (MQ=255)
                                                            tGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAAt                                                                                        <  2:636250/150‑1 (MQ=255)
                                                            tGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAAt                                                                                        <  1:1318836/150‑1 (MQ=255)
                                                               agcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATTGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATgag                                                                                     <  1:1087202/150‑1 (MQ=255)
                                                                   aaaaGGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaa                                                                                 <  2:3079164/150‑1 (MQ=255)
                                                                    aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                >  2:3399981/1‑150 (MQ=255)
                                                                    aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                <  2:3641390/150‑1 (MQ=255)
                                                                    aaaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATaaa                                                                                <  2:1012027/150‑1 (MQ=255)
                                                                     aaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAg                                                                               <  1:1322772/150‑1 (MQ=255)
                                                                     aaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAg                                                                               <  1:642819/150‑1 (MQ=255)
                                                                     aaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAg                                                                               <  1:1254993/150‑1 (MQ=255)
                                                                     aaTGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAg                                                                               >  1:42738/1‑150 (MQ=255)
                                                                       tGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGcc                                                                             <  2:366338/150‑1 (MQ=255)
                                                                       tGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGcc                                                                             >  1:2289777/1‑150 (MQ=255)
                                                                        gCCAGGCACACCGAGACAACAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGATGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCa                                                                            <  2:4340126/150‑1 (MQ=255)
                                                                        gCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCa                                                                            >  2:2696942/1‑150 (MQ=255)
                                                                        gCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCa                                                                            >  2:2965424/1‑150 (MQ=255)
                                                                        gCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCa                                                                            <  1:3855380/150‑1 (MQ=255)
                                                                           aGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGgc                                                                         >  2:2559232/1‑150 (MQ=255)
                                                                               acacCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGt                                                                     >  1:3308733/1‑150 (MQ=255)
                                                                                cacCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGtt                                                                    >  2:541333/1‑150 (MQ=255)
                                                                                cacCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGtt                                                                    <  1:3665234/150‑1 (MQ=255)
                                                                                   cGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAg                                                                 >  1:4060398/1‑150 (MQ=255)
                                                                                     agaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAggg                                                               <  1:3275562/150‑1 (MQ=255)
                                                                                      gaCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGga                                                              <  1:652550/150‑1 (MQ=255)
                                                                                       aCAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGgag                                                             <  1:4386120/150‑1 (MQ=255)
                                                                                         aaaaGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGagaa                                                           >  1:3263915/1‑150 (MQ=255)
                                                                                             gCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGAcc                                                       <  2:2228980/150‑1 (MQ=255)
                                                                                                 ccaTCGGCCCAACGCTGACATCTATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGAGGCATCTCGAAGCACGTTCTTCTGGTTATTCAATGTGATAATGCCAGGCGCGTTCAGGGAGAAGACCAGAt                                                   <  1:1875676/150‑1 (MQ=255)
                                                                                                 ccaTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGAt                                                   <  1:2148761/150‑1 (MQ=255)
                                                                                                   aTCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATaa                                                 >  2:417591/1‑150 (MQ=255)
                                                                                                     cGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAg                                               <  2:2199074/150‑1 (MQ=255)
                                                                                                      ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                              <  1:199679/150‑1 (MQ=255)
                                                                                                      ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                              <  1:3191371/150‑1 (MQ=255)
                                                                                                      ggCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGa                                              >  1:2418236/1‑150 (MQ=255)
                                                                                                        cccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAAt                                            <  1:3726004/150‑1 (MQ=255)
                                                                                                         ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                           <  1:1607163/150‑1 (MQ=255)
                                                                                                         ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATACAGCCAGGCGCGGTCAGGGAGAAGACCAGATAAAGAATg                                           >  2:2508490/1‑150 (MQ=255)
                                                                                                         ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGACTg                                           >  2:4362550/1‑150 (MQ=255)
                                                                                                         ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                           >  1:307201/1‑150 (MQ=255)
                                                                                                         ccAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATg                                           <  2:4333238/150‑1 (MQ=255)
                                                                                                           aaCGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGc                                         <  1:3001628/150‑1 (MQ=255)
                                                                                                             cGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCaa                                       >  2:1075695/1‑150 (MQ=255)
                                                                                                              gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                      <  1:827067/150‑1 (MQ=255)
                                                                                                              gCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAAt                                      >  1:151092/1‑150 (MQ=255)
                                                                                                               cTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATc                                     <  2:1370967/150‑1 (MQ=255)
                                                                                                                tGACATCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATca                                    >  2:2738199/1‑150 (MQ=255)
                                                                                                                  aCATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATcacc                                  >  2:1724882/1‑150 (MQ=255)
                                                                                                                     tCAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGAAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCaccacc                               >  2:773008/1‑150 (MQ=255)
                                                                                                                      cAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCAccacca                              >  1:29879/1‑150 (MQ=255)
                                                                                                                      cAATTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCAccacca                              <  2:2251136/150‑1 (MQ=255)
                                                                                                                       aaTTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAg                             <  1:2083431/150‑1 (MQ=255)
                                                                                                                       aaTTTCACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATTAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATAACAACCAg                             >  2:3751765/1‑150 (MQ=255)
                                                                                                                         tttCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCa                           <  1:3138314/150‑1 (MQ=255)
                                                                                                                            cACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAAc                        <  1:418905/150‑1 (MQ=255)
                                                                                                                            cACCGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAAc                        >  2:1788677/1‑150 (MQ=255)
                                                                                                                              ccGGAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCg                      >  1:4181625/1‑150 (MQ=255)
                                                                                                                                ggAGATAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAt                    <  2:1911010/150‑1 (MQ=255)
                                                                                                                                  agaTAATAATCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGAtct                  >  2:2535940/1‑150 (MQ=255)
                                                                                                                                    ataataATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCa                <  1:2784715/150‑1 (MQ=255)
                                                                                                                                      aataatCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAtt              <  1:1279036/150‑1 (MQ=255)
                                                                                                                                      aataatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCAtt              <  1:4300291/150‑1 (MQ=255)
                                                                                                                                        taatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTaa            <  1:4534845/150‑1 (MQ=255)
                                                                                                                                        taatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTaa            <  1:3449476/150‑1 (MQ=255)
                                                                                                                                        taatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTaa            >  2:2579119/1‑150 (MQ=255)
                                                                                                                                         aatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAat           >  2:3197718/1‑150 (MQ=255)
                                                                                                                                         aatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAat           <  2:537086/150‑1 (MQ=255)
                                                                                                                                         aatCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAat           >  1:2355586/1‑150 (MQ=255)
                                                                                                                                           tCAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAatat         <  1:1027254/150‑1 (MQ=255)
                                                                                                                                            cAGCGTCATAGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATg        >  1:2699/1‑150 (MQ=255)
                                                                                                                                               cGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCg     >  2:2875197/1‑150 (MQ=255)
                                                                                                                                                  cATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACt  >  2:240860/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCCATCAGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCATCTGGTTATTCAATGAGATAAAGCCAGGCGCGTTCAGGGAGAAGACCAGATAAAGAATGGCAATCACCACCAGCAAACCGATCTCATTAATATGGCGACT  >  CP009273/2669829‑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: