Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,048,928 C→A 100% E903* (GAA→TAA)  torS ← hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,048,9280CA99.6% 905.6 / ‑4.4 269E903* (GAA→TAA) torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (156/112);  total (156/113)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 269 total aligned reads displayed.

GCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAT  >  CP009273/1048782‑1049072
                                                                                                                                                  |                                                                                                                                                
 ccgccgCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACCATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACag                                                                                                                                               <  1:150317/149‑1 (MQ=255)
   gccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTTATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACagag                                                                                                                                             >  1:3463104/1‑149 (MQ=255)
   gccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACagag                                                                                                                                             >  2:3955767/1‑149 (MQ=255)
    ccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGc                                                                                                                                            >  2:1620117/1‑149 (MQ=255)
     cgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCg                                                                                                                                           >  2:2253530/1‑149 (MQ=255)
       cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACATAGCGGc                                                                                                                                         >  2:2540539/1‑149 (MQ=255)
          ccAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAac                                                                                                                                      <  1:1019597/149‑1 (MQ=255)
          ccAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAac                                                                                                                                      <  2:2577711/149‑1 (MQ=255)
           cAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAaca                                                                                                                                     >  1:3355901/1‑149 (MQ=255)
            aGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAacac                                                                                                                                    <  2:2099250/149‑1 (MQ=255)
               gggTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAAAACTGc                                                                                                                                 >  2:3679614/1‑149 (MQ=255)
                 gTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgt                                                                                                                               <  1:3497246/149‑1 (MQ=255)
                 gTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgt                                                                                                                               >  2:328452/1‑149 (MQ=255)
                 gTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgt                                                                                                                               >  2:1174126/1‑149 (MQ=255)
                     ttttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTaa                                                                                                                           <  2:2613990/149‑1 (MQ=255)
                     ttttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTaa                                                                                                                           >  2:911243/1‑149 (MQ=255)
                       ttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAAtt                                                                                                                         <  1:4024673/149‑1 (MQ=255)
                         aTTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttc                                                                                                                       <  1:3839401/149‑1 (MQ=255)
                         aTTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttc                                                                                                                       >  1:1693260/1‑149 (MQ=255)
                            ctctGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttcttc                                                                                                                    <  2:620320/149‑1 (MQ=255)
                            ctcgGCCATGATTTTGACATTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGTGCGGCAACTCTGCGTGTAATttcttc                                                                                                                    <  2:3024837/149‑1 (MQ=255)
                              ctGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCgt                                                                                                                  <  2:897495/149‑1 (MQ=255)
                               tGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCgtg                                                                                                                 >  1:4092624/1‑149 (MQ=255)
                                    tGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                            >  2:4031672/1‑149 (MQ=255)
                                    tGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                            <  1:1958454/149‑1 (MQ=255)
                                     gATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAg                                                                                                           >  1:755092/1‑149 (MQ=255)
                                     gATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAg                                                                                                           >  1:958120/1‑149 (MQ=255)
                                       ttttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAggg                                                                                                         <  1:2783846/149‑1 (MQ=255)
                                         ttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAggggg                                                                                                       >  1:3085529/1‑149 (MQ=255)
                                         ttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAggggg                                                                                                       <  2:3497697/149‑1 (MQ=255)
                                           gCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGc                                                                                                     >  1:600602/1‑149 (MQ=255)
                                           gCCACTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGc                                                                                                     <  1:2395725/149‑1 (MQ=255)
                                               gTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGAt                                                                                                 >  1:1130244/1‑149 (MQ=255)
                                               gTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGAt                                                                                                 <  2:3385151/149‑1 (MQ=255)
                                                tCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAACTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATa                                                                                                <  2:1073541/149‑1 (MQ=255)
                                                    ttAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATgg                                                                                            >  2:272119/1‑149 (MQ=255)
                                                        gTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgc                                                                                        >  2:3812130/1‑149 (MQ=255)
                                                         ttGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgct                                                                                       >  1:2231244/1‑149 (MQ=255)
                                                         ttGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgct                                                                                       <  2:3296359/149‑1 (MQ=255)
                                                          tGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgctg                                                                                      <  2:2201961/149‑1 (MQ=255)
                                                            tataCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgctgct                                                                                    >  2:2448594/1‑149 (MQ=255)
                                                              taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTcc                                                                                  <  2:3015738/149‑1 (MQ=255)
                                                              taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTcc                                                                                  >  1:327385/1‑149 (MQ=255)
                                                                    aTGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTtgt                                                                            >  1:3738032/1‑149 (MQ=255)
                                                                    aTGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTtgt                                                                            >  1:130569/1‑149 (MQ=255)
                                                                    aTGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTAAAATCGCTTTCAGGTCTTTCTTATGCAACAATGCTTACAGAGCGGCAACACTGCGTTTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGGtgg                                                                            >  2:3677882/1‑148 (MQ=255)
                                                                       ccACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                                         >  1:477820/1‑149 (MQ=255)
                                                                         aCTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                                       >  2:547636/1‑149 (MQ=255)
                                                                          cTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                                      >  2:2179400/1‑149 (MQ=255)
                                                                          cTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                                      >  2:3208128/1‑149 (MQ=255)
                                                                                tAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctgg                                                                >  2:1707984/1‑149 (MQ=255)
                                                                                 aGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                               <  2:3653692/149‑1 (MQ=255)
                                                                                 aGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                               >  1:3221956/1‑149 (MQ=255)
                                                                                  gttgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                              >  2:1549084/1‑149 (MQ=255)
                                                                                  gttgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                              >  1:3167298/1‑149 (MQ=255)
                                                                                   ttgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggctg                                                             >  1:1418449/1‑149 (MQ=255)
                                                                                    tgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                            >  2:2879055/1‑149 (MQ=255)
                                                                                       tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                                         <  1:394276/149‑1 (MQ=255)
                                                                                       tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                                         <  1:2214516/149‑1 (MQ=255)
                                                                                        cTTAACTGGCTGATATTGATTCAACTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                                        <  2:2976464/149‑1 (MQ=255)
                                                                                        cTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                                        >  1:1039901/1‑149 (MQ=255)
                                                                                          tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                                      <  1:1174126/149‑1 (MQ=255)
                                                                                            aCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGc                                                    >  2:1435427/1‑149 (MQ=255)
                                                                                             cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCa                                                   >  2:1108486/1‑149 (MQ=255)
                                                                                              tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAt                                                  <  2:3355901/149‑1 (MQ=255)
                                                                                               ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                                 <  1:3236973/149‑1 (MQ=255)
                                                                                               ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                                 >  2:2497565/1‑149 (MQ=255)
                                                                                                 cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                               <  1:1877589/149‑1 (MQ=255)
                                                                                                   gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                             >  2:73388/1‑149 (MQ=255)
                                                                                                    atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                            >  1:3330779/1‑149 (MQ=255)
                                                                                                    atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                            >  2:2869903/1‑149 (MQ=255)
                                                                                                     tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGa                                           >  1:2834832/1‑149 (MQ=255)
                                                                                                      atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGAc                                          <  1:803179/149‑1 (MQ=255)
                                                                                                           ttCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                                     <  2:1871086/149‑1 (MQ=255)
                                                                                                              aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCaa                                  >  1:3768543/1‑149 (MQ=255)
                                                                                                                 tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                               >  1:1677052/1‑149 (MQ=255)
                                                                                                                 tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                               >  2:1754997/1‑149 (MQ=255)
                                                                                                                 tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                               >  1:2719341/1‑149 (MQ=255)
                                                                                                                 tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                               >  2:3795111/1‑149 (MQ=255)
                                                                                                                      ttCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                          >  2:1919025/1‑149 (MQ=255)
                                                                                                                      ttCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                          <  1:1622591/149‑1 (MQ=255)
                                                                                                                        cAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTa                        <  1:1414509/149‑1 (MQ=255)
                                                                                                                         aGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCACCTGCTTTTTAg                       >  2:3920368/1‑149 (MQ=255)
                                                                                                                         aGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCATTTTAg                       <  2:3598664/149‑1 (MQ=255)
                                                                                                                          ggTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc                      >  1:1807484/1‑149 (MQ=255)
                                                                                                                          ggTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc                      >  1:1300003/1‑149 (MQ=255)
                                                                                                                           gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt                     >  1:924874/1‑149 (MQ=255)
                                                                                                                           gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt                     >  1:3559029/1‑149 (MQ=255)
                                                                                                                           gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt                     <  2:1306485/149‑1 (MQ=255)
                                                                                                                             ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGa                   <  1:3119147/149‑1 (MQ=255)
                                                                                                                               ttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg                 <  1:566168/149‑1 (MQ=255)
                                                                                                                               ttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg                 >  2:2940419/1‑149 (MQ=255)
                                                                                                                               ttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg                 >  2:358396/1‑149 (MQ=255)
                                                                                                                                  ttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcgc              <  2:1989537/149‑1 (MQ=255)
                                                                                                                                  ttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcgc              >  1:2156011/1‑149 (MQ=255)
                                                                                                                                     tGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGc           >  1:2936092/1‑149 (MQ=255)
                                                                                                                                       cAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCAc         <  1:3186683/149‑1 (MQ=255)
                                                                                                                                         aCCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGc       >  1:264240/1‑149 (MQ=255)
                                                                                                                                         aCCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGc       >  1:2460933/1‑149 (MQ=255)
                                                                                                                                          ccATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGAGGCACGCt      >  1:3435600/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAT  >  CP009273/1048782‑1049072

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: