Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,366,254 C→T 15.1% D228N (GAT→AAT)  nanA ← N‑acetylneuraminate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,366,2540CT15.1% 814.0 / 114.7 292D228N (GAT→AAT) nanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base C (116/131);  new base T (27/17);  total (143/149)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 305 total aligned reads displayed.

TACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGAC  >  CP009273/3366105‑3366403
                                                                                                                                                     |                                                                                                                                                     
  ccGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGc                                                                                                                                                     >  1:448919/1‑150 (MQ=255)
    ggTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCct                                                                                                                                                   <  2:2117190/150‑1 (MQ=255)
    ggTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCct                                                                                                                                                   >  2:3400983/1‑150 (MQ=255)
    ggTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCct                                                                                                                                                   <  2:636825/150‑1 (MQ=255)
         aaaCGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCttt                                                                                                                                              >  1:541665/1‑150 (MQ=255)
                  gcggcacagcggcacaGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCtt                                                                                                                                     <  2:2753019/150‑1 (MQ=255)
                   cggcacagcggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTa                                                                                                                                    >  1:4676489/1‑150 (MQ=255)
                        cagcggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGAt                                                                                                                               <  2:1446091/150‑1 (MQ=255)
                        cagcggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGAt                                                                                                                               >  2:592160/1‑150 (MQ=255)
                        cagcggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGAt                                                                                                                               >  2:4535804/1‑150 (MQ=255)
                         agcggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATc                                                                                                                              >  1:3673915/1‑150 (MQ=255)
                           cggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATccc                                                                                                                            <  2:3673305/150‑1 (MQ=255)
                            ggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATcccc                                                                                                                           >  2:1976821/1‑150 (MQ=255)
                            ggcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATcccc                                                                                                                           >  1:140997/1‑150 (MQ=255)
                             gcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCt                                                                                                                          <  2:1063125/150‑1 (MQ=255)
                             gcacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCt                                                                                                                          <  1:3363047/150‑1 (MQ=255)
                              cacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTg                                                                                                                         <  2:1251367/150‑1 (MQ=255)
                              cacagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTg                                                                                                                         <  2:2483846/150‑1 (MQ=255)
                                cagAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGAt                                                                                                                       <  1:4118684/150‑1 (MQ=255)
                                 agAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATa                                                                                                                      >  1:3835418/1‑150 (MQ=255)
                                      cGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGcc                                                                                                                 <  1:2056937/150‑1 (MQ=255)
                                       gACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCa                                                                                                                <  1:509679/150‑1 (MQ=255)
                                          aTCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGcc                                                                                                             >  1:3409749/1‑150 (MQ=255)
                                          aTCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGcc                                                                                                             >  1:1591719/1‑150 (MQ=255)
                                            ccATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCa                                                                                                           <  2:3821320/150‑1 (MQ=255)
                                             cATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCat                                                                                                          >  1:3243979/1‑150 (MQ=255)
                                              atatAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCatg                                                                                                         >  2:674203/1‑150 (MQ=255)
                                               tataATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCatga                                                                                                        >  1:871441/1‑150 (MQ=255)
                                               tataATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCatga                                                                                                        <  2:4798248/150‑1 (MQ=255)
                                                 taATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCatgatg                                                                                                      >  2:885749/1‑150 (MQ=255)
                                                 taATGGAGGACAGTTTTCAGGCCGCGGAATACGCCAGTTTTGATCATTAAATCAATGACTTTATTGCATTCAGTTTGCATTTTATGCTCGGTCTGGATATTGACTTCTTTCAGCGCCTTAACGATCACCTGATAGCGCCAGCCCatgatg                                                                                                      <  2:2274159/150‑1 (MQ=255)
                                                   aTGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGAtgtt                                                                                                    >  1:2399948/1‑150 (MQ=255)
                                                      gaggaCAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTa                                                                                                 <  1:3096585/150‑1 (MQ=255)
                                                          aCAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTa                                                                                             <  2:1222137/150‑1 (MQ=255)
                                                               tttCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGcc                                                                                        >  2:2649944/1‑150 (MQ=255)
                                                                  cAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGAt                                                                                     >  2:756437/1‑150 (MQ=255)
                                                                   aGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATa                                                                                    <  1:958832/150‑1 (MQ=255)
                                                                   aGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTAATGCCGATa                                                                                    >  1:1640724/1‑150 (MQ=255)
                                                                    ggCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATac                                                                                   <  1:1594639/150‑1 (MQ=255)
                                                                        gcgGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATaccacc                                                                               <  1:947260/150‑1 (MQ=255)
                                                                         cgGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATAccacca                                                                              >  2:4604031/1‑150 (MQ=255)
                                                                         cgGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATAccacca                                                                              >  2:3001645/1‑150 (MQ=255)
                                                                          ggAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCAt                                                                             >  2:4226221/1‑150 (MQ=255)
                                                                           gAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATc                                                                            >  1:121673/1‑150 (MQ=255)
                                                                             aTACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAg                                                                          <  1:3639367/150‑1 (MQ=255)
                                                                             aTACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAg                                                                          <  2:2080499/150‑1 (MQ=255)
                                                                               aCGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAgcg                                                                        <  1:2482500/150‑1 (MQ=255)
                                                                                 gCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGcgcc                                                                      >  2:3377587/1‑150 (MQ=255)
                                                                                 gCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGcgcc                                                                      <  2:2554627/150‑1 (MQ=255)
                                                                                  cccGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGcgccc                                                                     <  2:4072755/150‑1 (MQ=255)
                                                                                   ccGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGcgcccg                                                                    >  1:747285/1‑150 (MQ=255)
                                                                                   ccGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGcgcccg                                                                    >  1:750077/1‑150 (MQ=255)
                                                                                      ttttGATCAGTAAATCACTGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTTGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCca                                                                 <  1:3446866/150‑1 (MQ=255)
                                                                                      ttttGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCca                                                                 <  2:541645/150‑1 (MQ=255)
                                                                                       tttGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCcag                                                                <  2:2355142/150‑1 (MQ=255)
                                                                                       tttGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCcag                                                                <  1:2113421/150‑1 (MQ=255)
                                                                                          gATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCcagcag                                                             <  1:3552359/150‑1 (MQ=255)
                                                                                            tCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGAc                                                           >  2:4138614/1‑150 (MQ=255)
                                                                                                tAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCaga                                                       >  2:808916/1‑150 (MQ=255)
                                                                                                 aaaTCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCagag                                                      >  2:1666810/1‑150 (MQ=255)
                                                                                                  aaTCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAgg                                                     >  2:2173828/1‑150 (MQ=255)
                                                                                                  aaTCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAgg                                                     >  1:1311194/1‑150 (MQ=255)
                                                                                                     cAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGa                                                  >  1:1276519/1‑150 (MQ=255)
                                                                                                      aaTGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGaa                                                 <  1:1747253/150‑1 (MQ=255)
                                                                                                       aTGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAg                                                >  2:917073/1‑150 (MQ=255)
                                                                                                           cTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGAttt                                            <  2:3233209/150‑1 (MQ=255)
                                                                                                            tttATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTtc                                           >  1:4633383/1‑150 (MQ=255)
                                                                                                               aTTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTtcgtc                                        <  1:2684037/150‑1 (MQ=255)
                                                                                                               aTTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTtcgtc                                        >  1:3982896/1‑150 (MQ=255)
                                                                                                               aTTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTtcgtc                                        <  1:1948244/150‑1 (MQ=255)
                                                                                                                ttGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTtcgtcg                                       <  2:4756649/150‑1 (MQ=255)
                                                                                                                  gCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTa                                     >  1:3870572/1‑150 (MQ=255)
                                                                                                                    aTTCAGTTTGCAGTTTCTGCGCGGTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAAc                                   >  2:4115898/1‑150 (MQ=255)
                                                                                                                     ttCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAAcc                                  <  1:4172305/150‑1 (MQ=255)
                                                                                                                     ttCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAAcc                                  >  2:3108104/1‑150 (MQ=255)
                                                                                                                       cAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGt                                <  1:1043175/150‑1 (MQ=255)
                                                                                                                       cAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGt                                >  1:3573480/1‑150 (MQ=255)
                                                                                                                         gTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTta                              <  1:2840304/150‑1 (MQ=255)
                                                                                                                         gTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTta                              >  2:1900944/1‑150 (MQ=255)
                                                                                                                             gCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATaga                          <  1:376011/150‑1 (MQ=255)
                                                                                                                                gTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGca                       >  2:4326392/1‑150 (MQ=255)
                                                                                                                                         cgGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAg              <  2:2825789/150‑1 (MQ=255)
                                                                                                                                          ggTCTGGATATTGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAgg             >  2:16155/1‑150 (MQ=255)
                                                                                                                                          ggTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAgg             >  2:1387456/1‑150 (MQ=255)
                                                                                                                                           gTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATAAGGa            <  1:698031/150‑1 (MQ=255)
                                                                                                                                           gTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATAAGGa            <  1:3130291/150‑1 (MQ=255)
                                                                                                                                              tGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGt         >  1:4246006/1‑150 (MQ=255)
                                                                                                                                              tGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGt         >  2:2217898/1‑150 (MQ=255)
                                                                                                                                               ggATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGtt        <  2:1490854/150‑1 (MQ=255)
                                                                                                                                               ggATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGtt        <  1:1845059/150‑1 (MQ=255)
                                                                                                                                                 atctCGCCTTCCGTCAGCGCCTTAACGATCCCCTCAGAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCTCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCa      <  2:4694947/147‑1 (MQ=255)
                                                                                                                                                 atatCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCa      <  2:1817599/150‑1 (MQ=255)
                                                                                                                                                  tatCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCac     <  2:444567/150‑1 (MQ=255)
                                                                                                                                                  tatCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCac     <  2:3041449/150‑1 (MQ=255)
                                                                                                                                                  tatCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCac     <  2:2877163/150‑1 (MQ=255)
                                                                                                                                                   atCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCacg    <  1:2788836/150‑1 (MQ=255)
                                                                                                                                                   atCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCacg    <  1:2788983/150‑1 (MQ=255)
                                                                                                                                                     cGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCacgac  <  1:2175691/150‑1 (MQ=255)
                                                                                                                                                     cGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCacgac  <  2:3572866/150‑1 (MQ=255)
                                                                                                                                                     cGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCacgac  >  2:1140057/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACCGGTCCAAACGGTTTGCGGCACAGCGGCACAGAAACGACATCCATATAATGGAGGACAGTTTTCAGGCCGCGGAATACGCCCGTTTTGATCAGTAAATCAATGACTTTATTGCATTCAGTTTGCAGTTTCTGCGCGGTCTGGATATCGCCTTCTTTCAGCGCCTTAACGATCCCCTGATAGCGCCAGCCCATGATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGAC  >  CP009273/3366105‑3366403

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: