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,3870CA95.4% 1225.4 / 41.6 415pseudogene (229/2037 nt)rhsEpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base C (15/4);  new base A (200/196);  total (215/200)
Fisher's exact test for biased strand distribution p-value = 1.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 415 total aligned reads displayed.

CGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTT  >  CP009273/1522239‑1522532
                                                                                                                                                    |                                                                                                                                                 
cGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCt                                                                                                                                                  >  1:1399254/1‑150 (MQ=255)
cGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCt                                                                                                                                                  >  1:2834149/1‑150 (MQ=255)
        aGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCa                                                                                                                                          >  2:3760243/1‑150 (MQ=255)
          aCGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATg                                                                                                                                        <  2:5111274/150‑1 (MQ=255)
          aCGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATg                                                                                                                                        <  2:2355137/150‑1 (MQ=255)
           cGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGa                                                                                                                                       >  1:4566156/1‑150 (MQ=255)
             gTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAg                                                                                                                                     >  2:1261351/1‑150 (MQ=255)
                 agaACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTgg                                                                                                                                 >  2:4572268/1‑150 (MQ=255)
                  gaACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTggg                                                                                                                                <  1:119619/150‑1 (MQ=255)
                    aCCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCg                                                                                                                              <  1:630422/150‑1 (MQ=255)
                    aCCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCg                                                                                                                              >  1:1006841/1‑150 (MQ=255)
                    aCCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCg                                                                                                                              <  1:3865019/150‑1 (MQ=255)
                    aCCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCg                                                                                                                              <  1:3212884/150‑1 (MQ=255)
                      ccGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCggg                                                                                                                            <  2:2669512/150‑1 (MQ=255)
                       cGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGa                                                                                                                           >  2:3679180/1‑150 (MQ=255)
                         gagaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcg                                                                                                                         >  2:1972938/1‑150 (MQ=255)
                           gaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgcc                                                                                                                       <  2:2927044/150‑1 (MQ=255)
                           gaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgcc                                                                                                                       <  1:1922861/150‑1 (MQ=255)
                             cGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGACAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgccgc                                                                                                                     >  2:1430446/1‑150 (MQ=255)
                              gggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgct                                                                                                                    >  1:4746004/1‑150 (MQ=255)
                              gggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgct                                                                                                                    <  2:2799086/150‑1 (MQ=255)
                               ggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACATCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgctg                                                                                                                   <  1:4320635/150‑1 (MQ=255)
                               ggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCAGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGActccgcag                                                                                                                   >  2:2844419/1‑148 (MQ=255)
                                gggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgctgc                                                                                                                  <  1:3758101/150‑1 (MQ=255)
                                gggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgctgc                                                                                                                  <  2:4241678/150‑1 (MQ=255)
                                          gtgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGTCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAc                                                                                                        >  2:4834379/1‑150 (MQ=255)
                                          gtgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAc                                                                                                        >  2:1063845/1‑150 (MQ=255)
                                          gtgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAc                                                                                                        >  2:1787221/1‑150 (MQ=255)
                                          gtgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAc                                                                                                        <  2:203679/150‑1 (MQ=255)
                                           tgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAcg                                                                                                       >  1:945158/1‑150 (MQ=255)
                                            ggCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAcgc                                                                                                      >  1:2605373/1‑150 (MQ=255)
                                                ccATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGa                                                                                                  <  1:3339305/150‑1 (MQ=255)
                                                 cATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                                 >  2:1401719/1‑150 (MQ=255)
                                                 cATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                                 <  1:2695810/150‑1 (MQ=255)
                                                  aTGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAcc                                                                                                >  2:3911685/1‑150 (MQ=255)
                                                   tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATAGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCg                                                                                               >  2:4909114/1‑150 (MQ=255)
                                                    gagaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGc                                                                                              <  1:2069357/150‑1 (MQ=255)
                                                    gagaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGc                                                                                              >  1:4143868/1‑150 (MQ=255)
                                                      gaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCt                                                                                            <  1:2096613/150‑1 (MQ=255)
                                                        ccGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGc                                                                                          >  1:680833/1‑150 (MQ=255)
                                                        ccGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGc                                                                                          >  1:152829/1‑150 (MQ=255)
                                                         cGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCa                                                                                         >  2:3161594/1‑150 (MQ=255)
                                                         cGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCa                                                                                         <  2:153447/150‑1 (MQ=255)
                                                           gACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCAcc                                                                                       <  1:3702183/150‑1 (MQ=255)
                                                             cATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGt                                                                                     <  1:464684/150‑1 (MQ=255)
                                                               tGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTga                                                                                   >  2:2506277/1‑150 (MQ=255)
                                                                gCGTACAACGAACAGGGGCTGGCAAACCGCGTCACACCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTgag                                                                                  <  1:2394813/150‑1 (MQ=255)
                                                                   tACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACg                                                                               >  2:1573208/1‑150 (MQ=255)
                                                                    aCAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                              >  2:3074841/1‑150 (MQ=255)
                                                                    aCAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                              <  1:706826/150‑1 (MQ=255)
                                                                    aCAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                              <  2:4175694/150‑1 (MQ=255)
                                                                     cAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGt                                                                             >  1:395714/1‑150 (MQ=255)
                                                                        cGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTgcg                                                                          <  1:3041619/150‑1 (MQ=255)
                                                                        cGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTgcg                                                                          <  1:2326564/150‑1 (MQ=255)
                                                                         gAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTgcgc                                                                         >  2:872946/1‑150 (MQ=255)
                                                                         gAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTgcgc                                                                         >  2:3238468/1‑150 (MQ=255)
                                                                         gAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTgcgc                                                                         >  2:739379/1‑150 (MQ=255)
                                                                          aaCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCa                                                                        >  2:4391840/1‑150 (MQ=255)
                                                                            cAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGc                                                                      <  1:84900/150‑1 (MQ=255)
                                                                             aGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                     >  1:495574/1‑150 (MQ=255)
                                                                                 gCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTACACCGTGAGACGGTGCGCAGCTTCgg                                                                 <  2:4519352/150‑1 (MQ=255)
                                                                                   tGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                               >  1:3068713/1‑150 (MQ=255)
                                                                                   tGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                               >  1:3446508/1‑150 (MQ=255)
                                                                                      cAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgcagca                                                            >  2:1042059/1‑150 (MQ=255)
                                                                                       aaaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                           >  1:3159008/1‑150 (MQ=255)
                                                                                       aaaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                           >  1:1606240/1‑150 (MQ=255)
                                                                                       aaaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                           <  1:3621724/150‑1 (MQ=255)
                                                                                        aaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACTGCAGCGGCTAAATTTCGGGCATGAAGCTGGGCGCGACGCCGCTGCTGGATTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTACGGCAGCATg                                                          >  1:2990826/1‑150 (MQ=18)
                                                                                           cgcgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                       >  2:4104306/1‑150 (MQ=255)
                                                                                           cgcgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                       <  2:3019765/150‑1 (MQ=255)
                                                                                              gTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                    >  1:2165569/1‑150 (MQ=255)
                                                                                               tCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggca                                                   <  1:1430446/150‑1 (MQ=255)
                                                                                               tCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGCCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggca                                                   >  1:1803906/1‑150 (MQ=255)
                                                                                                cACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGgcaggcag                                                  >  1:3124298/1‑150 (MQ=255)
                                                                                                 aCGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                                 >  1:1688464/1‑150 (MQ=255)
                                                                                                  cGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCACTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAACATGGCAGGCAGTa                                                >  2:1866318/1‑150 (MQ=255)
                                                                                                       gACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGcc                                           <  2:5058890/150‑1 (MQ=255)
                                                                                                          aGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCa                                        >  1:3987443/1‑150 (MQ=255)
                                                                                                          aGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCa                                        <  2:3988786/150‑1 (MQ=255)
                                                                                                             cTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCatat                                     <  1:975936/150‑1 (MQ=255)
                                                                                                               gccgccGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaa                                   <  2:3415122/150‑1 (MQ=255)
                                                                                                                  gccgGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACt                                <  1:4390607/150‑1 (MQ=255)
                                                                                                                   ccgGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTg                               >  2:4737695/1‑150 (MQ=255)
                                                                                                                     ggTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAc                             >  1:4110919/1‑150 (MQ=255)
                                                                                                                      gTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAcc                            <  2:2222407/150‑1 (MQ=255)
                                                                                                                      gTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAca                            <  1:2792410/150‑2 (MQ=255)
                                                                                                                       tGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCa                           >  2:811020/1‑150 (MQ=255)
                                                                                                                        ggAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAg                          >  1:2713990/1‑150 (MQ=255)
                                                                                                                         gAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGc                         >  1:4237248/1‑150 (MQ=255)
                                                                                                                           gTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcac                       <  2:495574/150‑1 (MQ=255)
                                                                                                                           gTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcac                       <  2:656831/150‑1 (MQ=255)
                                                                                                                            tGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcaca                      >  2:2963645/1‑150 (MQ=255)
                                                                                                                               cTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATac                   >  1:1128888/1‑150 (MQ=255)
                                                                                                                                tGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATaca                  <  2:3509609/150‑1 (MQ=255)
                                                                                                                                 gACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATacac                 <  1:3679180/150‑1 (MQ=255)
                                                                                                                                  aCCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAcc                <  2:4657167/150‑1 (MQ=255)
                                                                                                                                         gCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGGTGCTGGAGTTCACGCGCGAGCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGGAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAg         <  1:4834379/150‑1 (MQ=255)
                                                                                                                                         gCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAg         <  1:4880735/150‑1 (MQ=255)
                                                                                                                                         gCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAg         <  1:5106985/150‑1 (MQ=255)
                                                                                                                                            gCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGcc      <  2:1615624/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTT  >  CP009273/1522239‑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: