Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,917,059 A→T 23.1% D56E (GAT→GAA pphA ← serine/threonine‑specific protein phosphatase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,917,0590AT23.1% 468.4 / 145.6 207D56E (GAT→GAApphAserine/threonine‑specific protein phosphatase 1
Reads supporting (aligned to +/- strand):  ref base A (80/79);  new base T (13/35);  total (93/114)
Fisher's exact test for biased strand distribution p-value = 4.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 230 total aligned reads displayed.

CAATAAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAAC  >  CP009273/1916911‑1917208
                                                                                                                                                    |                                                                                                                                                     
cAATAAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCAcgat                                                                                                                                                      <  2:4431883/150‑1 (MQ=255)
cAATAAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCAcgat                                                                                                                                                      <  2:2529227/150‑1 (MQ=255)
cAATAAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCAcgat                                                                                                                                                      <  2:2529465/150‑1 (MQ=255)
 aaTTAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCATTGTTGTTCCAGTAACTGCAGACAACGTTAACTTTGCGGCCCAcgatc                                                                                                                                                     <  1:2695247/150‑1 (MQ=255)
 aaTAAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCAcgatc                                                                                                                                                     >  1:1858423/1‑150 (MQ=255)
 aaTAAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCAcgatc                                                                                                                                                     <  2:3170479/150‑1 (MQ=255)
    aaaCCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCAcgatcgat                                                                                                                                                  >  2:1320110/1‑150 (MQ=255)
    aaaCCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCAcgatcgat                                                                                                                                                  >  2:1320758/1‑150 (MQ=255)
     aaCCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATa                                                                                                                                                 <  2:3814785/150‑1 (MQ=255)
         aGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGt                                                                                                                                             >  1:765957/1‑150 (MQ=255)
         aGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGt                                                                                                                                             >  1:4416863/1‑150 (MQ=255)
         aGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGt                                                                                                                                             >  1:150451/1‑150 (MQ=255)
              ccgccATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTAGCGGCCCACGATCGATAACGTCTccc                                                                                                                                        >  2:1620969/1‑150 (MQ=255)
                 ccATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACt                                                                                                                                     <  1:1975272/150‑1 (MQ=255)
                  cATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTg                                                                                                                                    <  2:4538099/150‑1 (MQ=255)
                    ttCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTgag                                                                                                                                  >  1:1588998/1‑150 (MQ=255)
                       atcaACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATa                                                                                                                               <  1:991298/150‑1 (MQ=255)
                          aaCCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGt                                                                                                                            >  1:1571839/1‑150 (MQ=255)
                          aaCCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGt                                                                                                                            <  2:2974775/150‑1 (MQ=255)
                          aaCCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGt                                                                                                                            <  2:542544/150‑1 (MQ=255)
                          aaCCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGt                                                                                                                            >  2:3564763/1‑150 (MQ=255)
                            ccACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTaa                                                                                                                          >  2:122654/1‑150 (MQ=255)
                             cacaAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTaaa                                                                                                                         >  1:317603/1‑150 (MQ=255)
                             cacaAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTaaa                                                                                                                         >  1:4021013/1‑150 (MQ=255)
                              acaAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAAt                                                                                                                        <  1:4727160/150‑1 (MQ=255)
                                aaaGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAAtct                                                                                                                      <  1:2121398/150‑1 (MQ=255)
                                aaaGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAAtct                                                                                                                      <  2:1862213/150‑1 (MQ=255)
                                aaaGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAAtct                                                                                                                      <  1:830674/150‑1 (MQ=255)
                                aaaGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCAAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAAtct                                                                                                                      >  2:1768700/1‑150 (MQ=255)
                                 aaGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAAtctc                                                                                                                     <  2:1620078/150‑1 (MQ=255)
                                 aaGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAAtctc                                                                                                                     <  2:3275695/150‑1 (MQ=255)
                                  aGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCg                                                                                                                    >  2:2016276/1‑150 (MQ=255)
                                     cATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCa                                                                                                                 <  2:4274498/150‑1 (MQ=255)
                                      aTCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGAGTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCAc                                                                                                                <  2:2646191/150‑1 (MQ=255)
                                         tgctgGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGa                                                                                                             <  2:2837922/150‑1 (MQ=255)
                                          gctgGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGAt                                                                                                            <  2:3228872/150‑1 (MQ=255)
                                           ctgGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATc                                                                                                           <  2:3859562/150‑1 (MQ=255)
                                           ctgGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATc                                                                                                           <  2:1396529/150‑1 (MQ=255)
                                           ctgGGATGCCAGCGCATCCATCGCCATCTGTTCAGGATTGCCTCGTACCGCACAAACCCAAGGGTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGAGAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATc                                                                                                           <  2:3954833/150‑1 (MQ=255)
                                             gggATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCaa                                                                                                         <  1:2373801/150‑1 (MQ=255)
                                             gggATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCaa                                                                                                         <  1:2373704/150‑1 (MQ=255)
                                              ggATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCaaa                                                                                                        <  1:4374146/150‑1 (MQ=255)
                                                aTGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAAc                                                                                                      >  1:2545175/1‑150 (MQ=255)
                                                aTGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAAc                                                                                                      <  1:79619/150‑1 (MQ=255)
                                                   ccAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGAc                                                                                                   <  1:191018/150‑1 (MQ=255)
                                                   ccAGCGCATCCATCGCCATCTGTTCACGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGAc                                                                                                   >  2:1988218/1‑150 (MQ=255)
                                                     aGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACaa                                                                                                 >  1:116238/1‑150 (MQ=255)
                                                          atccatcGCAATGTGTTCATGTTGGCCTCTTACCGCACAAACCCAATGTTGTTCCATTAACTGATTACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGTTAATTAAATCTCGCCACGGATAAAAACGACAATGCCa                                                                                            <  2:916651/150‑1 (MQ=255)
                                                            ccatcGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATa                                                                                          >  2:4749874/1‑150 (MQ=255)
                                                             catcGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATaa                                                                                         >  2:3727297/1‑150 (MQ=255)
                                                                  ccATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTgc                                                                                    >  2:1553340/1‑150 (MQ=255)
                                                                  ccATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTgc                                                                                    >  1:3929243/1‑150 (MQ=255)
                                                                    aTCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTgcgg                                                                                  >  1:4489414/1‑150 (MQ=255)
                                                                     tCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACCATGCCATAATTTgcggc                                                                                 >  1:1033077/1‑150 (MQ=255)
                                                                        gTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCATTTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCTATAACGTCTCCCCCTTAGATAATTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCa                                                                              <  1:3107145/150‑1 (MQ=255)
                                                                          tCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAAc                                                                            >  1:4313677/1‑150 (MQ=255)
                                                                          tCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAAc                                                                            >  1:3510686/1‑150 (MQ=255)
                                                                            aTGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTGGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTg                                                                          <  1:4627258/150‑1 (MQ=255)
                                                                              gATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCt                                                                        <  2:1474986/150‑1 (MQ=255)
                                                                               aTTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTc                                                                       >  2:3437778/1‑150 (MQ=255)
                                                                                 tGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGATATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCaa                                                                     <  1:4855529/150‑1 (MQ=255)
                                                                                     tctTACCGCACAAACCCAATGTTTTTCCAGTAACTGCATACAACTTAAACTTTGCGGCCCACTATCGATAACGTCTCCCACTGATATATTTAAATCTCGCCACGTATAAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACa                                                                 <  1:128263/150‑1 (MQ=255)
                                                                                       ttACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAAc                                                               >  1:4818292/1‑150 (MQ=255)
                                                                                         aCCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCg                                                             <  2:3900696/150‑1 (MQ=255)
                                                                                         aCCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCg                                                             >  1:731718/1‑150 (MQ=255)
                                                                                         aCCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCg                                                             >  1:3427992/1‑150 (MQ=255)
                                                                                                aaaCCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTatat                                                      <  2:2625577/150‑1 (MQ=255)
                                                                                                 aaCCCATTGTTGTTCCATTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCTATAACGTCTCCCACTGATATAAGTAAATCTCGCCACGTATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATc                                                     <  1:2457563/150‑1 (MQ=255)
                                                                                                   cccAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGc                                                   >  2:1113705/1‑150 (MQ=255)
                                                                                                      attgttgtTCCTGTAACTGCAGACAACTTAAACTTTGCGGCCCACTATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATAAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAg                                                <  1:186974/148‑1 (MQ=255)
                                                                                                         gttgttCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGaaa                                             >  1:205356/1‑150 (MQ=255)
                                                                                                          ttgttCCAGTAACTGCAGACAACGTAAACTTGGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAg                                            <  1:2661602/150‑1 (MQ=255)
                                                                                                           tgttCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACTTCTCCCACTGATATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGc                                           <  1:2325457/150‑1 (MQ=255)
                                                                                                            gttCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGcc                                          >  2:4787058/1‑150 (MQ=255)
                                                                                                               ccAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGa                                       <  2:2765008/150‑1 (MQ=255)
                                                                                                                cAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGat                                      <  2:1054786/150‑1 (MQ=255)
                                                                                                                   tAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATg                                   >  1:3129658/1‑150 (MQ=255)
                                                                                                                   tAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATg                                   >  1:1266337/1‑150 (MQ=255)
                                                                                                                    aaCTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGt                                  >  2:3384494/1‑150 (MQ=255)
                                                                                                                          aGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAt                            <  2:665369/150‑1 (MQ=255)
                                                                                                                          aGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAt                            >  1:2486792/1‑150 (MQ=255)
                                                                                                                            aCAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtg                          <  2:1426253/150‑1 (MQ=255)
                                                                                                                            aCAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAGtg                          >  2:2316495/1‑150 (MQ=255)
                                                                                                                              aaCGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgat                        <  2:519823/150‑1 (MQ=255)
                                                                                                                              aaCGTAAACTTTGCGGCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgat                        <  2:4159501/150‑1 (MQ=255)
                                                                                                                              aaCGTAAACTTTGCGGCCCACGATCGATAATGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgat                        >  1:3733427/1‑150 (MQ=255)
                                                                                                                                 gTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAc                     >  1:546842/1‑150 (MQ=255)
                                                                                                                                 gTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAc                     >  1:2980149/1‑150 (MQ=255)
                                                                                                                                      cTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCa                >  1:200693/1‑150 (MQ=255)
                                                                                                                                      cTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCa                >  2:4641861/1‑150 (MQ=255)
                                                                                                                                      cTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCa                >  1:274839/1‑150 (MQ=255)
                                                                                                                                        ttGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt              <  2:150451/150‑1 (MQ=255)
                                                                                                                                           cGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtct           <  1:2788675/150‑1 (MQ=255)
                                                                                                                                             gCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtctc          >  2:363564/1‑149 (MQ=255)
                                                                                                                                             gCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtctc          >  2:361404/1‑149 (MQ=255)
                                                                                                                                             gCCCACGTTCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtctc          <  1:363564/149‑1 (MQ=255)
                                                                                                                                              cccACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTg        >  1:2284323/1‑150 (MQ=255)
                                                                                                                                                  cgatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATaa    <  1:1772515/150‑1 (MQ=255)
                                                                                                                                                  cgatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATaa    <  2:4485552/150‑1 (MQ=255)
                                                                                                                                                   gatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATaaa   >  1:2796963/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATAAACCAGTCGCCGCCATTCATCAACCACAAAGACATCTGCTGGGATGCCAGCGCATCCATCGCCATCTGTTCATGATTGCCTCTTACCGCACAAACCCAATGTTGTTCCAGTAACTGCAGACAACGTAAACTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAAC  >  CP009273/1916911‑1917208

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: