Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,983,298 G→T 9.1% R277S (CGC→AGC)  hemC ← hydroxymethylbilane synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,983,2980GT 9.1% 697.6 / 30.4 196R277S (CGC→AGC) hemChydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base G (84/94);  new base T (8/10);  total (92/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 212 total aligned reads displayed.

TCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAA  >  CP009273/3983151‑3983444
                                                                                                                                                   |                                                                                                                                                  
tctcCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACcgcg                                                                                                                                                    >  1:788811/1‑148 (MQ=255)
   ccAGCGGGAGACGGGCGTGTGACAAGGATACTCATGCCTGTGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCgcggcg                                                                                                                                                 <  1:3093249/148‑1 (MQ=255)
   ccAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCgcggcg                                                                                                                                                 <  2:478260/148‑1 (MQ=255)
   ccAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCgcggcg                                                                                                                                                 >  2:3326781/1‑148 (MQ=255)
     aGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGtt                                                                                                                                               <  1:2226784/148‑1 (MQ=255)
           agaCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcgc                                                                                                                                         >  1:1049495/1‑148 (MQ=255)
             aCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGa                                                                                                                                       <  1:1627692/148‑1 (MQ=255)
                ggcgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaata                                                                                                                                    >  2:2835005/1‑148 (MQ=255)
                ggcgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaata                                                                                                                                    <  2:2445704/148‑1 (MQ=255)
                     gTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTg                                                                                                                               <  2:2775849/148‑1 (MQ=255)
                     gTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTg                                                                                                                               <  2:2584150/148‑1 (MQ=255)
                     gTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTGCACCGCGAATAATCTg                                                                                                                               >  1:450356/1‑148 (MQ=255)
                        aCAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGa                                                                                                                            <  2:488285/148‑1 (MQ=255)
                        aCAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGa                                                                                                                            >  1:2176263/1‑148 (MQ=255)
                        aCAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGa                                                                                                                            <  1:486540/148‑1 (MQ=255)
                         cAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGaa                                                                                                                           >  1:2816095/1‑148 (MQ=255)
                         cAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGaa                                                                                                                           >  1:98579/1‑148 (MQ=255)
                         cAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGaa                                                                                                                           <  2:3205335/148‑1 (MQ=255)
                         cAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGaa                                                                                                                           >  2:3156081/1‑148 (MQ=255)
                         cAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGaa                                                                                                                           >  2:1775168/1‑148 (MQ=255)
                          aaGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAAc                                                                                                                          <  1:2386750/148‑1 (MQ=255)
                          aaGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAAc                                                                                                                          <  1:1405650/148‑1 (MQ=255)
                            ggATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCg                                                                                                                        >  2:521564/1‑148 (MQ=255)
                            ggATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCg                                                                                                                        >  1:2510781/1‑148 (MQ=255)
                               tACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc                                                                                                                     <  1:1313170/148‑1 (MQ=255)
                               tACTCATGCCGGGGCGTCTCAGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc                                                                                                                     <  1:3076690/148‑1 (MQ=255)
                                aCTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCg                                                                                                                    <  2:2817616/148‑1 (MQ=255)
                                   cATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcg                                                                                                                 <  1:596345/148‑1 (MQ=255)
                                    aTGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcgc                                                                                                                >  1:3346824/1‑148 (MQ=255)
                                     tGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcgcg                                                                                                               >  1:1111441/1‑148 (MQ=255)
                                        cGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCg                                                                                                            <  1:3326501/148‑1 (MQ=255)
                                          gggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                                          <  1:1724511/148‑1 (MQ=255)
                                          gggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                                          <  2:910121/148‑1 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc                                                                                                    >  1:3538946/1‑148 (MQ=255)
                                                  ccGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                  >  1:3483249/1‑148 (MQ=255)
                                                  ccGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                  >  1:1162691/1‑148 (MQ=255)
                                                       ataGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                                             <  2:1664108/148‑1 (MQ=255)
                                                       ataGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                                             <  2:1331306/148‑1 (MQ=255)
                                                       ataGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                                             <  2:1925041/148‑1 (MQ=255)
                                                        taGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGcc                                                                                            <  1:1123010/148‑1 (MQ=255)
                                                         aGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                                           >  1:3174919/1‑148 (MQ=255)
                                                         aGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                                           >  2:803017/1‑148 (MQ=255)
                                                           aCTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                                         <  2:788851/148‑1 (MQ=255)
                                                            cTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAt                                                                                        <  2:1162691/148‑1 (MQ=255)
                                                            cTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAt                                                                                        <  2:1657117/148‑1 (MQ=255)
                                                            cTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAt                                                                                        >  2:3470896/1‑148 (MQ=255)
                                                              tCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAttt                                                                                      <  2:2128180/148‑1 (MQ=255)
                                                                  cGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCATATTTCGcc                                                                                  <  1:1209935/148‑1 (MQ=255)
                                                                  cGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGcc                                                                                  <  1:1209907/148‑1 (MQ=255)
                                                                    aGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                                                <  1:682070/148‑1 (MQ=255)
                                                                    aGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                                                <  2:100821/148‑1 (MQ=255)
                                                                    aGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                                                <  2:1901738/148‑1 (MQ=255)
                                                                     ggATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATc                                                                               >  1:450675/1‑148 (MQ=255)
                                                                     ggATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATc                                                                               >  2:1696392/1‑148 (MQ=255)
                                                                      gATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCa                                                                              <  2:1630502/148‑1 (MQ=255)
                                                                       aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCCTTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                                             <  2:3394589/148‑1 (MQ=255)
                                                                         cacgcccgcgCCGTTATTCAGTAGATCTTCTGCCAGCGAAATCCCCATTAGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACAGTCAGCCGCGCCGACAAGCGCACGCACCCAGATTTCGCCATCAATa                                                                           >  2:1943022/3‑148 (MQ=255)
                                                                           cgcgcgcgCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAg                                                                         >  1:1813129/1‑148 (MQ=255)
                                                                                    cGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                                                <  2:1087977/148‑1 (MQ=255)
                                                                                     gTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg                                                               >  1:2470189/1‑148 (MQ=255)
                                                                                     gTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg                                                               <  1:3530870/148‑1 (MQ=255)
                                                                                      ttATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                                              <  2:1428102/148‑1 (MQ=255)
                                                                                       tATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                                             >  1:2997949/1‑148 (MQ=255)
                                                                                          tCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAcc                                                          >  1:2165691/1‑148 (MQ=255)
                                                                                           cAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCTGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCa                                                         <  1:1659902/148‑1 (MQ=255)
                                                                                           cAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCa                                                         <  2:485690/148‑1 (MQ=255)
                                                                                                  tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCa                                                  >  2:1506178/1‑148 (MQ=255)
                                                                                                  tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCa                                                  >  2:2407304/1‑148 (MQ=255)
                                                                                                   cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                                                 >  1:3425350/1‑148 (MQ=255)
                                                                                                   cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                                                 >  1:2344789/1‑148 (MQ=255)
                                                                                                       tGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGa                                             <  1:1078454/148‑1 (MQ=255)
                                                                                                           aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCATCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                                         >  1:3286403/1‑148 (MQ=255)
                                                                                                              gAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGc                                      <  2:41543/148‑1 (MQ=255)
                                                                                                               aaaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGcc                                     <  2:711195/148‑1 (MQ=255)
                                                                                                                aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt                                    >  1:933620/1‑148 (MQ=255)
                                                                                                                aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt                                    <  1:204733/148‑1 (MQ=255)
                                                                                                                 aTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt                                   <  2:2644595/148‑1 (MQ=255)
                                                                                                                 aTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt                                   >  2:2330427/1‑148 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                  >  2:2706696/1‑148 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                  >  2:2706735/1‑148 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                  <  2:2815239/148‑1 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                  >  2:622242/1‑148 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                  <  2:3246827/148‑1 (MQ=255)
                                                                                                                   ccccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCTCCTGACATCCGCCTTCg                                 >  2:783734/1‑148 (MQ=255)
                                                                                                                      cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga                              <  1:2787296/148‑1 (MQ=255)
                                                                                                                        tttGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg                            <  2:642329/148‑1 (MQ=255)
                                                                                                                         ttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg                           >  2:2186077/1‑148 (MQ=255)
                                                                                                                         ttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg                           >  2:3011261/1‑148 (MQ=255)
                                                                                                                            ttCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTa                        >  1:2663270/1‑148 (MQ=255)
                                                                                                                                gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCa                    >  1:1947196/1‑148 (MQ=255)
                                                                                                                                  aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATg                  <  2:383829/148‑1 (MQ=255)
                                                                                                                                  aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATg                  >  1:1496977/1‑148 (MQ=255)
                                                                                                                                   tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATgg                 >  2:883881/1‑148 (MQ=255)
                                                                                                                                    cTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTCTTCATGgc                >  2:3321880/1‑148 (MQ=255)
                                                                                                                                     ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGgcg               <  1:2835005/148‑1 (MQ=255)
                                                                                                                                      tGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGgcgc              >  2:1892144/1‑148 (MQ=255)
                                                                                                                                         ggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGtt           <  2:1394988/148‑1 (MQ=255)
                                                                                                                                           cgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGGATTCATGGCGCGTTCt         >  1:1143451/1‑148 (MQ=255)
                                                                                                                                            gcACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTg        <  2:1562435/148‑1 (MQ=255)
                                                                                                                                                 gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGAGTTCTGCGGTa   <  1:3132195/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAA  >  CP009273/3983151‑3983444

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