Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 923,339 G→T 29.0% S438* (TCG→TAG)  cydC ← fused glutathione, cysteine exporter subunits of ABC superfamily: membrane component/ATP‑binding component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273923,3390GT29.0% 718.2 / 293.2 341S438* (TCG→TAG) cydCfused glutathione, cysteine exporter subunits of ABC superfamily: membrane component/ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base G (120/122);  new base T (44/55);  total (164/177)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 363 total aligned reads displayed.

CGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGCCCTGTTGCGGGTCCCATGC  >  CP009273/923192‑923487
                                                                                                                                                   |                                                                                                                                                    
  ggCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCTACAGAGCCTCATCACTACTGCCaggcgagg                                                                                                                                                   <  1:2017414/149‑1 (MQ=255)
  ggCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCaggcgagg                                                                                                                                                   >  2:4296141/1‑149 (MQ=255)
   gCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCaggcgaggc                                                                                                                                                  <  2:1248970/149‑1 (MQ=255)
     gATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAggcgaggcga                                                                                                                                                <  1:2509762/149‑1 (MQ=255)
      aTAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGgcgaggcgag                                                                                                                                               >  2:2791439/1‑149 (MQ=255)
         gCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTaa                                                                                                                                            >  2:4593316/1‑149 (MQ=255)
            aGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAg                                                                                                                                         >  2:807375/1‑149 (MQ=255)
             gACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGa                                                                                                                                        >  2:4979631/1‑149 (MQ=255)
             gACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGa                                                                                                                                        >  1:4103720/1‑149 (MQ=255)
                 gcgcAGTTCTCCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTAt                                                                                                                                    <  2:3975915/149‑1 (MQ=255)
                 gcgcAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTAt                                                                                                                                    >  2:857282/1‑149 (MQ=255)
                   gcAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCa                                                                                                                                  <  2:3523566/149‑1 (MQ=255)
                   gcAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCa                                                                                                                                  <  2:103829/149‑1 (MQ=255)
                   gcAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCa                                                                                                                                  <  1:712476/149‑1 (MQ=255)
                          accaccGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCg                                                                                                                           >  2:2907680/1‑149 (MQ=255)
                           ccaccGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGt                                                                                                                          <  2:1816654/149‑1 (MQ=255)
                             accGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTgg                                                                                                                        >  1:5189021/1‑149 (MQ=255)
                              ccGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGgc                                                                                                                       <  2:47822/149‑1 (MQ=255)
                                 gagaGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCt                                                                                                                    <  2:4114888/149‑1 (MQ=255)
                                 gagaGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCt                                                                                                                    >  2:5138413/1‑149 (MQ=255)
                                   gagCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTaa                                                                                                                  >  1:1627613/1‑149 (MQ=255)
                                     gCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCTACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAAc                                                                                                                <  2:563358/149‑1 (MQ=255)
                                      cTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACa                                                                                                               <  1:481243/149‑1 (MQ=255)
                                      cTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACa                                                                                                               <  1:2342178/149‑1 (MQ=255)
                                        ggCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGa                                                                                                             >  1:4303801/1‑149 (MQ=255)
                                        ggCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGa                                                                                                             <  2:1216073/149‑1 (MQ=255)
                                              cGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACt                                                                                                       >  2:4363291/1‑149 (MQ=255)
                                               gCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTc                                                                                                      >  1:3381943/1‑149 (MQ=255)
                                               gCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTc                                                                                                      >  1:3247264/1‑149 (MQ=255)
                                                     aCCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCTACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAg                                                                                                >  1:5110430/1‑149 (MQ=255)
                                                       cTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGa                                                                                              >  1:3674722/1‑149 (MQ=255)
                                                         aaCCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGaac                                                                                            <  2:4080867/149‑1 (MQ=255)
                                                           ccAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGaacaa                                                                                          >  1:2367901/1‑149 (MQ=255)
                                                              aCTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGc                                                                                       >  2:370237/1‑149 (MQ=255)
                                                                 gTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGa                                                                                    >  1:657681/1‑149 (MQ=255)
                                                                 gTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGa                                                                                    <  1:933514/149‑1 (MQ=255)
                                                                  ttGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGAt                                                                                   >  1:5020188/1‑149 (MQ=255)
                                                                    gagaCCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATgg                                                                                 >  1:278904/1‑149 (MQ=255)
                                                                      gaCCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTc                                                                               >  2:1646927/1‑149 (MQ=255)
                                                                      gaCCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTc                                                                               <  1:1264620/149‑1 (MQ=255)
                                                                        ccTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCTACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTg                                                                             >  1:2922671/1‑149 (MQ=255)
                                                                          tGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCTACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTc                                                                           <  1:2454476/149‑1 (MQ=255)
                                                                           gCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCTACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCg                                                                          <  2:3618709/149‑1 (MQ=255)
                                                                             aTCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTa                                                                        >  1:885807/1‑149 (MQ=255)
                                                                             aTCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTa                                                                        <  1:1025211/149‑1 (MQ=255)
                                                                              tCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTag                                                                       >  2:4451197/1‑149 (MQ=255)
                                                                                    agcagcTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGc                                                                 <  1:1762265/149‑1 (MQ=255)
                                                                                    agcagcTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGc                                                                 >  1:4631180/1‑149 (MQ=255)
                                                                                      cagcTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCtt                                                               >  2:3202193/1‑149 (MQ=255)
                                                                                      cagcTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCtt                                                               <  1:3850939/149‑1 (MQ=255)
                                                                                         cTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCttcat                                                            >  2:3783961/1‑149 (MQ=255)
                                                                                          ttttCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCg                                                                          <  2:5141996/134‑1 (MQ=255)
                                                                                          ttttCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCg                                                                          >  1:5141891/1‑134 (MQ=255)
                                                                                          ttttCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCttcatt                                                           <  2:2613624/149‑1 (MQ=255)
                                                                                           tttCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTTAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCttcattc                                                          <  2:1893942/149‑1 (MQ=255)
                                                                                           tttCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCttcattc                                                          <  1:153247/149‑1 (MQ=255)
                                                                                            ttCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCttcattca                                                         <  1:588213/149‑1 (MQ=255)
                                                                                            ttCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCttcattca                                                         <  1:4194942/149‑1 (MQ=255)
                                                                                               cAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAggc                                                      >  2:2003998/1‑149 (MQ=255)
                                                                                               cAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAggc                                                      <  1:25167/149‑1 (MQ=255)
                                                                                                aGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAggct                                                     >  2:3414783/1‑149 (MQ=255)
                                                                                                aGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAggct                                                     <  2:2549056/149‑1 (MQ=255)
                                                                                                aGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAggct                                                     <  2:3191201/149‑1 (MQ=255)
                                                                                                 ggCCAACGCGACGCAAGATCTCCTACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAggctg                                                    <  2:2072726/149‑1 (MQ=255)
                                                                                                    cAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAggctggct                                                 >  1:4419462/1‑149 (MQ=255)
                                                                                                     aaCGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTa                                                <  2:234900/149‑1 (MQ=255)
                                                                                                     aaCGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTa                                                <  1:3047340/149‑1 (MQ=255)
                                                                                                     aaCGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTa                                                >  2:2143913/1‑149 (MQ=255)
                                                                                                       cgcgACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCCCGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATg                                              <  2:1927914/149‑1 (MQ=255)
                                                                                                       cgcgACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATg                                              <  1:2715219/149‑1 (MQ=255)
                                                                                                         cgACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATggg                                            <  2:4035617/149‑1 (MQ=255)
                                                                                                           aCGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGc                                          <  2:2753129/149‑1 (MQ=255)
                                                                                                              cAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTAt                                       <  2:3717219/149‑1 (MQ=255)
                                                                                                                aGATCTCCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCGTCATTCAGGCTGGCGTTGGGGCGATCg                                     >  1:3569185/1‑149 (MQ=255)
                                                                                                                aGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCg                                     >  2:3876647/1‑149 (MQ=255)
                                                                                                                   tctcCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTa                                  >  1:107000/1‑149 (MQ=255)
                                                                                                                   tctcCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTa                                  <  2:1745529/149‑1 (MQ=255)
                                                                                                                   tctcCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTa                                  >  2:5081151/1‑149 (MQ=255)
                                                                                                                    ctcCTACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTaa                                 >  1:4535376/1‑149 (MQ=255)
                                                                                                                     tcCGACAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAg                                >  2:1294432/1‑149 (MQ=255)
                                                                                                                     tcCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAg                                >  2:3273997/1‑149 (MQ=255)
                                                                                                                          cAGAGCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCaaaa                           <  1:4761042/149‑1 (MQ=255)
                                                                                                                           agagCCTCATCACTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTAAGGCTGGCTATGGGGCTATAGTTAAGCAAAAt                          >  1:4872001/1‑149 (MQ=255)
                                                                                                                           agagCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAAt                          <  2:2215091/149‑1 (MQ=255)
                                                                                                                           agagCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAAt                          <  2:1075581/149‑1 (MQ=255)
                                                                                                                            gagCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAAtc                         <  1:4451169/149‑1 (MQ=255)
                                                                                                                            gagCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAAtc                         <  2:652089/149‑1 (MQ=255)
                                                                                                                              gCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAAtctc                       <  1:3127002/149‑1 (MQ=255)
                                                                                                                              gCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAAtctc                       <  1:1171415/149‑1 (MQ=255)
                                                                                                                              gCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAAtcta                       >  2:5191856/1‑148 (MQ=255)
                                                                                                                               ccTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCg                      >  1:1436563/1‑149 (MQ=255)
                                                                                                                                cTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGc                     <  2:4222322/149‑1 (MQ=255)
                                                                                                                                 tcatcaCTACTGCCAGGCTAGGCGAGTCAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGcc                    >  2:88338/1‑149 (MQ=255)
                                                                                                                                   atcaCTACTGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCGATGGGGCGATCGTTAAGCAAAATCTCGCCCt                  >  1:3298005/1‑149 (MQ=255)
                                                                                                                                    tcaCTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCGGGCTATTGGGCTATCGTTAAGCAAAATCACGCCCTg                 >  2:1463801/1‑149 (MQ=255)
                                                                                                                                     caCTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGCCCTGt                <  2:4835208/149‑1 (MQ=255)
                                                                                                                                        tactGCCAGGCTAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGCCCTGTTGc             <  1:1963499/149‑1 (MQ=255)
                                                                                                                                          ctGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGCCCTGTTGCgg           >  2:1573525/1‑149 (MQ=255)
                                                                                                                                            gCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGCCCTGTTGCGGGt         <  2:3587882/149‑1 (MQ=255)
                                                                                                                                                ggcgaggcgaGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGCCCTGTTGCGGGTCCCa     >  2:601004/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGGCGATAGCCAGACGGCGCAGTTCACCACCGGAGAGCTGGCGTCCGCCTTCACCTAACCAACTGTTGAGACCTGCATCCTCGAGCAGCTTTTCCAGGCCAACGCGACGCAAGATCTCCGACAGAGCCTCATCACTACTGCCAGGCGAGGCGAGTAAAAGATTATCACGCAGCGTGGCGCTAAACAGATGCACTCGCTGAGGAACAACGCTGATGGTCTGTCGTAGAGCCGCTTCATTCAGGCTGGCTATGGGGCTATCGTTAAGCAAAATCTCGCCCTGTTGCGGGTCCCATGC  >  CP009273/923192‑923487

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: