Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 104,561 G→T 100% 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,5610GT98.3% 869.9 / 3.0 250A124S (GCG→TCG) secMregulator of secA translation
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (120/125);  total (124/126)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 252 total aligned reads displayed.

AACTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAA  >  CP009273/104414‑104707
                                                                                                                                                   |                                                                                                                                                  
aaCTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAgctcg                                                                                                                                                  <  1:4040586/150‑1 (MQ=255)
aaCTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAgctcg                                                                                                                                                  >  1:1406544/1‑150 (MQ=255)
  cTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCgct                                                                                                                                                <  2:1019066/150‑1 (MQ=255)
  cTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCgct                                                                                                                                                <  2:1586753/150‑1 (MQ=255)
    aTTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCTTTACTGGATACGCTCAGCTCgctgc                                                                                                                                              <  2:2110325/150‑1 (MQ=255)
      tCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGctgctg                                                                                                                                            <  2:1847720/150‑1 (MQ=255)
      tCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGctgctg                                                                                                                                            <  2:1847591/150‑1 (MQ=255)
         gTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGAcc                                                                                                                                         >  1:2172775/1‑150 (MQ=255)
              ttACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGa                                                                                                                                    <  1:608801/150‑1 (MQ=255)
                aCTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAg                                                                                                                                  >  1:3742256/1‑150 (MQ=255)
                    gCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGctc                                                                                                                              <  1:1319733/150‑3 (MQ=255)
                    gCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGctc                                                                                                                              <  2:1639968/150‑3 (MQ=255)
                    gCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCAc                                                                                                                              <  2:3152541/150‑1 (MQ=255)
                      aTCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGc                                                                                                                            >  2:604599/1‑150 (MQ=255)
                       tCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGcc                                                                                                                           <  2:2553157/150‑1 (MQ=255)
                       tCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGcc                                                                                                                           <  1:536437/150‑1 (MQ=255)
                         aaCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGt                                                                                                                         <  1:3759882/150‑1 (MQ=255)
                         aaCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGt                                                                                                                         <  1:842927/150‑1 (MQ=255)
                          aCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTc                                                                                                                        <  1:2218990/150‑1 (MQ=255)
                          aCATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTc                                                                                                                        <  1:2001859/150‑1 (MQ=255)
                               ccATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaaa                                                                                                                   <  1:1449647/150‑1 (MQ=255)
                               ccATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGaaa                                                                                                                   >  2:3767887/1‑150 (MQ=255)
                                 aTTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAg                                                                                                                 >  1:1811192/1‑150 (MQ=255)
                                 aTTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAg                                                                                                                 >  1:2542806/1‑150 (MQ=255)
                                  ttCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAgg                                                                                                                <  2:2487910/150‑1 (MQ=255)
                                   tCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAggg                                                                                                               <  1:4137167/150‑1 (MQ=255)
                                    cGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGt                                                                                                              <  2:1760263/150‑1 (MQ=255)
                                     gCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGtt                                                                                                             <  2:3456730/150‑1 (MQ=255)
                                          gTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTTCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCAGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGc                                                                                                        >  1:54433/1‑150 (MQ=255)
                                           tAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCa                                                                                                       <  2:1426039/150‑1 (MQ=255)
                                            aaTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGNTGCTGACCCNGGAAGGCACGCCGTCTGAAAAGGGTTNTCGCAt                                                                                                      >  1:2847281/1‑150 (MQ=255)
                                            aaTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCAt                                                                                                      >  2:2629910/1‑150 (MQ=255)
                                             aTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAtt                                                                                                     <  2:880408/150‑1 (MQ=255)
                                             aTCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCAtt                                                                                                     <  1:1116035/150‑1 (MQ=255)
                                              tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                    >  1:2876144/1‑150 (MQ=255)
                                              tCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTg                                                                                                    >  2:3863601/1‑150 (MQ=255)
                                                 gTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGAtt                                                                                                 <  2:807482/150‑1 (MQ=255)
                                                 gTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGACTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGAtt                                                                                                 >  1:1148209/1‑150 (MQ=255)
                                                   cATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTAt                                                                                               <  1:1145217/150‑1 (MQ=255)
                                                   cATCTTACTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTAt                                                                                               >  1:1764511/1‑150 (MQ=255)
                                                    ctctttcttTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATg                                                                                              <  2:671411/149‑1 (MQ=255)
                                                     tctttcttTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATgc                                                                                             >  2:30367/1‑150 (MQ=255)
                                                      ctttctttCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTAGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATgcg                                                                                            >  2:142001/1‑150 (MQ=255)
                                                       tttctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATgcgc                                                                                           >  1:3790604/1‑150 (MQ=255)
                                                        ttctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCa                                                                                          >  2:1256184/1‑150 (MQ=255)
                                                         tctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAt                                                                                         <  2:3888303/150‑1 (MQ=255)
                                                         tctttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAt                                                                                         <  2:400785/150‑1 (MQ=255)
                                                           tttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAttt                                                                                       >  1:1279353/1‑150 (MQ=255)
                                                           tttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAttt                                                                                       <  2:1140607/150‑1 (MQ=255)
                                                            ttcGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCAtttt                                                                                      <  2:3143785/150‑1 (MQ=255)
                                                              cGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAc                                                                                    >  1:3827367/1‑150 (MQ=255)
                                                              cGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAc                                                                                    >  1:3251610/1‑150 (MQ=255)
                                                              cGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAc                                                                                    <  2:1097697/150‑1 (MQ=255)
                                                               gCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAcc                                                                                   >  2:3299480/1‑150 (MQ=255)
                                                               gCAATGGCACCGCAAACACGGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAcc                                                                                   <  1:2909413/150‑1 (MQ=255)
                                                                cAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTAccc                                                                                  >  1:1912773/1‑150 (MQ=255)
                                                                    ggCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCcaca                                                                              <  1:1345898/150‑1 (MQ=255)
                                                                    ggCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCcaca                                                                              >  2:2345827/1‑150 (MQ=255)
                                                                     gCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACaa                                                                             >  1:1888150/1‑150 (MQ=255)
                                                                       aCCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGc                                                                           >  2:131426/1‑150 (MQ=255)
                                                                           cAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCaaaa                                                                       >  1:3720247/1‑150 (MQ=255)
                                                                            aaaCACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACCTCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGAGTATGCGCATTTTACCCCACAAGCAAAAt                                                                      >  2:1056715/1‑150 (MQ=255)
                                                                               cacTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCa                                                                   >  1:3426975/1‑150 (MQ=255)
                                                                                acTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAg                                                                  <  1:2576870/150‑1 (MQ=255)
                                                                                   gCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACAGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCTCGCCGTCAGAAAAGGGTTATCGCAGTGATTATGCGCATTTTACACCACAAGCAAAATTCAGCAc                                                               >  1:333660/1‑150 (MQ=255)
                                                                                    cccGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACg                                                              <  2:2750708/150‑1 (MQ=255)
                                                                                       gTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGccc                                                           <  1:2693257/150‑1 (MQ=255)
                                                                                       gTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGccc                                                           >  1:2202080/1‑150 (MQ=255)
                                                                                         tGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGt                                                         >  1:1412745/1‑150 (MQ=255)
                                                                                           cTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCt                                                       >  2:432945/1‑150 (MQ=255)
                                                                                           cTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCt                                                       >  1:1584342/1‑150 (MQ=255)
                                                                                            tGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTg                                                      <  2:2801784/150‑1 (MQ=255)
                                                                                               agaaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGAt                                                   <  1:76789/150‑1 (MQ=255)
                                                                                                 aaTCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATaa                                                 >  1:1465606/1‑150 (MQ=255)
                                                                                                   tCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGc                                               <  2:4096134/150‑1 (MQ=255)
                                                                                                   tCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGc                                               >  2:1599879/1‑150 (MQ=255)
                                                                                                   tCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGc                                               <  2:554813/150‑1 (MQ=255)
                                                                                                        gCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGgc                                          <  2:2298788/150‑1 (MQ=255)
                                                                                                         ccTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGgcg                                         >  1:968501/1‑150 (MQ=255)
                                                                                                           tctTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCa                                       <  2:3244155/150‑1 (MQ=255)
                                                                                                              tCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAgg                                    <  1:2327109/150‑1 (MQ=255)
                                                                                                                aGGCGCAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCa                                  >  1:1177329/1‑150 (MQ=255)
                                                                                                                   cgcAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATcc                               >  2:3618841/1‑150 (MQ=255)
                                                                                                                   cgcAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATcc                               <  2:2944334/150‑1 (MQ=255)
                                                                                                                   cgcAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATcc                               <  1:2289239/150‑1 (MQ=255)
                                                                                                                    gcAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCg                              >  1:2609134/1‑150 (MQ=255)
                                                                                                                     cAACATCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGt                             >  2:3686198/1‑150 (MQ=255)
                                                                                                                          tCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTgg                        >  1:2441095/1‑150 (MQ=255)
                                                                                                                          tCTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTgg                        >  1:2441171/1‑150 (MQ=255)
                                                                                                                           cTTGCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGc                       <  1:2005641/150‑1 (MQ=255)
                                                                                                                              gCATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCt                    >  1:2639888/1‑150 (MQ=255)
                                                                                                                               cATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTc                   >  2:1732693/1‑150 (MQ=255)
                                                                                                                               cATTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTc                   <  1:2759415/150‑1 (MQ=255)
                                                                                                                                aTTACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCa                  <  1:1665945/150‑1 (MQ=255)
                                                                                                                                  tACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAAc                <  1:3615268/150‑1 (MQ=255)
                                                                                                                                  tACTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAAc                >  2:509103/1‑150 (MQ=255)
                                                                                                                                   aCTGGATACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACg               >  1:3092108/1‑150 (MQ=255)
                                                                                                                                        aTACGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCa          >  2:2105510/1‑150 (MQ=255)
                                                                                                                                          aCGCTCAGCTCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCAcc        <  1:2420188/150‑1 (MQ=255)
                                                                                                                                           cGCTCAGCTCGCTGCTGACCCAGGAAGGCTCGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCt       <  1:1022422/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACTATTCCGTTGATTACTGGCATCAACATGCCATTCGCACGGTAATCCGTCATCTTTCTTTCGCAATGGCACCGCAAACACTGCCCGTTGCTGAAGAATCTTTGCCTCTTCAGGCGCAACATCTTGCATTACTGGATACGCTCAGCGCGCTGCTGACCCAGGAAGGCACGCCGTCTGAAAAGGGTTATCGCATTGATTATGCGCATTTTACCCCACAAGCAAAATTCAGCACGCCCGTCTGGATAAGCCAGGCGCAAGGCATCCGTGCTGGCCCTCAACGCCTCACCTAACAA  >  CP009273/104414‑104707

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: