Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,695 C→A 12.7% G24C (GGT→TGT)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,851,6950CA12.7% 707.5 / 49.4 221G24C (GGT→TGT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (96/95);  new base A (7/22);  total (103/118)
Fisher's exact test for biased strand distribution p-value = 9.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

TGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACA  >  CP009273/3851548‑3851842
                                                                                                                                                   |                                                                                                                                                   
   tCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCt                                                                                                                                                  <  1:939160/148‑1 (MQ=255)
         ggCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCa                                                                                                                                            >  1:2299642/1‑148 (MQ=255)
          gCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCaa                                                                                                                                           <  2:1152978/148‑1 (MQ=255)
                  ggTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCat                                                                                                                                   <  1:223143/148‑1 (MQ=255)
                  ggTTAGCGTAGATACCGGTAAATAGTCCGGCCCGCGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCat                                                                                                                                   <  2:2952941/148‑1 (MQ=255)
                   gTTAGCTTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCata                                                                                                                                  <  1:846484/148‑1 (MQ=255)
                   gTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCata                                                                                                                                  >  1:2745655/1‑148 (MQ=255)
                     tAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATAtt                                                                                                                                >  1:957317/1‑148 (MQ=255)
                     tAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATAtt                                                                                                                                >  1:2852376/1‑148 (MQ=255)
                     tAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATAtt                                                                                                                                >  1:1654470/1‑148 (MQ=255)
                      aGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtg                                                                                                                               <  1:357965/148‑1 (MQ=255)
                       gCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtgg                                                                                                                              <  2:2612350/148‑1 (MQ=255)
                          tAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtggtgg                                                                                                                           >  1:1110871/1‑148 (MQ=255)
                            gATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGcc                                                                                                                         >  2:846454/1‑148 (MQ=255)
                             aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCt                                                                                                                        >  1:2585444/1‑148 (MQ=255)
                                   gTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTAt                                                                                                                  <  1:1964331/148‑1 (MQ=255)
                                    tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGACGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATc                                                                                                                 <  1:1766535/148‑1 (MQ=255)
                                       aTAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGt                                                                                                              >  1:742956/1‑148 (MQ=255)
                                        tAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCCTATTGGTGGCCTGGGTCTCGGTc                                                                                                             >  2:3046036/1‑148 (MQ=255)
                                        tAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTc                                                                                                             <  2:56377/148‑1 (MQ=255)
                                          gTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCAt                                                                                                           >  2:2911743/1‑148 (MQ=255)
                                           tCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATg                                                                                                          >  2:72732/1‑148 (MQ=255)
                                           tCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATg                                                                                                          <  2:2281156/148‑1 (MQ=255)
                                           tCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATg                                                                                                          <  2:1251254/148‑1 (MQ=255)
                                           tCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATg                                                                                                          >  2:1210544/1‑148 (MQ=255)
                                               gCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAAATAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGa                                                                                                      <  2:1153967/148‑1 (MQ=255)
                                                cccGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGaa                                                                                                     >  1:1504250/1‑148 (MQ=255)
                                                 ccGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAAc                                                                                                    >  2:1941378/1‑148 (MQ=255)
                                                  cGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACa                                                                                                   >  1:3039981/1‑148 (MQ=255)
                                                    tGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGa                                                                                                 >  1:2145496/1‑148 (MQ=255)
                                                      cAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGaaa                                                                                               <  1:842293/148‑1 (MQ=255)
                                                        ggCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAt                                                                                             >  1:3163144/1‑148 (MQ=255)
                                                        ggCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAt                                                                                             >  2:1356963/1‑148 (MQ=255)
                                                         gCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAtt                                                                                            >  2:322181/1‑148 (MQ=255)
                                                         gCGTACAAACCGGAGAACAGGTGTAGGCGGAATAAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAACtt                                                                                            >  2:2539879/1‑148 (MQ=255)
                                                          cGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTg                                                                                           <  2:631312/148‑1 (MQ=255)
                                                           gTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTgg                                                                                          <  1:693113/148‑1 (MQ=255)
                                                             aCAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTgggg                                                                                        >  2:1140674/1‑148 (MQ=255)
                                                               aaaCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCg                                                                                      <  2:616358/148‑1 (MQ=255)
                                                                 aCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGtt                                                                                    >  2:1660652/1‑148 (MQ=255)
                                                                  ccGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGttt                                                                                   >  1:3234026/1‑148 (MQ=255)
                                                                  ccGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTAGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGttt                                                                                   <  1:3063286/148‑1 (MQ=255)
                                                                   cGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTc                                                                                  <  1:424209/148‑1 (MQ=255)
                                                                      tGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATa                                                                               >  1:2529612/1‑148 (MQ=255)
                                                                      tGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATa                                                                               <  2:2469889/148‑1 (MQ=255)
                                                                       gAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAg                                                                              >  1:1105630/1‑148 (MQ=255)
                                                                         tcAGGTTTAGGCGTAATTAATGAGAATAACTTCCGCCGACAGTCTTTCAATATTTTGCGTATTCTTCGGTTTACCGATATAGAAACCGTCCTTATTGGTGGCCTGGGTATCGGTCATGACGAACTGTATATTGGGGCGTTTGATATGc                                                                            <  2:1867187/147‑1 (MQ=255)
                                                                         aCAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGc                                                                            <  1:488506/148‑1 (MQ=255)
                                                                         aCAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCAAGACGAACAGAAAATTGGGGCGTTTCATAGGc                                                                            <  1:1332256/148‑1 (MQ=255)
                                                                          cAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCt                                                                           <  1:2243167/148‑1 (MQ=255)
                                                                          cAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCt                                                                           >  2:2216299/1‑148 (MQ=255)
                                                                            ggTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTAt                                                                         <  1:1730756/148‑1 (MQ=255)
                                                                                 aGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTc                                                                    <  2:2813878/148‑1 (MQ=255)
                                                                                 aGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTc                                                                    <  1:2382554/148‑1 (MQ=255)
                                                                                 aGGCGGAATTAAAGACTATCCCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATTTTGGTGGCCTTGGTATCGTTCATGACGAACAGACAATTGGGGCGTTTCATAGGCTATCCTTc                                                                    <  1:587220/148‑1 (MQ=255)
                                                                                  ggCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTcc                                                                   <  2:425779/148‑1 (MQ=255)
                                                                                       aaTTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCg                                                              <  1:1323374/148‑1 (MQ=255)
                                                                                        aTTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCgg                                                             >  1:3118892/1‑148 (MQ=255)
                                                                                           aaaGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAg                                                          <  2:1109204/148‑1 (MQ=255)
                                                                                              gCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTaa                                                       <  1:229103/148‑1 (MQ=255)
                                                                                               cGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTaaa                                                      <  1:1232980/148‑1 (MQ=255)
                                                                                                 aaTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAgg                                                    >  1:1055122/1‑148 (MQ=255)
                                                                                                  aTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGt                                                   >  2:463972/1‑148 (MQ=255)
                                                                                                   tACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGtt                                                  >  2:1737459/1‑148 (MQ=255)
                                                                                                    aCCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTa                                                 <  1:671073/148‑1 (MQ=255)
                                                                                                      cTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACg                                               >  2:3049609/1‑148 (MQ=255)
                                                                                                      cTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACg                                               <  1:101808/148‑1 (MQ=255)
                                                                                                       ttCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACgg                                              <  1:309698/148‑1 (MQ=255)
                                                                                                       ttCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACgg                                              <  1:1976283/148‑1 (MQ=255)
                                                                                                       ttCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACgg                                              >  2:1729898/1‑148 (MQ=255)
                                                                                                       ttCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACgg                                              <  1:2442118/148‑1 (MQ=255)
                                                                                                        tCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATCGCACCCGACCAGCGGGGTGGCCTGGGTAGCGGTCATGACGCACAGAAGGTGGGGGCGGGTCAGCGGCGCGGCTTCCCGGCGGGGGTAGAgggtggggg                                             >  2:1903827/1‑141 (MQ=255)
                                                                                                         ccgccgCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCt                                            >  2:821374/1‑148 (MQ=255)
                                                                                                          cgccgcCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTg                                           <  1:2109909/148‑1 (MQ=255)
                                                                                                               ccAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCg                                      >  1:1691953/1‑148 (MQ=255)
                                                                                                                cAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGt                                     <  1:1210544/148‑1 (MQ=255)
                                                                                                                 aGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTg                                    <  2:355547/148‑1 (MQ=255)
                                                                                                                  gACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGACCAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGa                                   <  2:1778102/148‑1 (MQ=255)
                                                                                                                  gACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGa                                   <  1:2257450/148‑1 (MQ=255)
                                                                                                                    cTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACg                                 <  1:2730915/148‑1 (MQ=255)
                                                                                                                      aTCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGAt                               >  2:509799/1‑148 (MQ=255)
                                                                                                                      aTCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGAt                               >  1:2814438/1‑148 (MQ=255)
                                                                                                                       tCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATg                              >  2:3182978/1‑148 (MQ=255)
                                                                                                                       tCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATg                              >  2:361818/1‑148 (MQ=255)
                                                                                                                        cAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgc                             <  2:2724289/148‑1 (MQ=255)
                                                                                                                        cAATATTTTGCGTATTCAGCGGTTTACCGCGATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGACCAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGGTgc                             >  2:1262624/1‑148 (MQ=255)
                                                                                                                           tatTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgccgc                          >  1:1926127/1‑148 (MQ=255)
                                                                                                                           tatTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgccgc                          <  2:2080429/148‑1 (MQ=255)
                                                                                                                              tttGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAg                       <  1:3184020/148‑1 (MQ=255)
                                                                                                                               ttGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGc                      <  1:959741/148‑1 (MQ=255)
                                                                                                                                    tATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCa                 >  1:2961733/1‑148 (MQ=255)
                                                                                                                                    tATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCa                 >  1:2245375/1‑148 (MQ=255)
                                                                                                                                       tCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCaaaa              >  1:1480628/1‑148 (MQ=255)
                                                                                                                                         aGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAAt            <  1:1851420/148‑1 (MQ=255)
                                                                                                                                            ggTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAg         <  1:2911743/148‑1 (MQ=255)
                                                                                                                                            ggTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAg         >  2:1876635/1‑148 (MQ=255)
                                                                                                                                                tACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAg     <  2:2663251/148‑1 (MQ=255)
                                                                                                                                                 aCCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGATAGa    <  1:509799/148‑1 (MQ=255)
                                                                                                                                                  ccGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGAc   >  2:2201283/1‑148 (MQ=255)
                                                                                                                                                   cGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACa  <  2:742956/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGACA  >  CP009273/3851548‑3851842

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: ‑