Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 551,499 G→T 100% R171S (CGC→AGC)  ybcI ← DUF457 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273551,4990GT99.5% 626.6 / ‑4.3 191R171S (CGC→AGC) ybcIDUF457 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (104/86);  total (104/87)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 191 total aligned reads displayed.

AAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCACCTGCCAGGGAGCGAAAAAGC  >  CP009273/551355‑551642
                                                                                                                                                |                                                                                                                                               
 aGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTccac                                                                                                                                             >  2:1649345/1‑148 (MQ=255)
 aGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTccac                                                                                                                                             >  2:378170/1‑148 (MQ=255)
  gACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTccacc                                                                                                                                            >  2:3126203/1‑148 (MQ=255)
     tGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACcaca                                                                                                                                         >  1:1185107/1‑148 (MQ=255)
     tGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACcaca                                                                                                                                         >  1:2634791/1‑148 (MQ=255)
       ggcggcggcGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCacaac                                                                                                                                       >  1:2484551/1‑148 (MQ=255)
       ggcggcggcGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCacaac                                                                                                                                       >  1:1547862/1‑148 (MQ=255)
         cggcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACAt                                                                                                                                     <  2:3167835/148‑1 (MQ=255)
         cggcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACAt                                                                                                                                     <  2:2271950/148‑1 (MQ=255)
         cggcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACAt                                                                                                                                     >  1:810796/1‑148 (MQ=255)
           gcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTc                                                                                                                                   <  1:333116/148‑1 (MQ=255)
           gcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTc                                                                                                                                   <  1:203508/148‑1 (MQ=255)
            cggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTcc                                                                                                                                  <  2:2681324/148‑1 (MQ=255)
              gcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCa                                                                                                                                <  2:1163686/148‑1 (MQ=255)
               cggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCAt                                                                                                                               <  1:3385817/148‑1 (MQ=255)
               cggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCAt                                                                                                                               <  1:822967/148‑1 (MQ=255)
               cggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCAt                                                                                                                               <  2:90896/148‑1 (MQ=255)
                ggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATc                                                                                                                              <  1:2680456/148‑1 (MQ=255)
                        cgcgTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCAt                                                                                                                      <  2:240449/148‑1 (MQ=255)
                         gcgTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATc                                                                                                                     <  1:762604/148‑1 (MQ=255)
                            tGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATcccc                                                                                                                  <  2:2118762/148‑1 (MQ=255)
                                  tCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAc                                                                                                            >  1:1219808/1‑148 (MQ=255)
                                  tCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAc                                                                                                            >  2:3305714/1‑148 (MQ=255)
                                   cTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCAGCAGTAGCATCCCCGGCAAcc                                                                                                           >  2:1022503/1‑148 (MQ=255)
                                    ttAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACca                                                                                                          <  1:860135/148‑1 (MQ=255)
                                     tAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACcac                                                                                                         >  1:529337/1‑148 (MQ=255)
                                     tAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACcac                                                                                                         >  1:3124436/1‑148 (MQ=255)
                                       aCGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCacac                                                                                                       >  1:1985875/1‑148 (MQ=255)
                                       aCGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCacac                                                                                                       >  2:460927/1‑148 (MQ=255)
                                        cGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACAcc                                                                                                      >  2:2816868/1‑148 (MQ=255)
                                         gagaTGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACAccc                                                                                                     <  1:2268585/148‑1 (MQ=255)
                                               gggATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCa                                                                                               >  2:2507849/1‑148 (MQ=255)
                                                ggATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCaa                                                                                              >  1:1493800/1‑148 (MQ=255)
                                                 gATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAAt                                                                                             >  2:1399285/1‑148 (MQ=255)
                                                          ggAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaata                                                                                    <  1:768211/148‑1 (MQ=255)
                                                          ggAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACAGCAATTCAGAaata                                                                                    >  2:755450/1‑148 (MQ=255)
                                                            aaGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaataat                                                                                  >  2:2003165/1‑148 (MQ=255)
                                                            aaGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaataat                                                                                  >  1:2603898/1‑148 (MQ=255)
                                                             aGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATc                                                                                 <  2:1479850/148‑1 (MQ=255)
                                                             aGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATc                                                                                 >  1:278466/1‑148 (MQ=255)
                                                              gATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCTCTATTGCCTGATGCGAGTTTTCTCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCAGGCAACCACACCCACATCAATTCAGAAATAATCa                                                                                <  2:2933478/148‑1 (MQ=255)
                                                               aTGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCAc                                                                               >  2:1449804/1‑148 (MQ=255)
                                                                  gCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGACATAATCACCTg                                                                            >  2:1275194/1‑148 (MQ=255)
                                                                  gCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTg                                                                            >  2:916183/1‑148 (MQ=255)
                                                                   cccGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGa                                                                           >  2:522140/1‑148 (MQ=255)
                                                                       cAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGcc                                                                       >  1:103025/1‑148 (MQ=255)
                                                                       cAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGcc                                                                       >  1:103091/1‑148 (MQ=255)
                                                                       cAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGcc                                                                       <  2:2603898/148‑1 (MQ=255)
                                                                         aGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCg                                                                     >  2:1407670/1‑148 (MQ=255)
                                                                              ccaccTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACgg                                                                >  1:1568408/1‑148 (MQ=255)
                                                                                accTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCg                                                              >  1:3234642/1‑148 (MQ=255)
                                                                                accTGGCGTCGTAAGTAATTGCGCTAGTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCg                                                              <  2:3299338/148‑1 (MQ=255)
                                                                                   tGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtgg                                                           <  1:1290548/148‑1 (MQ=255)
                                                                                    ggCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggt                                                          <  1:544674/148‑1 (MQ=255)
                                                                                     gCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggtg                                                         <  1:1896037/148‑1 (MQ=255)
                                                                                          gtaagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACg                                                    <  2:174735/148‑1 (MQ=255)
                                                                                           taagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACga                                                   >  2:819754/1‑148 (MQ=255)
                                                                                            aagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgag                                                  >  1:2162253/1‑148 (MQ=255)
                                                                                            aagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgag                                                  >  1:785237/1‑148 (MQ=255)
                                                                                             agtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                                                 <  1:1645491/148‑1 (MQ=255)
                                                                                             agtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                                                 <  2:2021737/148‑1 (MQ=255)
                                                                                               taaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTCACGAGACa                                               >  2:3093708/1‑148 (MQ=255)
                                                                                                  ttGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATg                                            >  1:304028/1‑148 (MQ=255)
                                                                                                     cgcTATTGCCGGATGCGAGTTTTCGCATCCGGTTTTCGTCTGCTCCACCACAACCTTCCCATCAGTAGCATCCCCGGCAACCACACCCTCATCAATTCAGAAATAATCACCTGATGCCCGTACGGCTTGGTGTAACGATACAATGcta                                         <  2:1189991/148‑3 (MQ=255)
                                                                                                     cgcTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGcta                                         <  2:1189874/148‑3 (MQ=255)
                                                                                                         aTTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACgg                                     >  1:703634/1‑148 (MQ=255)
                                                                                                         aTTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACgg                                     <  1:3397895/148‑1 (MQ=255)
                                                                                                           tGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGgcg                                   >  2:1239526/1‑148 (MQ=255)
                                                                                                              cGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGa                                <  1:460927/148‑1 (MQ=255)
                                                                                                                 aTGCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACtt                             <  2:29843/148‑1 (MQ=255)
                                                                                                                   gCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACtttt                           <  1:1467600/148‑1 (MQ=255)
                                                                                                                   gCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACtttt                           <  2:1219808/148‑1 (MQ=255)
                                                                                                                   gCGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACtttt                           >  2:2512719/1‑148 (MQ=255)
                                                                                                                    cGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTa                          >  2:662866/1‑148 (MQ=255)
                                                                                                                    cGAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTa                          <  1:3126203/148‑1 (MQ=255)
                                                                                                                     gAGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTAt                         <  1:3354131/148‑1 (MQ=255)
                                                                                                                      aGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCATTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGGGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATc                        >  1:2883982/1‑148 (MQ=255)
                                                                                                                      aGTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGACATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATc                        >  2:3088722/1‑148 (MQ=255)
                                                                                                                       gTTTTCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCa                       <  1:192898/148‑1 (MQ=255)
                                                                                                                           tCGCATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTg                   >  1:2188736/1‑148 (MQ=255)
                                                                                                                              cATCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCa                <  1:675872/148‑1 (MQ=255)
                                                                                                                               aTCCGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAg               >  2:1690412/1‑148 (MQ=255)
                                                                                                                                 ccGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAggg             >  1:1690367/1‑148 (MQ=255)
                                                                                                                                 ccGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAggg             >  1:1110474/1‑148 (MQ=255)
                                                                                                                                 ccGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAggg             <  1:284253/148‑1 (MQ=255)
                                                                                                                                 ccGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAggg             <  1:284251/148‑1 (MQ=255)
                                                                                                                                 ccGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAggg             >  1:1689546/1‑148 (MQ=255)
                                                                                                                                 ccGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAggg             <  1:2086644/148‑1 (MQ=255)
                                                                                                                                 ccGGTTATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATAAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAggg             <  1:608083/148‑1 (MQ=255)
                                                                                                                                     ttATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCg         >  2:2168929/1‑148 (MQ=255)
                                                                                                                                     ttATCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCg         >  1:2536819/1‑148 (MQ=255)
                                                                                                                                       aTCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGCCTTGTATCACCTGCCAGGGAGCGaa       >  2:2564278/1‑148 (MQ=255)
                                                                                                                                       aTCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCGaa       >  2:221761/1‑148 (MQ=255)
                                                                                                                                       aTCGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCGaa       >  2:806098/1‑148 (MQ=255)
                                                                                                                                         cGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCGaaaa     >  2:1887687/1‑148 (MQ=255)
                                                                                                                                         cGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCGaaaa     <  2:3007735/148‑1 (MQ=255)
                                                                                                                                         cGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCGaaaa     <  2:1126986/148‑1 (MQ=255)
                                                                                                                                         cGTCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCGaaaa     <  1:819754/148‑1 (MQ=255)
                                                                                                                                           tCTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTTTATCACCTGCCAGGGAGCGAAAAAg   <  2:1060922/148‑1 (MQ=255)
                                                                                                                                            cTGCTCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCTAACGGCGCGACTTGTATCACCTGCCAGGGAGCGAAAAAGc  >  2:2796965/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCACCTGCCAGGGAGCGAAAAAGC  >  CP009273/551355‑551642

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: