Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,549,180:1 +C 5.5% intergenic (+52/‑84) pydB → / → pydX bifunctional D‑hydantoinase/dihydropyrimidinase/NADP‑dependent dihydropyrimidine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,549,1801.C 5.5% 301.2 / 10.5 73intergenic (+52/‑84)pydB/pydXbifunctional D‑hydantoinase/dihydropyrimidinase/NADP‑dependent dihydropyrimidine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base . (37/32);  new base C (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

CGCGGGGCGCTATGTGGAACGGCCGGCGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGC  >  NC_002947/4549032‑4549326
                                                                                                                                                      |                                                                                                                                                   
cgcgGGGCGCTATGTGGAACGGCCGGCGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAg                                                                                                                                                      <  7:80852/149‑1 (MQ=255)
    gggCGCTATGTGGAACGGCCGGCGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCATGCCCATCGCGACGCAAGGCCGCTCCCACAAGGAc                                                                                                                                                  <  1:53678/149‑1 (MQ=255)
            tgtgGAACGGCCGGCGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATc                                                                                                                                           <  5:83420/147‑1 (MQ=255)
              tgGAACGGCCGGCGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCt                                                                                                                                       <  6:275376/149‑1 (MQ=255)
                  aCGGCCGGCGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTgg                                                                                                                                   <  8:178837/149‑1 (MQ=255)
                    ggccggcGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTgg                                                                                                                                   >  5:124899/1‑147 (MQ=255)
                    ggccggcGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGa                                                                                                                                 <  7:109419/149‑1 (MQ=255)
                      ccggcGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGGCGATCGCGGTCCTTGTGGGAGc                                                                                                                               <  5:324684/149‑1 (MQ=38)
                      ccggcGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGc                                                                                                                               <  5:100404/149‑1 (MQ=255)
                       cggcGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGAc                                                                                                                                                  >  8:244496/1‑129 (MQ=255)
                       cggcGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGAc                                                                                                                                                  <  7:244496/129‑1 (MQ=255)
                        ggcGTACCCGGTGGTGTATGAGGTGGTGGGGCGCCGGGCCGAGCATCAGCGGGCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCGACAAGGACCGCGTTCGCCTGTGGGAGCgg                                                                                                                             <  2:324265/149‑1 (MQ=39)
                        ggcGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCgg                                                                                                                             <  4:283683/149‑1 (MQ=255)
                            tACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGGCGATCGCGGTCCTTGTGGGAGCGGcctt                                                                                                                         <  7:395/149‑3 (MQ=38)
                             aCCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCg                                                                                                                              >  2:169720/1‑143 (MQ=255)
                             aCCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCg                                                                                                                              <  1:169720/143‑1 (MQ=255)
                                cggcggTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGaa                                                                                                                     <  7:72844/149‑1 (MQ=255)
                                  gcggTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATc                                                                                                                   <  3:365328/149‑1 (MQ=255)
                                     gtgtATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGAc                                                                                                                <  6:26708/149‑1 (MQ=255)
                                         aTGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGc                                                                                                             <  2:175125/148‑1 (MQ=255)
                                          tGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCaa                                                                                                           >  3:83847/1‑149 (MQ=255)
                                           gAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCaa                                                                                                           >  8:161398/1‑148 (MQ=255)
                                             ggTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATAc                                                                                                        <  4:49495/149‑1 (MQ=255)
                                              gTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATAcc                                                                                                       >  4:389224/1‑149 (MQ=255)
                                              gTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATAcc                                                                                                       >  6:15556/1‑149 (MQ=255)
                                              gTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATCCCTTATACGAGCTGAAGTCATTGGGGCTGATGTGAAGCCCATCGCGACGCAAGCCCGCTCCAAAAAAGAACGAGATCCCCTGTGGGAGCGGCCCCAGGATCGACCGCAATAcc                                                                                                       >  4:238079/1‑149 (MQ=255)
                                                    gggCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGt                                                                                                    >  2:252224/1‑146 (MQ=255)
                                                    gggCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGt                                                                                                    <  1:252224/146‑1 (MQ=255)
                                                        gcCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGAGCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGc                                                                                             <  3:196954/149‑1 (MQ=255)
                                                          cGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCAc                                                                                           >  2:64813/1‑149 (MQ=255)
                                                          cGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCAc                                                                                           >  2:59128/1‑149 (MQ=255)
                                                                   cATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCaaaaaaaa                                                                                 <  5:15556/149‑1 (MQ=255)
                                                                       aGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATa                                                                               >  6:325036/1‑147 (MQ=255)
                                                                         cGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACGGCAATACCTGTCGCCTGCACCAAAAAAAAATagag                                                                            <  1:64813/149‑1 (MQ=255)
                                                                         cGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAAATagag                                                                            <  1:59128/149‑1 (MQ=255)
                                                                                   cTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAAcc                                                                 >  3:204426/1‑149 (MQ=255)
                                                                                         cgcgCTGAAGTCATTGTGGCTGCTGTGCAGCCCATCGCTACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGAc                                                            <  3:238079/148‑1 (MQ=255)
                                                                                          gcgcTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGAcc                                                           >  1:195594/1‑148 (MQ=255)
                                                                                           cgcTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTc                                                                                                   <  3:161515/108‑1 (MQ=255)
                                                                                           cgcTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTc                                                                                                   >  4:161515/1‑108 (MQ=255)
                                                                                                aaGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGccc                                                    >  1:327488/1‑149 (MQ=255)
                                                                                                aaGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTCTGGCAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATACAGCGGCTACAACCGTGCCCAACGccc                                                    >  2:344972/1‑149 (MQ=39)
                                                                                                     aTTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAAc                                               >  1:388835/1‑149 (MQ=255)
                                                                                                       atgggCTGCACAGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCCACAACCGTGACCAACGccc                                                    >  2:88036/3‑142 (MQ=255)
                                                                                                       atgggCTGCACAGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCCACAACCGTGACCAACGccc                                                    <  1:88036/140‑1 (MQ=255)
                                                                                                         gggCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCAtt                                           >  3:28424/1‑149 (MQ=255)
                                                                                                         gggCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCAtt                                           >  6:305621/1‑149 (MQ=255)
                                                                                                         gggCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCAtt                                           <  2:211331/149‑1 (MQ=255)
                                                                                                             tgctgTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCaaaaaaa                                                                                  >  6:24317/1‑106 (MQ=255)
                                                                                                             tgctgTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCaaaaaaa                                                                                  <  5:24317/106‑1 (MQ=255)
                                                                                                             tgctgTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGcc                                       >  8:327601/1‑149 (MQ=255)
                                                                                                                   gCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGccc                                 >  5:38588/1‑149 (MQ=255)
                                                                                                                    cAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGcccc                                >  5:228990/1‑149 (MQ=255)
                                                                                                                     aGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGccccc                               >  8:19603/1‑149 (MQ=255)
                                                                                                                        cccATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCgg                             <  3:389224/149‑1 (MQ=255)
                                                                                                                           aTCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGcc                          >  1:276661/1‑149 (MQ=255)
                                                                                                                            tCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCg                         >  5:302958/1‑149 (MQ=255)
                                                                                                                            tCGCGACGCAAGGCCGCTCCCACAAGGACCGCGAGCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTCTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCg                         <  2:271375/149‑1 (MQ=255)
                                                                                                                              gcgACGCAATGCCGCTCCCACAAGGCCCGCGATCGCCTGTGGGAGCGGCCCCAGCATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGgc                       <  8:358600/149‑1 (MQ=255)
                                                                                                                                    tggggCCGCTCCCACAGGCGATCGCGGTCCTTGTAGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGCCCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGAc                 >  1:52129/4‑149 (MQ=37)
                                                                                                                                     ggggCCGCTCCCACAGGCGATCGCGGTCCTTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGCCCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGAcc                >  7:324262/3‑149 (MQ=37)
                                                                                                                                     aaGGCCGCTCCCACAAGGACCGCTATCGCGTGTGGGAGCGGCCCCAGAAGCGACCGCAATACCTGTCGCCGGCACCAAAAAAAATAGAGAGGCTACCACCGTGACCACCGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGa                  >  7:60653/1‑147 (MQ=255)
                                                                                                                                     aaGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCACAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCACCGCCCGGAACCATTTGACGCGCCCCCGGGCCGGCGCCGAcc                >  7:339623/1‑149 (MQ=255)
                                                                                                                                      aGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCa               <  6:38588/149‑1 (MQ=255)
                                                                                                                                        gCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGt             >  5:82493/1‑149 (MQ=255)
                                                                                                                                         ccGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGtt            <  6:82493/149‑1 (MQ=255)
                                                                                                                                          cGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCAGCGCCCCCGGGCCCGCGCCGACCAGtt            >  2:158661/1‑149 (MQ=255)
                                                                                                                                           gCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTgg          >  2:211905/1‑149 (MQ=255)
                                                                                                                                               ccACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGcc        >  6:261383/1‑147 (MQ=255)
                                                                                                                                                 acaAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAgc    >  7:261882/1‑149 (MQ=255)
                                                                                                                                                 acaAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAgc    <  7:327601/149‑1 (MQ=255)
                                                                                                                                                 acaAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAgc    >  2:118378/1‑149 (MQ=255)
                                                                                                                                                 acaAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAgc    <  5:325036/149‑1 (MQ=255)
                                                                                                                                                   aaGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAAAGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGcgcagctc  >  7:151073/1‑142 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGCGGGGCGCTATGTGGAACGGCCGGCGTACCCGGCGGTGTATGAGGTGCTGGGGCGCCGGGCCGAGCATCAGCGGCCGATGCCTGTTCCGCGCTGAAGTCATTGGGGCTGCTGTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAAGGACCGCGATCGCCTGTGGGAGCGGCCCCAGAATCGACCGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGC  >  NC_002947/4549032‑4549326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: