Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 104,561 G→T 93.6% A124S (GCG→TCG)  secM → regulator of secA translation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273104,5610GT93.6% 944.9 / 45.2 281A124S (GCG→TCG) secMregulator of secA translation
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base T (138/125);  total (148/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

AACTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAAC  >  CP009273/104414‑104705
                                                                                                                                                   |                                                                                                                                                
aaCTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAgctcg                                                                                                                                                <  1:1070294/150‑1 (MQ=255)
  cTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCgct                                                                                                                                              >  1:3314712/1‑150 (MQ=255)
  cTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCgct                                                                                                                                              <  1:4665780/150‑1 (MQ=255)
    aTTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCgctgc                                                                                                                                            >  2:3748725/1‑150 (MQ=255)
     ttCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCgctgct                                                                                                                                           >  2:2411433/1‑150 (MQ=255)
     ttCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCgctgct                                                                                                                                           <  1:406266/150‑1 (MQ=255)
       ccGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGa                                                                                                                                         <  1:4060061/150‑1 (MQ=255)
        cGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGAc                                                                                                                                        <  1:817672/150‑1 (MQ=255)
        cGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGAc                                                                                                                                        <  2:2590961/150‑1 (MQ=255)
           tGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCa                                                                                                                                     >  2:3606400/1‑150 (MQ=255)
           tGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCACAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCa                                                                                                                                     <  2:1422167/150‑1 (MQ=255)
             aTTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAgg                                                                                                                                   >  1:390060/1‑150 (MQ=255)
              ttACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGa                                                                                                                                  <  1:1696349/150‑1 (MQ=255)
              ttACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGa                                                                                                                                  <  1:369633/150‑1 (MQ=255)
                aCTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAg                                                                                                                                <  2:2018340/150‑1 (MQ=255)
                 cTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAgg                                                                                                                               <  2:1302077/150‑1 (MQ=255)
                 cTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAgg                                                                                                                               >  2:2139651/1‑150 (MQ=255)
                     cATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACg                                                                                                                           >  1:2983228/1‑150 (MQ=255)
                     cATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACg                                                                                                                           >  2:3474910/1‑150 (MQ=255)
                     cATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACg                                                                                                                           <  1:4632737/150‑1 (MQ=255)
                       tCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGcc                                                                                                                         >  2:250424/1‑150 (MQ=255)
                       tCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGcc                                                                                                                         <  1:427380/150‑1 (MQ=255)
                          aCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTc                                                                                                                      >  1:1760039/1‑150 (MQ=255)
                           cATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCt                                                                                                                     <  1:2074609/150‑1 (MQ=255)
                             tGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGa                                                                                                                   <  1:4092049/150‑1 (MQ=255)
                             tGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGa                                                                                                                   >  1:2873789/1‑150 (MQ=255)
                             tGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGa                                                                                                                   <  2:4648954/150‑1 (MQ=255)
                             tGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGa                                                                                                                   <  1:1555873/150‑1 (MQ=255)
                               ccATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaaa                                                                                                                 >  2:3204142/1‑150 (MQ=255)
                                cATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaaaa                                                                                                                <  1:3459137/150‑1 (MQ=255)
                                cATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaaaa                                                                                                                >  2:3457900/1‑150 (MQ=255)
                                 aTTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAg                                                                                                               >  2:3804690/1‑150 (MQ=255)
                                 aTTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAg                                                                                                               >  1:2749961/1‑150 (MQ=255)
                                  ttCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAgg                                                                                                              >  1:2438283/1‑150 (MQ=255)
                                    cGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGt                                                                                                            <  1:255192/150‑1 (MQ=255)
                                     gCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGtt                                                                                                           <  2:3753893/150‑1 (MQ=255)
                                       aCGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTAt                                                                                                         >  1:1154936/1‑150 (MQ=255)
                                        cGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATc                                                                                                        >  2:1401789/1‑150 (MQ=255)
                                            aaTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAt                                                                                                    <  2:4360982/150‑1 (MQ=255)
                                             aTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAtt                                                                                                   <  2:4003421/150‑1 (MQ=255)
                                             aTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAtt                                                                                                   <  2:3071228/150‑1 (MQ=255)
                                             aTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAtt                                                                                                   <  1:4407607/150‑1 (MQ=255)
                                              tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                  >  2:14866/1‑150 (MQ=255)
                                              tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                  >  2:3320249/1‑150 (MQ=255)
                                               ccGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGa                                                                                                 <  2:2448044/150‑1 (MQ=255)
                                                 gTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGAtt                                                                                               >  1:2005720/1‑150 (MQ=255)
                                                  tCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTa                                                                                              <  1:4577436/150‑1 (MQ=255)
                                                   cATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTAt                                                                                             <  1:592568/150‑1 (MQ=255)
                                                    aTCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAGTGAGTATg                                                                                            >  1:4382979/1‑150 (MQ=255)
                                                       tttctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATgcgc                                                                                         >  2:2385445/1‑150 (MQ=255)
                                                       tttctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAGTGATGATgcgc                                                                                         >  1:2360026/1‑150 (MQ=255)
                                                         tctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAt                                                                                       >  1:3099296/1‑150 (MQ=255)
                                                            ttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtttt                                                                                    <  2:817378/150‑1 (MQ=255)
                                                            ttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtttt                                                                                    <  1:1000642/150‑1 (MQ=255)
                                                               gCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAcc                                                                                 >  2:4315695/1‑150 (MQ=255)
                                                                cAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAccc                                                                                >  2:952673/1‑150 (MQ=255)
                                                                cAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAccc                                                                                >  1:3647353/1‑150 (MQ=255)
                                                                    ggCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCcaca                                                                            >  2:3532811/1‑150 (MQ=255)
                                                                      cACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAg                                                                          >  1:1009170/1‑150 (MQ=255)
                                                                      cACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAg                                                                          >  2:2121868/1‑150 (MQ=255)
                                                                      cACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAg                                                                          <  1:3737118/150‑1 (MQ=255)
                                                                       aCCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGc                                                                         >  2:86477/1‑150 (MQ=255)
                                                                         cGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaa                                                                       <  1:2946151/150‑1 (MQ=255)
                                                                         cGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaa                                                                       >  2:4129309/1‑150 (MQ=255)
                                                                         cGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaa                                                                       >  1:661158/1‑150 (MQ=255)
                                                                          gCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaaa                                                                      <  1:4338687/150‑1 (MQ=255)
                                                                           cAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaaaa                                                                     >  2:4565664/1‑150 (MQ=255)
                                                                           cAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaaaa                                                                     >  2:168902/1‑150 (MQ=255)
                                                                               cacTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCa                                                                 <  2:1230387/150‑1 (MQ=255)
                                                                                 cTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGc                                                               <  1:4024195/150‑1 (MQ=255)
                                                                                    cccGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACg                                                            >  1:2912874/1‑150 (MQ=255)
                                                                                      cGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGcc                                                          >  2:641952/1‑150 (MQ=255)
                                                                                       gTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGccc                                                         >  2:265057/1‑150 (MQ=255)
                                                                                       gTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGccc                                                         <  2:370902/150‑1 (MQ=255)
                                                                                       gTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGccc                                                         <  2:1132880/150‑1 (MQ=255)
                                                                                       gTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGccc                                                         >  1:924016/1‑150 (MQ=255)
                                                                                             gaagaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCGGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTgg                                                   <  2:3040351/150‑1 (MQ=255)
                                                                                              aagaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGa                                                  >  2:2945482/1‑150 (MQ=255)
                                                                                                gaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATa                                                >  2:2876510/1‑150 (MQ=255)
                                                                                                  aTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAg                                              <  1:1257696/150‑1 (MQ=255)
                                                                                                     tttGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCa                                           <  2:2837434/150‑1 (MQ=255)
                                                                                                       tGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAgg                                         <  1:13438/150‑1 (MQ=255)
                                                                                                        gCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGgc                                        >  2:4447461/1‑150 (MQ=255)
                                                                                                         ccTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGgcg                                       >  2:239381/1‑150 (MQ=255)
                                                                                                             ttCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAg                                   <  1:1914955/150‑1 (MQ=255)
                                                                                                                    gcAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCg                            >  1:1412124/1‑150 (MQ=255)
                                                                                                                      aaCATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTg                          >  2:2401291/1‑150 (MQ=255)
                                                                                                                      aaCATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTg                          >  1:2300557/1‑150 (MQ=255)
                                                                                                                         aTCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTg                       <  2:4108809/150‑1 (MQ=255)
                                                                                                                          tCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTgg                      <  2:3186809/150‑1 (MQ=255)
                                                                                                                               cATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTc                 <  2:3663308/150‑1 (MQ=255)
                                                                                                                                aTTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCa                <  2:3832168/150‑1 (MQ=255)
                                                                                                                                 ttACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCaa               >  2:4305620/1‑150 (MQ=255)
                                                                                                                                     tGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGcc           >  2:2695418/1‑150 (MQ=255)
                                                                                                                                       gATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTc         <  1:2397426/150‑1 (MQ=255)
                                                                                                                                       gATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTc         <  2:4675557/150‑1 (MQ=255)
                                                                                                                                       gATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTc         <  1:14866/150‑1 (MQ=255)
                                                                                                                                        aTACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCa        <  2:2334660/150‑1 (MQ=255)
                                                                                                                                          aCGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCAcc      <  1:3802925/150‑1 (MQ=255)
                                                                                                                                              tCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTaac  <  2:780067/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AACTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAAC  >  CP009273/104414‑104705

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: