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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,983,2980GT10.6% 607.9 / 33.5 178R277S (CGC→AGC) hemChydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base G (86/73);  new base T (13/6);  total (99/79)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 181 total aligned reads displayed.

TCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTA  >  CP009273/3983151‑3983443
                                                                                                                                                   |                                                                                                                                                 
tctcCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACcgcg                                                                                                                                                   >  1:2208667/1‑148 (MQ=255)
    cAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGt                                                                                                                                               <  1:1851012/148‑1 (MQ=255)
      gCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGCCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAAACCACATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTc                                                                                                                                             >  2:853048/1‑148 (MQ=255)
         ggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCAcc                                                                                                                                          >  1:2087060/1‑148 (MQ=255)
         ggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCAcc                                                                                                                                          >  2:761578/1‑148 (MQ=255)
          gagaCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcg                                                                                                                                         <  2:177368/148‑1 (MQ=255)
                ggcgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaata                                                                                                                                   <  2:613566/148‑1 (MQ=255)
                  cgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGaataat                                                                                                                                 <  1:728506/148‑1 (MQ=255)
                    ggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCt                                                                                                                               >  1:2679891/1‑148 (MQ=255)
                            ggagtcTCATGCCGGGGCGTCTCCGTTATAGTCTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCg                                                                                                                       <  2:301706/143‑1 (MQ=255)
                            ggATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCg                                                                                                                       >  1:873421/1‑148 (MQ=255)
                            ggATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCg                                                                                                                       >  1:3078392/1‑148 (MQ=255)
                             gATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGt                                                                                                                      >  2:2831021/1‑148 (MQ=255)
                              aTACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCCTTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTc                                                                                                                     >  1:2951521/1‑148 (MQ=255)
                              aTACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTc                                                                                                                     <  1:2446790/148‑1 (MQ=255)
                               tACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc                                                                                                                    >  1:762832/1‑148 (MQ=255)
                                aCTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCg                                                                                                                   >  1:460139/1‑148 (MQ=255)
                                  tCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgc                                                                                                                 >  2:726558/1‑148 (MQ=255)
                                  tCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgc                                                                                                                 <  2:24521/148‑1 (MQ=255)
                                          gggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                                         <  2:2876201/148‑1 (MQ=255)
                                          gggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                                         <  1:82919/148‑1 (MQ=255)
                                           ggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                                        <  2:758342/148‑1 (MQ=255)
                                           ggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                                        <  2:1426870/148‑1 (MQ=255)
                                           ggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                                        <  2:1427712/148‑1 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCCTCTTGCGGCGCACCGCGGCGTTCTCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc                                                                                                   >  2:1977537/1‑148 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc                                                                                                   >  1:2657682/1‑148 (MQ=255)
                                                  ccGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                 >  2:42994/1‑148 (MQ=255)
                                                   cGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacg                                                                                                <  1:1109466/148‑1 (MQ=255)
                                                   cGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacg                                                                                                <  1:1236166/148‑1 (MQ=255)
                                                         aGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                                          >  1:686748/1‑148 (MQ=255)
                                                           aCTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                                        <  1:1072108/148‑1 (MQ=255)
                                                           aCTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                                        >  1:1438835/1‑148 (MQ=255)
                                                            cTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAt                                                                                       <  2:767348/148‑1 (MQ=255)
                                                             ttCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAtt                                                                                      >  2:2666315/1‑148 (MQ=255)
                                                              tCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAttt                                                                                     >  1:290334/1‑148 (MQ=255)
                                                                aGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCg                                                                                   <  1:1867739/148‑1 (MQ=255)
                                                                  cGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGcc                                                                                 >  1:187435/1‑148 (MQ=255)
                                                                  cGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGcc                                                                                 >  2:1205042/1‑148 (MQ=255)
                                                                  cGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGcc                                                                                 >  2:1761616/1‑148 (MQ=255)
                                                                    aGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                                               <  2:362568/148‑1 (MQ=255)
                                                                    aGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                                               >  2:1258520/1‑148 (MQ=255)
                                                                    aGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                                               <  2:1075100/148‑1 (MQ=255)
                                                                      gATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCa                                                                             >  2:1535162/1‑148 (MQ=255)
                                                                      gATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCa                                                                             >  2:2716474/1‑148 (MQ=255)
                                                                       aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                                            <  2:3084627/148‑1 (MQ=255)
                                                                       aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                                            <  1:2423413/148‑1 (MQ=255)
                                                                       aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                                            <  1:763128/148‑1 (MQ=255)
                                                                       aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                                            <  1:1112106/148‑1 (MQ=255)
                                                                          tcGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATaa                                                                         >  2:1885870/1‑148 (MQ=255)
                                                                          tcGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATaa                                                                         <  2:1047621/148‑1 (MQ=255)
                                                                              gcgcgcCGTTATTCAGTAGCTCTTCTGCCAGCGATATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCTCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCATATTTCGCCATCAATAACCTc                                                                     <  2:684482/148‑1 (MQ=255)
                                                                                   ccGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                                                <  2:203315/148‑1 (MQ=255)
                                                                                   ccGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                                                >  1:3126875/1‑148 (MQ=255)
                                                                                     gTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGGCCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg                                                              <  1:686073/148‑1 (MQ=255)
                                                                                     gTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg                                                              <  1:2373546/148‑1 (MQ=255)
                                                                                     gTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg                                                              <  1:745202/148‑1 (MQ=255)
                                                                                       tagtCAGTAGCTCTTCTGCCAGCGAAATCCCCCGATGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                                            <  2:2550104/145‑1 (MQ=255)
                                                                                       tATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                                            >  1:229456/1‑148 (MQ=255)
                                                                                            aGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa                                                       >  1:1914161/1‑148 (MQ=255)
                                                                                             gTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAt                                                      <  2:763143/148‑1 (MQ=255)
                                                                                                gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTCCCAATTgg                                                   >  2:906325/1‑148 (MQ=255)
                                                                                                 ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                                                  >  1:928981/1‑148 (MQ=255)
                                                                                                 ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                                                  >  1:146267/1‑148 (MQ=255)
                                                                                                 ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGc                                                  >  2:2139529/1‑148 (MQ=255)
                                                                                                   cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                                                >  1:2014493/1‑148 (MQ=255)
                                                                                                      ctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                                             >  2:1106992/1‑148 (MQ=255)
                                                                                                           aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                                        >  1:774260/1‑148 (MQ=255)
                                                                                                             cGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCg                                      >  1:769138/1‑148 (MQ=255)
                                                                                                             cGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCg                                      >  1:2671677/1‑148 (MQ=255)
                                                                                                                 aTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt                                  >  2:380146/1‑148 (MQ=255)
                                                                                                                 aTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt                                  >  1:901776/1‑148 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                 >  2:2719562/1‑148 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                 >  1:2754498/1‑148 (MQ=255)
                                                                                                                   ccccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCg                                >  2:2321131/1‑148 (MQ=255)
                                                                                                                     ccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag                              <  2:1320018/148‑1 (MQ=255)
                                                                                                                      cATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga                             >  2:206445/1‑148 (MQ=255)
                                                                                                                      cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga                             <  1:1828139/148‑1 (MQ=255)
                                                                                                                       aTTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAc                            <  2:2581183/148‑1 (MQ=255)
                                                                                                                        tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg                           <  2:1456267/148‑1 (MQ=255)
                                                                                                                        tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg                           >  2:3184950/1‑148 (MQ=255)
                                                                                                                         ttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg                          >  1:2573405/1‑148 (MQ=255)
                                                                                                                         ttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg                          <  2:1634734/148‑1 (MQ=255)
                                                                                                                           gTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGt                        <  2:2087060/148‑1 (MQ=255)
                                                                                                                           gTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGt                        <  1:1065780/148‑1 (MQ=255)
                                                                                                                             tCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTAt                      <  2:2956283/148‑1 (MQ=255)
                                                                                                                               ggCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATGTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTc                    <  2:259373/148‑1 (MQ=255)
                                                                                                                                gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCa                   <  2:2924579/148‑1 (MQ=255)
                                                                                                                                gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCa                   <  1:3270291/148‑1 (MQ=255)
                                                                                                                                 cATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTCCCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCAt                  >  2:2896715/1‑148 (MQ=255)
                                                                                                                                 cATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCAt                  >  2:464307/1‑148 (MQ=255)
                                                                                                                                  aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATg                 <  2:1517581/148‑1 (MQ=255)
                                                                                                                                   tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATgg                >  1:2635181/1‑148 (MQ=255)
                                                                                                                                   tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATgg                >  1:577208/1‑148 (MQ=255)
                                                                                                                                   tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCGTCGAGCCGGGTATTCATgg                >  1:1733520/1‑148 (MQ=255)
                                                                                                                                      tGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGgcgc             >  1:2109043/1‑148 (MQ=255)
                                                                                                                                      tGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGgcgc             >  2:1142000/1‑148 (MQ=255)
                                                                                                                                        cggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGt           <  1:2444008/148‑1 (MQ=255)
                                                                                                                                         ggcgCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGtt          <  1:380146/148‑1 (MQ=255)
                                                                                                                                           cgcACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCt        >  1:1565672/1‑148 (MQ=255)
                                                                                                                                                 gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTa  <  1:2854688/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTA  >  CP009273/3983151‑3983443

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