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,3870CA98.2% 667.5 / 1.7 217pseudogene (229/2037 nt)rhsEpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base A (84/128);  total (84/133)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 219 total aligned reads displayed.

ACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGT  >  CP009273/1522242‑1522531
                                                                                                                                                 |                                                                                                                                                
   ccAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCg                                                                                                                                             <  1:1180362/148‑1 (MQ=255)
       aCGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCa                                                                                                                                         <  1:451558/148‑1 (MQ=255)
        cGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCAt                                                                                                                                        <  1:1183483/148‑1 (MQ=255)
        cGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCAt                                                                                                                                        <  1:2961065/148‑1 (MQ=255)
           tGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGaa                                                                                                                                     <  2:685081/148‑1 (MQ=255)
           tGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGaa                                                                                                                                     <  2:3508161/148‑1 (MQ=255)
           tGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGaa                                                                                                                                     <  2:3508115/148‑1 (MQ=255)
            ggAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAg                                                                                                                                    <  2:2534656/148‑1 (MQ=255)
                aaCCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTgg                                                                                                                                >  1:3134895/1‑148 (MQ=255)
                aaCCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTgg                                                                                                                                >  2:3540807/1‑148 (MQ=255)
                     ggAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCggg                                                                                                                           <  1:2412486/148‑1 (MQ=255)
                      gagaCGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGa                                                                                                                          <  1:435117/148‑1 (MQ=255)
                              ggAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgctg                                                                                                                  <  1:1609022/148‑1 (MQ=255)
                                aaCTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCgctgct                                                                                                                <  1:2348326/148‑1 (MQ=255)
                                    gctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAg                                                                                                            <  1:2058573/148‑1 (MQ=255)
                                     ctgTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGt                                                                                                           <  2:1053784/148‑1 (MQ=255)
                                      tgtgGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGtt                                                                                                          <  1:1269538/148‑1 (MQ=255)
                                          gCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAcg                                                                                                      <  1:3057836/148‑1 (MQ=255)
                                            aCCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCAcgcg                                                                                                    <  2:1779383/148‑1 (MQ=255)
                                               aTGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGa                                                                                                 >  1:1307353/1‑148 (MQ=255)
                                                tGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAc                                                                                                >  1:1772886/1‑148 (MQ=255)
                                                 gagaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAcc                                                                                               >  2:372643/1‑148 (MQ=255)
                                                 gagaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGAcc                                                                                               >  2:2368081/1‑148 (MQ=255)
                                                   gaCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGc                                                                                             >  2:3134785/1‑148 (MQ=255)
                                                    aCCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGcc                                                                                            >  1:1068282/1‑148 (MQ=255)
                                                        gACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCa                                                                                        <  1:2027498/148‑1 (MQ=255)
                                                         aCATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCAc                                                                                       >  2:2499790/1‑148 (MQ=255)
                                                          cATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCAcc                                                                                      <  1:716917/148‑1 (MQ=255)
                                                          cATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCAcc                                                                                      >  2:2748408/1‑148 (MQ=255)
                                                            tGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGt                                                                                    >  2:3036685/1‑148 (MQ=255)
                                                            tGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGt                                                                                    >  1:1399751/1‑148 (MQ=255)
                                                            tGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGt                                                                                    <  1:1119470/148‑1 (MQ=255)
                                                             gCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTg                                                                                   <  1:3536746/148‑1 (MQ=255)
                                                                 aCAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGAc                                                                               <  2:1990155/148‑1 (MQ=255)
                                                                   aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCGGGTGGTGTGGCTGACCTACGGCAGCGGCTATTTTGCGGGCATTATGCTGGGCGGGACGCCGCTGCGGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                             <  2:188667/148‑1 (MQ=18)
                                                                   aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                             <  2:887417/148‑1 (MQ=255)
                                                                   aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                             <  1:1457107/148‑1 (MQ=255)
                                                                   aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACgg                                                                             <  1:1107644/148‑1 (MQ=255)
                                                                     cGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTg                                                                           >  1:1732195/1‑148 (MQ=255)
                                                                        aCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGTGACGGTgcgc                                                                        >  1:33629/1‑148 (MQ=255)
                                                                        aCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTgcgc                                                                        >  1:1838415/1‑148 (MQ=255)
                                                                         cAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCa                                                                       >  2:3281856/1‑148 (MQ=255)
                                                                         cAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCa                                                                       >  2:1174270/1‑148 (MQ=255)
                                                                          aGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAg                                                                      <  1:871604/148‑1 (MQ=255)
                                                                           ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTg                                                                                          >  1:3023503/1‑127 (MQ=255)
                                                                           ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGc                                                                     >  2:2770739/1‑148 (MQ=255)
                                                                           ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGc                                                                     >  2:658956/1‑148 (MQ=255)
                                                                           ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATAGTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGc                                                                     <  2:2932320/148‑1 (MQ=255)
                                                                            gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                    <  2:1631904/148‑1 (MQ=255)
                                                                            gggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCt                                                                    >  2:818395/1‑148 (MQ=255)
                                                                             ggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCtt                                                                   >  1:51203/1‑148 (MQ=255)
                                                                               cTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                 <  2:2562678/148‑1 (MQ=255)
                                                                               cTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                 <  1:1672936/148‑1 (MQ=255)
                                                                               cTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                 >  1:1409307/1‑148 (MQ=255)
                                                                                tGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCgg                                                                <  2:209434/148‑1 (MQ=255)
                                                                                  gCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                              >  1:3339551/1‑148 (MQ=255)
                                                                                  gCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                              >  1:1499665/1‑148 (MQ=255)
                                                                                     aaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgcagca                                                           <  2:355303/148‑1 (MQ=255)
                                                                                      aCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                          <  1:680069/148‑1 (MQ=255)
                                                                                       ccGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATg                                                         >  2:2236229/1‑148 (MQ=255)
                                                                                       ccGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATg                                                         <  2:282056/148‑1 (MQ=255)
                                                                                       ccGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATg                                                         >  2:1688803/1‑148 (MQ=255)
                                                                                        cgcgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATgg                                                        <  1:2779919/148‑1 (MQ=255)
                                                                                         gcgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATAGTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                       <  2:2452501/148‑1 (MQ=255)
                                                                                          cgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                      <  2:1823176/148‑1 (MQ=255)
                                                                                            tCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcagg                                                    <  2:1254312/148‑1 (MQ=255)
                                                                                            tCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcagg                                                    <  1:897688/148‑1 (MQ=255)
                                                                                             cACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                   >  2:1448687/1‑148 (MQ=255)
                                                                                             cACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                   <  2:621389/148‑1 (MQ=255)
                                                                                                     aCAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGc                                           >  1:618052/1‑148 (MQ=255)
                                                                                                        gTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGc                                        <  2:3306627/148‑1 (MQ=255)
                                                                                                          cTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCat                                      <  2:1230397/148‑1 (MQ=255)
                                                                                                           tGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCata                                     <  2:1502117/148‑1 (MQ=255)
                                                                                                           tGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCata                                     <  1:1049821/148‑1 (MQ=255)
                                                                                                              cgccgGTGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaa                                  <  2:1405083/148‑1 (MQ=255)
                                                                                                              cgccgGTGGAGTGGCTGACCTACGGCAGAGGCTATATTTCGGGCATGAAGATGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACAGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaa                                  <  1:808809/148‑1 (MQ=18)
                                                                                                               gccgGTGGAGTGGCTGACCTACGGCCGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaaa                                 >  1:1355065/1‑148 (MQ=255)
                                                                                                                    tGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAc                            >  1:3264219/1‑148 (MQ=255)
                                                                                                                    tGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAc                            >  1:1491547/1‑148 (MQ=255)
                                                                                                                    tGGAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAc                            >  1:2128209/1‑148 (MQ=255)
                                                                                                                      gAGTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCa                          >  2:3213198/1‑148 (MQ=255)
                                                                                                                        gTGGCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGc                        <  2:2287800/148‑1 (MQ=255)
                                                                                                                           gCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTTGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcaca                     <  1:3397273/148‑1 (MQ=255)
                                                                                                                           gCTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcaca                     <  1:2014784/148‑1 (MQ=255)
                                                                                                                            cTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAt                    <  1:1571923/148‑1 (MQ=255)
                                                                                                                            cTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAt                    <  1:378138/148‑1 (MQ=255)
                                                                                                                            cTGACCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAt                    <  1:1893142/148‑1 (MQ=255)
                                                                                                                              gACCTACGGCAGCGGCTATAGTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATac                  <  1:42844/148‑1 (MQ=255)
                                                                                                                               aCCTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATaca                 <  2:1143345/148‑1 (MQ=255)
                                                                                                                                ccTACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCTCATATAAACTGACCAGCACATacac                <  2:219437/148‑1 (MQ=255)
                                                                                                                                  tACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAccc              <  1:2236229/148‑1 (MQ=255)
                                                                                                                                  tACGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAccc              >  2:2616900/1‑148 (MQ=255)
                                                                                                                                    cGGCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAccccc            >  1:2512400/1‑148 (MQ=255)
                                                                                                                                      gCAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGc          <  2:105958/148‑1 (MQ=255)
                                                                                                                                       cAGCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCa         <  1:1485361/148‑1 (MQ=255)
                                                                                                                                         gCGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAgg       <  2:1307353/148‑1 (MQ=255)
                                                                                                                                          cGGCTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGc      <  2:2888028/148‑1 (MQ=255)
                                                                                                                                            gCTATATTTCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCa    >  1:762250/1‑148 (MQ=255)
                                                                                                                                             cTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAg   <  1:2368081/148‑1 (MQ=255)
                                                                                                                                             cTATATTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAg   >  1:2514659/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGT  >  CP009273/1522242‑1522531

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: