Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,328,313 C→A 71.3% R68S (CGT→AGT)  cysB → N‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,328,3130CA71.3% 209.1 / 159.5 177R68S (CGT→AGT) cysBN‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base A (65/61);  total (86/91)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

GTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCTTTATT  >  CP009273/1328166‑1328459
                                                                                                                                                   |                                                                                                                                                  
  ctcATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGt                                                                                                                                                  <  2:2792025/148‑1 (MQ=255)
  ctcATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGt                                                                                                                                                  >  2:2133407/1‑148 (MQ=255)
  ctcATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGt                                                                                                                                                  <  1:1065953/148‑1 (MQ=255)
   tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTa                                                                                                                                                 <  2:221398/148‑1 (MQ=255)
   tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTa                                                                                                                                                 <  1:1745013/148‑1 (MQ=255)
    catcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTAt                                                                                                                                                <  1:2582858/148‑1 (MQ=255)
       caACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTAtcgc                                                                                                                                             >  1:2978438/1‑148 (MQ=255)
        aaCAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTAtcgct                                                                                                                                            <  2:823211/148‑1 (MQ=255)
          cAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTAtcgctcg                                                                                                                                          <  1:798282/148‑1 (MQ=255)
            gCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTcgcg                                                                                                                                        >  2:500476/1‑148 (MQ=255)
               gAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCACGAAATACTTCGTATCGCTCGCGAAg                                                                                                                                     >  2:2374653/1‑148 (MQ=255)
               gAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGACATAATTCGTATCGCTCGCGAAg                                                                                                                                     >  2:3382311/1‑148 (MQ=255)
                 aGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTc                                                                                                                                   <  2:570782/148‑1 (MQ=255)
                 aGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTc                                                                                                                                   >  2:2048977/1‑148 (MQ=255)
                  ggACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTcc                                                                                                                                  <  1:167413/148‑1 (MQ=255)
                   gACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATATTTTCCCGAAGCGGCTAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCt                                                                                                                                 <  2:1632747/148‑1 (MQ=255)
                    aCTTTACACATCACCACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTg                                                                                                                                <  2:2494200/148‑1 (MQ=255)
                    aCTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTg                                                                                                                                <  2:1453362/148‑1 (MQ=255)
                    aCTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTg                                                                                                                                <  2:3275509/148‑1 (MQ=255)
                          cacaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGaaa                                                                                                                          <  1:1313941/148‑1 (MQ=255)
                           acaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAg                                                                                                                         <  2:119053/148‑1 (MQ=255)
                            caTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGt                                                                                                                        >  1:3501356/1‑148 (MQ=255)
                              tCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCg                                                                                                                      <  2:1899234/148‑1 (MQ=255)
                              tCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCg                                                                                                                      <  2:2956604/148‑1 (MQ=255)
                               cacaACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGa                                                                                                                     >  2:2620248/1‑148 (MQ=255)
                                    cccGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCa                                                                                                                >  1:2900672/1‑148 (MQ=255)
                                     ccGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCAt                                                                                                               <  1:405646/148‑1 (MQ=255)
                                      cGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATa                                                                                                              <  2:1124998/148‑1 (MQ=255)
                                        ggATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATaaa                                                                                                            <  2:667396/148‑1 (MQ=255)
                                           tCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATc                                                                                                         <  1:772444/148‑1 (MQ=255)
                                            cAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCg                                                                                                        <  2:1313374/148‑1 (MQ=255)
                                            cAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCg                                                                                                        <  2:76331/148‑1 (MQ=255)
                                            cAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCg                                                                                                        <  2:2978438/148‑1 (MQ=255)
                                             aGTAAAAAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCgg                                                                                                       <  1:720641/148‑1 (MQ=255)
                                              gTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGt                                                                                                      >  2:752480/1‑148 (MQ=255)
                                                aaaCAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTg                                                                                                    <  1:2668806/148‑1 (MQ=255)
                                                aaaCAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTg                                                                                                    >  1:3469833/1‑148 (MQ=255)
                                                 aaCAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGc                                                                                                   <  1:2138941/148‑1 (MQ=255)
                                                  aCAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGcc                                                                                                  >  2:2239676/1‑148 (MQ=255)
                                                   cAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCg                                                                                                 >  2:1930500/1‑148 (MQ=255)
                                                      gTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGgag                                                                                              <  1:2951618/148‑1 (MQ=255)
                                                          gAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGca                                                                                          <  2:78567/148‑1 (MQ=255)
                                                           aaTGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGcac                                                                                         >  1:2234293/1‑148 (MQ=255)
                                                            aTGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGcaca                                                                                        <  1:476416/148‑1 (MQ=255)
                                                             tGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCacac                                                                                       <  1:3208233/148‑1 (MQ=255)
                                                                 ggAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTgg                                                                                   <  1:1927140/148‑1 (MQ=255)
                                                                  gAAGACGAGCTAGGCCTTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGc                                                                                  >  2:2157820/1‑148 (MQ=255)
                                                                  gAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGc                                                                                  >  2:2157814/1‑148 (MQ=255)
                                                                  gAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGc                                                                                  >  2:1509570/1‑148 (MQ=255)
                                                                   aaGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTggtc                                                                                 <  2:290896/148‑3 (MQ=255)
                                                                   aaGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTggtc                                                                                 <  2:2063318/148‑3 (MQ=255)
                                                                    aGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCg                                                                                >  1:963050/1‑148 (MQ=255)
                                                                      aCGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGa                                                                              >  1:2717920/1‑148 (MQ=255)
                                                                            tAGGCATACAAATTTTGTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAgg                                                                        <  2:2588887/148‑1 (MQ=255)
                                                                               gCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTc                                                                     >  2:228631/1‑148 (MQ=255)
                                                                               gCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTc                                                                     >  1:238322/1‑148 (MQ=255)
                                                                                cATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCa                                                                    >  1:144977/1‑148 (MQ=255)
                                                                                  ttCAAACTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACt                                                                  >  1:593713/1‑148 (MQ=255)
                                                                                   tCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTg                                                                 >  1:220961/1‑148 (MQ=255)
                                                                                    cAATTTTTTTCCCGAAGCGGCAAGCACCTGCCGCAGGTACCGCCCGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAATGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGt                                                                <  2:2128182/148‑1 (MQ=255)
                                                                                      aaTTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGtat                                                              <  2:941185/148‑1 (MQ=255)
                                                                                        ttttttCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTatat                                                            <  2:1266479/148‑1 (MQ=255)
                                                                                                cGAAGCGGCAAGCACCCGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGccacca                                                    <  2:149716/148‑1 (MQ=255)
                                                                                                  aaGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACg                                                  <  2:2911671/148‑1 (MQ=255)
                                                                                                  aaGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACg                                                  >  2:3418416/1‑148 (MQ=255)
                                                                                                   aGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGc                                                 >  2:758518/1‑148 (MQ=255)
                                                                                                    gCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCa                                                >  1:448314/1‑148 (MQ=255)
                                                                                                     cGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCAt                                               <  1:544697/148‑1 (MQ=255)
                                                                                                           gCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAg                                         >  1:2879612/1‑148 (MQ=255)
                                                                                                            cACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCTCCACCACGCATACCCAgg                                        >  1:324278/1‑148 (MQ=255)
                                                                                                                tGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACg                                    >  1:1675005/1‑148 (MQ=255)
                                                                                                                 gACGCAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGc                                   >  1:2875189/1‑148 (MQ=255)
                                                                                                                     cAGGTAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACg                               >  2:2115466/1‑148 (MQ=255)
                                                                                                                     cAGGAAACGACAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACg                               >  2:964984/1‑148 (MQ=255)
                                                                                                                       ggTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCa                             <  1:2032846/148‑1 (MQ=255)
                                                                                                                       ggTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCa                             <  2:3344246/148‑1 (MQ=255)
                                                                                                                       ggTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCa                             <  1:858074/148‑1 (MQ=255)
                                                                                                                         tAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCAtt                           >  2:2504171/1‑148 (MQ=255)
                                                                                                                         tAACGCCAGCAGGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCAtt                           >  1:295710/1‑148 (MQ=255)
                                                                                                                          aaCGCCAGCAGGGCAAGAAATAATTAGTCTCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTa                          <  2:3399550/148‑1 (MQ=255)
                                                                                                                                agcagGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAAc                    >  2:3365503/1‑148 (MQ=255)
                                                                                                                                 gcagGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACg                   <  2:719574/148‑1 (MQ=255)
                                                                                                                                 gcagGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACg                   >  1:2970163/1‑148 (MQ=255)
                                                                                                                                  cagGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGt                  <  2:2027023/148‑1 (MQ=255)
                                                                                                                                  cagGGCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGt                  >  2:3009951/1‑148 (MQ=255)
                                                                                                                                    gggCAAGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGtca                >  2:2607892/1‑148 (MQ=255)
                                                                                                                                       cAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGtcatca             >  1:209715/1‑148 (MQ=255)
                                                                                                                                        aaGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCaa            <  2:3148574/148‑1 (MQ=255)
                                                                                                                                         aGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCaaa           <  1:1153483/148‑1 (MQ=255)
                                                                                                                                         aGAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCaaa           >  2:238204/1‑148 (MQ=255)
                                                                                                                                          gAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAg          >  2:1965542/1‑148 (MQ=255)
                                                                                                                                          gAAATAATTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAg          <  1:2040099/148‑1 (MQ=255)
                                                                                                                                            aataatTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGc        >  1:1938260/1‑148 (MQ=255)
                                                                                                                                             ataatTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCt       >  2:1188139/1‑148 (MQ=255)
                                                                                                                                             ataatTAGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCt       <  1:2248278/148‑1 (MQ=255)
                                                                                                                                              taatTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCtt      <  1:475264/148‑1 (MQ=255)
                                                                                                                                                atTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCTTTa    <  2:2689962/148‑1 (MQ=255)
                                                                                                                                                 ttCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCTTTAt   <  1:2620248/148‑1 (MQ=255)
                                                                                                                                                 ttCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCTTTAt   <  1:2620208/148‑1 (MQ=255)
                                                                                                                                                  tCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCTTTAtt  <  2:1919418/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCTGTCGAAAGTCGATGCCATAAAATCGGTTGCCGGAGAGCACACCTGGCCGGATAAAGGTTCACTGTATATCGCCACCACGCATACCCAGGCACGCTACGCATTACCAAACGTCATCAAAGGCTTTATT  >  CP009273/1328166‑1328459

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: