Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,522,387 C→A 100% pseudogene (229/2037 nt) rhsE → pseudogene, Rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,522,3870CA97.3% 756.3 / 8.8 247pseudogene (229/2037 nt)rhsEpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (75/164);  total (78/169)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 250 total aligned reads displayed.

AACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTT  >  CP009273/1522241‑1522532
                                                                                                                                                  |                                                                                                                                                 
aaCGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATAtt                                                                                                                                                 >  2:293838/1‑149 (MQ=255)
 aCGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTg                                                                                                                                                >  1:3710951/1‑149 (MQ=255)
 aCGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTg                                                                                                                                                >  1:1922378/1‑149 (MQ=255)
  cGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGc                                                                                                                                               <  1:1058040/149‑1 (MQ=255)
  cGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGc                                                                                                                                               <  2:2764639/149‑1 (MQ=255)
  cGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGc                                                                                                                                               <  2:513569/149‑1 (MQ=255)
  cGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGc                                                                                                                                               <  1:2117854/149‑1 (MQ=255)
   gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCg                                                                                                                                              >  1:2532377/1‑149 (MQ=255)
    ccAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCgg                                                                                                                                             <  2:3312197/149‑1 (MQ=255)
     cAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCggg                                                                                                                                            <  2:108323/149‑1 (MQ=255)
      aGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGc                                                                                                                                           >  2:2771861/1‑149 (MQ=255)
           gTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGTGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGaa                                                                                                                                      <  2:310067/149‑1 (MQ=255)
           gTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGaa                                                                                                                                      <  2:2879884/149‑1 (MQ=255)
             ggAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGc                                                                                                                                    <  2:3733224/149‑1 (MQ=255)
             ggAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGc                                                                                                                                    <  2:550270/149‑1 (MQ=255)
              gagaACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCt                                                                                                                                   >  1:2213863/1‑149 (MQ=255)
               agaACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGTCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATGTTGCGGGCATGAAGCTg                                                                                                                                  <  2:1026692/149‑1 (MQ=255)
               agaACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTg                                                                                                                                  <  2:1315029/149‑1 (MQ=255)
                      ggAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGa                                                                                                                           <  1:914954/149‑1 (MQ=255)
                       gagaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAc                                                                                                                          >  2:1735705/1‑149 (MQ=255)
                         gaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgc                                                                                                                        <  1:65708/149‑1 (MQ=255)
                         gaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgc                                                                                                                        >  1:4056290/1‑149 (MQ=255)
                             ggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgct                                                                                                                    <  1:56277/149‑1 (MQ=255)
                              gggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgctg                                                                                                                   <  1:2290511/149‑1 (MQ=255)
                                 aaCTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGctgctg                                                                                                                >  2:3474328/1‑149 (MQ=255)
                                     gctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGt                                                                                                            >  1:42519/1‑149 (MQ=255)
                                     gctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGt                                                                                                            <  1:4073756/149‑1 (MQ=255)
                                      ctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGtt                                                                                                           <  1:2385311/149‑1 (MQ=255)
                                        gtgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCa                                                                                                         <  2:2522212/149‑1 (MQ=255)
                                        gtgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCa                                                                                                         <  1:800197/149‑1 (MQ=255)
                                         tgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAc                                                                                                        <  1:623471/149‑1 (MQ=255)
                                         tgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAc                                                                                                        <  1:2687570/149‑1 (MQ=255)
                                                aTGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                                 <  1:3505783/149‑1 (MQ=255)
                                                aTGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                                 <  2:16440/149‑1 (MQ=255)
                                                aTGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                                 <  1:768997/149‑1 (MQ=255)
                                                aTGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                                 <  2:16688/149‑1 (MQ=255)
                                                 tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAcc                                                                                                >  1:1226779/1‑149 (MQ=255)
                                                 tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAcc                                                                                                <  1:1225211/149‑1 (MQ=255)
                                                 tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAcc                                                                                                >  1:3273165/1‑149 (MQ=255)
                                                 tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAcc                                                                                                >  1:836588/1‑149 (MQ=255)
                                                  gagaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCg                                                                                               >  1:3035912/1‑149 (MQ=255)
                                                    gaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGcc                                                                                             <  2:1172080/149‑1 (MQ=255)
                                                    gaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGcc                                                                                             <  2:3795750/149‑1 (MQ=255)
                                                     aCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCt                                                                                            <  2:4050880/149‑1 (MQ=255)
                                                     aCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCt                                                                                            <  1:825460/149‑1 (MQ=255)
                                                     aCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCt                                                                                            <  2:827885/149‑1 (MQ=255)
                                                            aTGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGt                                                                                     >  2:1623268/1‑149 (MQ=255)
                                                            aTGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGt                                                                                     <  1:136526/149‑1 (MQ=255)
                                                            aTGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGt                                                                                     <  2:1117350/149‑1 (MQ=255)
                                                             tGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTg                                                                                    >  1:1677125/1‑149 (MQ=255)
                                                              gCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTga                                                                                   <  1:1286944/149‑1 (MQ=255)
                                                              gCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTga                                                                                   <  1:3019494/149‑1 (MQ=255)
                                                                gTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTgaga                                                                                 <  2:652387/149‑1 (MQ=255)
                                                                 tACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGAc                                                                                <  2:235484/149‑1 (MQ=255)
                                                                  aCAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACg                                                                               <  2:889414/149‑1 (MQ=255)
                                                                  aCAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACg                                                                               <  2:3328387/149‑1 (MQ=255)
                                                                   cAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                              >  1:3341725/1‑149 (MQ=255)
                                                                    aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGt                                                                             <  2:1400177/149‑1 (MQ=255)
                                                                    aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGt                                                                             <  1:2887343/149‑1 (MQ=255)
                                                                            ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                     >  1:981002/1‑149 (MQ=255)
                                                                            ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                     >  2:1534548/1‑149 (MQ=255)
                                                                             gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCtt                                                                    <  1:1318966/149‑1 (MQ=255)
                                                                             gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACAGTGAGACGGTGCGCAGCtt                                                                    >  1:2797148/1‑149 (MQ=255)
                                                                                 tGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgc                                                                <  2:3484158/149‑1 (MQ=255)
                                                                                  ggCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                               <  2:1965731/149‑1 (MQ=255)
                                                                                     aaaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgcagca                                                            <  1:2171644/149‑1 (MQ=255)
                                                                                               aCGCCGGACAGTCTGCCGCCGGTGGTGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGCGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGgcaggcag                                                  >  2:3679737/1‑149 (MQ=255)
                                                                                                cGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                                 >  1:1768849/1‑149 (MQ=255)
                                                                                                 gCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTa                                                <  1:109700/149‑1 (MQ=255)
                                                                                                 gCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTa                                                >  1:148495/1‑149 (MQ=255)
                                                                                                  ccGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTaa                                               <  2:3170652/149‑1 (MQ=255)
                                                                                                   cGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAAt                                              <  1:2225992/149‑1 (MQ=255)
                                                                                                        aGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGc                                         >  1:3502669/1‑149 (MQ=255)
                                                                                                         gTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCa                                        >  1:2601472/1‑149 (MQ=255)
                                                                                                         gTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCa                                        >  1:3172632/1‑149 (MQ=255)
                                                                                                         gTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCa                                        <  1:1578876/149‑1 (MQ=255)
                                                                                                           cTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCata                                      <  2:1098716/149‑1 (MQ=255)
                                                                                                             gccgccGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAtata                                    >  2:1504390/1‑149 (MQ=255)
                                                                                                             gccgccGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAtata                                    >  2:1357472/1‑149 (MQ=255)
                                                                                                             gccgccGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAgata                                    >  2:2185150/1‑149 (MQ=255)
                                                                                                                   ggTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGa                              <  1:2452653/149‑1 (MQ=255)
                                                                                                                      ggAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCa                           >  2:393269/1‑149 (MQ=255)
                                                                                                                      ggAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCa                           >  2:2461755/1‑149 (MQ=255)
                                                                                                                         gtgtCTGACCTACGGCATCGGCTATATTGCGGGCATGTAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCTCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGca                        <  2:1892236/149‑1 (MQ=18)
                                                                                                                         gTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGca                        <  2:3213100/149‑1 (MQ=255)
                                                                                                                           ggCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcaca                      <  2:4056420/149‑1 (MQ=255)
                                                                                                                            gCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAt                     <  1:522723/149‑1 (MQ=255)
                                                                                                                            gCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATAGAAACTGACCAGCACAt                     >  2:2386093/1‑149 (MQ=255)
                                                                                                                             cTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATa                    <  2:2134775/149‑1 (MQ=255)
                                                                                                                             cTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATa                    <  2:248298/149‑1 (MQ=255)
                                                                                                                              tGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATac                   >  2:914683/1‑149 (MQ=255)
                                                                                                                               gACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATaca                  <  1:140631/149‑1 (MQ=255)
                                                                                                                                aCCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATacac                 <  1:3302135/149‑1 (MQ=255)
                                                                                                                                 ccTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAcc                <  1:1895532/149‑1 (MQ=255)
                                                                                                                                   tACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAcccc              <  1:3608136/149‑1 (MQ=255)
                                                                                                                                       gCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCa          <  2:3699286/149‑1 (MQ=255)
                                                                                                                                          gCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGc       >  2:1951314/1‑149 (MQ=255)
                                                                                                                                             gCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAg    >  2:2566792/1‑149 (MQ=255)
                                                                                                                                              cTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGt   <  2:3035912/149‑1 (MQ=255)
                                                                                                                                               tataTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGtt  <  2:3536143/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTT  >  CP009273/1522241‑1522532

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: