Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,965,573 C→T 8.0% intergenic (+71/+29) had → / ← PP_5232 (S)‑2‑haloacid dehalogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,965,5730CT 8.0% 177.1 / 4.1 62intergenic (+71/+29)had/PP_5232(S)‑2‑haloacid dehalogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base T (3/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTGGCGATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAATTCTTCCATCGCCTTGGCCAG  >  NC_002947/5965427‑5965721
                                                                                                                                                   |                                                                                                                                                    
gACTGGCGATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGa                                                                                                                                                    <  3:111972/149‑1 (MQ=255)
   tGGCGATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGAcgg                                                                                                                                                 >  3:122627/1‑149 (MQ=255)
     gCGATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGAcggcg                                                                                                                                               >  8:160636/1‑149 (MQ=255)
       gATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGGTCATAAAAAAGCCCGCAGCGAcggcgg                                                                                                                                              <  1:3528/149‑1 (MQ=255)
       gATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCggg                                                                                                                                             <  7:160636/149‑1 (MQ=255)
       gATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCAGAAGCGACGGCggg                                                                                                                                             >  3:25368/1‑149 (MQ=255)
                  ggATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGa                                                                                                                                  >  2:157869/1‑149 (MQ=255)
                        cGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCTGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGc                                                                                                                             >  2:69522/1‑148 (MQ=255)
                        cGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGc                                                                                                                             >  1:291149/1‑148 (MQ=255)
                         gAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCa                                                                                                                           <  3:233605/149‑1 (MQ=255)
                               cAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCTGGCACTGGCAACACAAATCCCACAGG‑CATAAAAAAGCCCGCAGCGATGGCGGGCTTTTTTGCGATCAGCCACTCAGa                                                                                                                     <  8:201088/149‑1 (MQ=255)
                                aGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGAt                                                                                                                    <  8:197464/149‑1 (MQ=255)
                                 gcgcCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATa                                                                                                                   <  5:14569/149‑1 (MQ=255)
                                     cTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGggc                                                                                                               >  3:71363/1‑149 (MQ=255)
                                         cGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCt                                                                                                           <  4:122627/149‑1 (MQ=255)
                                            aGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGcc                                                                                                        <  6:196008/149‑1 (MQ=255)
                                                   cGACGTGCTCGCGCGGTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGGGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACttg                                                                                                 <  7:10654/149‑1 (MQ=255)
                                                   cGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACttg                                                                                                 >  5:23094/1‑149 (MQ=255)
                                                   cGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACttg                                                                                                 <  2:155715/149‑1 (MQ=255)
                                                         gCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCCGTCGCTGCGGGCtttt                                                                                                                                        >  8:221683/1‑104 (MQ=38)
                                                         gCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCCGTCGCTGCGGGCtttt                                                                                                                                        <  7:221683/104‑1 (MQ=38)
                                                             gcgcgcTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGc                                                                                        <  4:129807/148‑1 (MQ=255)
                                                                 gctggcgctgACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGcttctt                                                                                   >  5:229211/1‑149 (MQ=255)
                                                                  ctggcgctgACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCAGCCGTACTTGTTGTCTGGCTTCTTg                                                                                  >  3:82576/1‑149 (MQ=255)
                                                                        ctgACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTc                                                                                            <  4:97473/133‑1 (MQ=255)
                                                                        ctgACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTc                                                                                            >  3:97473/1‑133 (MQ=255)
                                                                        ctgACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCggg                                                                            <  1:157869/149‑1 (MQ=255)
                                                                        ctgACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCggg                                                                            >  2:26830/1‑149 (MQ=255)
                                                                                 gCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGAcc                                                                   <  6:23094/149‑1 (MQ=255)
                                                                                   cGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCAc                                                                 >  6:122751/1‑149 (MQ=255)
                                                                                     cTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGTGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGt                                                               >  6:116140/1‑149 (MQ=255)
                                                                                       ttGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTg                                                             >  8:22337/1‑149 (MQ=255)
                                                                                        tgctgcCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTg                                                             <  1:69522/148‑1 (MQ=255)
                                                                                        tGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGtt                                                              <  1:25301/147‑1 (MQ=255)
                                                                                            gCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGTGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTc                                                        <  5:116140/149‑1 (MQ=255)
                                                                                             cccATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTc                                                        >  8:147591/1‑148 (MQ=255)
                                                                                              ccATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCg                                                       >  3:295400/1‑148 (MQ=255)
                                                                                                   cgcgGGCACGGGCAACACAAAACCCACATG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTc                                                  <  6:229211/148‑1 (MQ=255)
                                                                                                      gggcacgggcaACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGc                                              <  5:122751/149‑1 (MQ=255)
                                                                                                      gggcacgggcaACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGc                                              <  5:79938/149‑1 (MQ=255)
                                                                                                                aaCACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTgccgc                                    >  4:245856/1‑149 (MQ=255)
                                                                                                                   acaAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCcgccgc                                 >  2:87328/1‑149 (MQ=255)
                                                                                                                    caAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCcgccgc                                 >  7:28988/1‑148 (MQ=255)
                                                                                                                     aaaaCCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCga                               <  7:147591/149‑1 (MQ=255)
                                                                                                                         cccACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTTCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGaaa                           <  4:295400/149‑1 (MQ=255)
                                                                                                                         cccACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGaaa                           <  3:245856/149‑1 (MQ=255)
                                                                                                                         cccACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGaaa                           <  7:22337/149‑1 (MQ=255)
                                                                                                                              aGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCtt                                         >  3:106089/1‑130 (MQ=255)
                                                                                                                              aGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCtt                                         <  4:106089/130‑1 (MQ=255)
                                                                                                                              aGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGaaa                           >  4:620/1‑144 (MQ=255)
                                                                                                                              aGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGaaa                           <  3:620/144‑1 (MQ=255)
                                                                                                                               gg‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCCGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAAtt                     <  4:82576/149‑1 (MQ=255)
                                                                                                                                g‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAAttc                    >  5:320875/1‑149 (MQ=255)
                                                                                                                                   aTAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAAttctt                  >  1:244744/1‑149 (MQ=255)
                                                                                                                                   aTAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAAttctt                  <  2:244744/149‑1 (MQ=255)
                                                                                                                                      aaaaaGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAATTCTTcc                <  6:320875/148‑1 (MQ=255)
                                                                                                                                       aaaaGCCCGCCGTCGCTGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAAttcttc                 >  2:70102/1‑146 (MQ=255)
                                                                                                                                       aaaaGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAATTCTTCCAt              >  2:121718/1‑149 (MQ=255)
                                                                                                                                             ccGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCcgccgc                                 >  7:29644/1‑124 (MQ=255)
                                                                                                                                             ccGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCcgccgc                                 <  8:29644/124‑1 (MQ=255)
                                                                                                                                             ccGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAATTCTTCCATCGCCtt        >  3:97189/1‑149 (MQ=255)
                                                                                                                                                  gCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAATTCTTCCATCGCCTTGGCCa   <  8:28988/149‑1 (MQ=255)
                                                                                                                                                   cGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAATTCTTCCATCGCCTTGGCCAg  >  5:83101/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACTGGCGATCAGGCGCCGGATGCCGAAATCCAGCGCCTGTCGCAGCTGCCCGACGTGCTCGCGCGCTGGCGCTGACGGGAGCCGCTTTGCGGCCCATTCGCGGGCACGGGCAACACAAAACCCACAGG‑CATAAAAAAGCCCGCAGCGACGGCGGGCTTTTTTGCGATCAGCCACTCAGATAGGGCGGCTGCCGTACTTGTTGTCTGGCTTCTTGGGCGGGTCAGCGACCACGTTGGCCTCGACTTCCTGCACCTTGCCGCCGCGAGAAAGGAATTCTTCCATCGCCTTGGCCAG  >  NC_002947/5965427‑5965721

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