Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,522,387 C→A 89.5% 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,3870CA89.5% 426.7 / 48.1 171pseudogene (229/2037 nt)rhsEpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base C (8/10);  new base A (69/84);  total (77/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 171 total aligned reads displayed.

AACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAG  >  CP009273/1522241‑1522530
                                                                                                                                                  |                                                                                                                                               
  cGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTg                                                                                                                                              <  1:1494676/148‑1 (MQ=255)
   gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGTGTGGCTGACCTACGGCAGCGGCTATCTTGc                                                                                                                                             <  1:1006564/148‑1 (MQ=255)
   gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGc                                                                                                                                             <  1:1006038/148‑1 (MQ=255)
   gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGc                                                                                                                                             <  1:2962202/148‑1 (MQ=255)
   gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGc                                                                                                                                             <  1:2962047/148‑1 (MQ=255)
   gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACAATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGc                                                                                                                                             <  1:618408/148‑1 (MQ=255)
    ccAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCg                                                                                                                                            >  1:1118721/1‑148 (MQ=255)
    ccAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCg                                                                                                                                            >  1:2020892/1‑148 (MQ=255)
     cAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCgg                                                                                                                                           <  1:163779/148‑1 (MQ=255)
           gTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGa                                                                                                                                     <  2:1283508/148‑1 (MQ=255)
           gTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGa                                                                                                                                     <  2:1283387/148‑1 (MQ=255)
            tGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGaa                                                                                                                                    >  1:2357854/1‑148 (MQ=255)
             ggAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGTAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAg                                                                                                                                   <  2:2073136/148‑1 (MQ=255)
             ggAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAg                                                                                                                                   >  2:918607/1‑148 (MQ=255)
             ggAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAg                                                                                                                                   >  2:2572930/1‑148 (MQ=255)
                  aCCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTggg                                                                                                                              <  1:2675942/148‑1 (MQ=255)
                   cccGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGc                                                                                                                             >  2:2731452/1‑148 (MQ=255)
                      ggAGACGGGGGAACTGCTGTGGCACCTTGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGAGAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCggg                                                                                                                          <  1:1391998/148‑1 (MQ=255)
                      ggAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCggg                                                                                                                          <  1:2164580/148‑1 (MQ=255)
                       ttgaCGGGGTAACTGCTGTGGCTCCATTAGACCGGACATGCGTACTACGATCAGGGGCTGGCAAACCTCGTCACGCCGGACAGTCTGTCGGCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCgggt                                                                                                                         <  2:121471/146‑2 (MQ=37)
                       gagaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGTCTATATTGCGGGCATGAAGCTGGGCGGGa                                                                                                                         >  2:2469296/1‑148 (MQ=255)
                       gagaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGa                                                                                                                         >  2:1600567/1‑148 (MQ=255)
                       gagaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGa                                                                                                                         <  1:2616752/148‑1 (MQ=255)
                           cGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGAcgcc                                                                                                                     >  1:2991985/1‑148 (MQ=255)
                            gggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgccg                                                                                                                    <  2:2196745/148‑1 (MQ=255)
                            gggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgccg                                                                                                                    <  2:719224/148‑1 (MQ=255)
                            gggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgccg                                                                                                                    <  1:571103/148‑1 (MQ=255)
                            gggggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGAcgccg                                                                                                                    <  1:369469/148‑1 (MQ=255)
                                    tgctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGa                                                                                                            >  2:3131764/1‑148 (MQ=255)
                                    tgctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGa                                                                                                            >  1:1244416/1‑148 (MQ=255)
                                     gctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAg                                                                                                           >  1:1219575/1‑148 (MQ=255)
                                      ctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGt                                                                                                          >  1:243119/1‑148 (MQ=255)
                                      ctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGt                                                                                                          >  1:1912947/1‑148 (MQ=255)
                                      ctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGt                                                                                                          >  1:2004809/1‑148 (MQ=255)
                                         tgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCa                                                                                                       <  1:610172/148‑1 (MQ=255)
                                          ggCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAc                                                                                                      >  2:2440698/1‑148 (MQ=255)
                                               cATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAcgcgcg                                                                                                 <  2:1554592/148‑1 (MQ=255)
                                                 tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                               >  1:2775660/1‑148 (MQ=255)
                                                 tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                               >  1:1984399/1‑148 (MQ=255)
                                                    gaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGc                                                                                            <  1:1941977/148‑1 (MQ=255)
                                                     aCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGcc                                                                                           >  1:593042/1‑148 (MQ=255)
                                                     aCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGcc                                                                                           <  1:434143/148‑1 (MQ=255)
                                                          aCATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCAc                                                                                      >  1:2269056/1‑148 (MQ=255)
                                                          aCATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCAc                                                                                      >  2:140669/1‑148 (MQ=255)
                                                          aCATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCAc                                                                                      >  1:2269179/1‑148 (MQ=255)
                                                                 tACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTgaga                                                                               <  2:2994544/148‑1 (MQ=255)
                                                                  aCAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGAc                                                                              >  2:26641/1‑148 (MQ=255)
                                                                    aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                            <  1:2111409/148‑1 (MQ=255)
                                                                     aCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGt                                                                           >  2:2330348/1‑148 (MQ=255)
                                                                          cAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCa                                                                      >  2:1850448/1‑148 (MQ=255)
                                                                          cAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCa                                                                      >  2:178679/1‑148 (MQ=255)
                                                                            ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGc                                                                    <  1:1600352/148‑1 (MQ=255)
                                                                            ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGc                                                                    <  2:1909819/148‑1 (MQ=255)
                                                                             gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                   >  2:3140118/1‑148 (MQ=255)
                                                                             gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                   <  1:2373183/148‑1 (MQ=255)
                                                                             gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                   >  1:2467156/1‑148 (MQ=255)
                                                                             gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                   >  2:2869122/1‑148 (MQ=255)
                                                                             gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                   >  1:2467090/1‑148 (MQ=255)
                                                                                 tGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCgg                                                               <  2:2224030/148‑1 (MQ=255)
                                                                                   gCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                             <  1:1051428/148‑1 (MQ=255)
                                                                                    cAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgcag                                                            >  2:2835498/1‑148 (MQ=255)
                                                                                      aaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgcagca                                                          <  2:3159580/148‑1 (MQ=255)
                                                                                       aCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                         <  1:963227/148‑1 (MQ=255)
                                                                                           cgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                     <  2:2761904/148‑1 (MQ=255)
                                                                                           cgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                     <  1:2945707/148‑1 (MQ=255)
                                                                                              cACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                  >  2:2233249/1‑148 (MQ=255)
                                                                                                 gCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                               <  1:2479165/148‑1 (MQ=255)
                                                                                                 gCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                               <  1:2479149/148‑1 (MQ=255)
                                                                                                     gACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATg                                           <  1:193001/148‑1 (MQ=255)
                                                                                                      aCAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGc                                          >  1:972011/1‑148 (MQ=255)
                                                                                                      aCAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGc                                          <  2:613293/148‑1 (MQ=255)
                                                                                                        aGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCg                                        >  1:275438/1‑148 (MQ=255)
                                                                                                          tCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTCATGCCGcc                                      >  1:130286/1‑147 (MQ=255)
                                                                                                          tCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCa                                      <  2:637668/148‑1 (MQ=255)
                                                                                                               cgccgGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaa                                 <  1:1352022/148‑1 (MQ=255)
                                                                                                                gccgGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaaa                                <  1:1652955/148‑1 (MQ=255)
                                                                                                                gccgGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaaa                                <  2:1440181/148‑1 (MQ=255)
                                                                                                                gccgGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaaa                                <  2:1223636/148‑1 (MQ=255)
                                                                                                                   ggTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTg                             <  2:1868709/148‑1 (MQ=255)
                                                                                                                    gTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGa                            <  2:1185006/148‑1 (MQ=255)
                                                                                                                     tGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAc                           <  2:2766578/148‑1 (MQ=255)
                                                                                                                      ggAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAcc                          >  2:2005352/1‑148 (MQ=255)
                                                                                                                      ggAGTGGCTGACCTACGGCAGCGGCTATAGTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAcc                          <  1:2049393/148‑1 (MQ=255)
                                                                                                                         gTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGc                       <  1:1276389/148‑1 (MQ=255)
                                                                                                                          tGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGca                      >  2:1692880/1‑148 (MQ=255)
                                                                                                                          tGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGca                      >  2:1823184/1‑148 (MQ=255)
                                                                                                                          tGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGca                      >  2:3246744/1‑148 (MQ=255)
                                                                                                                             cTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAt                   <  1:1832794/148‑1 (MQ=255)
                                                                                                                               gACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATac                 <  2:1105228/148‑1 (MQ=255)
                                                                                                                                aCCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTTCACCGTGAGACTGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATaca                <  2:2358797/148‑1 (MQ=255)
                                                                                                                                aCCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATaca                >  2:1116416/1‑148 (MQ=255)
                                                                                                                                aCCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATaca                >  2:450717/1‑148 (MQ=255)
                                                                                                                                aCCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATaca                >  2:450761/1‑148 (MQ=255)
                                                                                                                                   tACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAccc             >  1:2073830/1‑148 (MQ=255)
                                                                                                                                       gCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGc         >  2:676468/1‑148 (MQ=255)
                                                                                                                                        cAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCa        <  1:539549/148‑1 (MQ=255)
                                                                                                                                        cAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCa        >  1:1569398/1‑148 (MQ=255)
                                                                                                                                           cGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGc     >  2:2091238/1‑148 (MQ=255)
                                                                                                                                            ggCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGcc    <  1:1835099/148‑1 (MQ=255)
                                                                                                                                              cTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAg  <  2:2060225/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAG  >  CP009273/1522241‑1522530

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: