Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,768,057 G→C 5.6% P252R (CCT→CGT) ‡ PP_4218 ← lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,768,0570GC 5.6% 285.6 / 6.5 89P252R (CCT→CGT) ‡PP_4218lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base G (46/38);  new base C (3/2);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGG  >  NC_002947/4767909‑4768202
                                                                                                                                                    |                                                                                                                                                   
ctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAg                                                                                                                                                     <  8:139792/149‑1 (MQ=255)
    atgatCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGc                                                                                                                                               <  3:306091/149‑1 (MQ=255)
    atgatCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGc                                                                                                                                               <  2:54339/149‑1 (MQ=255)
           aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGTCAATCAGCATCAGTGCAGG‑‑TGCACTGTTg                                                                                                                                        <  1:32258/149‑1 (MQ=255)
               cATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTgg                                                                                                                                    >  6:142718/1‑149 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTggg                                                                                                                                   >  3:201558/1‑148 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACC‑‑TGCACTGATGCTgat                                                                                                                                   <  2:257738/149‑3 (MQ=255)
                  aCGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGa                                                                                                                                 >  2:425902/1‑149 (MQ=255)
                       tGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTg                                                                                                                            <  2:76471/149‑1 (MQ=255)
                         aaGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGcc                                                                                                                          >  8:315487/1‑149 (MQ=255)
                           gTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCa                                                                                                                        <  7:315487/149‑1 (MQ=255)
                              aTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGc                                                                                                                      >  5:18376/1‑148 (MQ=255)
                               tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                                    >  7:394134/1‑149 (MQ=255)
                                gggTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAAc                                                                                                                   >  7:338451/1‑149 (MQ=255)
                                 ggTCACTCCCGGGATCACCGCAAGCTGCCCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACg                                                                                                                  <  5:254000/149‑1 (MQ=255)
                                  gTCACTCCCGGGATCACCGCAAGCTGGAGCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACgg                                                                                                                 <  7:17450/149‑1 (MQ=255)
                                  gTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACgg                                                                                                                 >  1:9691/1‑149 (MQ=255)
                                      cTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGa                                                                                                              >  2:301431/1‑148 (MQ=255)
                                       tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaa                                                                                                             >  2:238414/1‑148 (MQ=255)
                                          cGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCCCTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                         >  5:298871/1‑149 (MQ=255)
                                            ggATCACCGCAAGCTGCACCGCCACCCCTGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgc                                                                                                       <  8:259033/149‑1 (MQ=255)
                                               tCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGa                                                                                                    <  6:120518/149‑1 (MQ=255)
                                                 aCCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                                  >  7:357644/1‑149 (MQ=255)
                                                  ccGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCg                                                                                                 >  8:158404/1‑149 (MQ=255)
                                                   cGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                                >  5:206260/1‑149 (MQ=255)
                                                     cAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                                    <  7:371867/127‑1 (MQ=255)
                                                     cAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                                    >  8:371867/1‑127 (MQ=255)
                                                        gCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACgg                                                                                           >  8:115904/1‑149 (MQ=255)
                                                          tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTc                                                                                         >  8:325240/1‑149 (MQ=255)
                                                            cACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACg                                                                                                                  <  5:429814/122‑1 (MQ=255)
                                                            cACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACg                                                                                                                  >  6:429814/1‑122 (MQ=255)
                                                               cGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCg                                                                                    <  8:357644/149‑1 (MQ=255)
                                                                   aCCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCAc                                                                                >  4:55002/1‑149 (MQ=255)
                                                                     cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGc                                                                                                                                   <  6:268644/96‑1 (MQ=255)
                                                                     cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGc                                                                                                                                   >  5:268644/1‑96 (MQ=255)
                                                                         gCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGc                                                                          <  6:123726/149‑1 (MQ=255)
                                                                           ttGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGtt                                                                                                                                         >  6:11039/1‑84 (MQ=255)
                                                                           ttGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGtt                                                                                                                                         <  5:11039/84‑1 (MQ=255)
                                                                              cTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTAtt                                                                     <  1:425902/149‑1 (MQ=255)
                                                                                 cAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                                  <  8:394134/149‑1 (MQ=255)
                                                                                 cAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                                  >  3:245360/1‑149 (MQ=255)
                                                                                   aGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                <  3:55002/149‑1 (MQ=255)
                                                                                   aGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                <  2:9691/149‑1 (MQ=255)
                                                                                    gTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGt                                                               <  6:206260/149‑1 (MQ=255)
                                                                                       ccGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGGATGCAATGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAACCGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGt                                                               >  1:130753/1‑148 (MQ=255)
                                                                                           gCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCa                                                        >  1:181280/1‑149 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGGACTCGGCAATCAGCATCTGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCAGTCCAg                                                  <  6:442569/142‑1 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGTCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgc                                           >  2:98884/1‑149 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                                  >  7:30671/1‑126 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                                  <  8:30671/126‑1 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                          >  4:46723/1‑134 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                          <  3:46723/134‑1 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgc                                           >  5:442571/1‑149 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgc                                           >  3:265607/1‑149 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGCCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGATGGCGCCGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCCTGGCGTATTGCTGGTAGAACCACTGCAg                                                  >  5:442569/1‑142 (MQ=255)
                                                                                                             tCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTc                                      <  7:372442/149‑1 (MQ=255)
                                                                                                               aTGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTc                                      >  3:190280/1‑147 (MQ=255)
                                                                                                                tGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGATGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAg                                   >  5:279071/1‑149 (MQ=255)
                                                                                                                   aGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                <  8:338451/149‑1 (MQ=255)
                                                                                                                            cACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCt                       >  2:167800/1‑149 (MQ=255)
                                                                                                                            cACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCt                       >  4:52067/1‑149 (MQ=255)
                                                                                                                             aCTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTg                      >  1:293838/1‑149 (MQ=255)
                                                                                                                             aCTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTg                      >  7:2346/1‑149 (MQ=255)
                                                                                                                              cTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGc                     <  7:325240/149‑1 (MQ=255)
                                                                                                                                 ggCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCGCCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTACCa                   >  5:140776/1‑148 (MQ=255)
                                                                                                                                 ggCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGcc                    <  1:167800/147‑1 (MQ=255)
                                                                                                                                    aaTCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCa                                 <  4:162886/131‑1 (MQ=255)
                                                                                                                                    aaTCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCa                                 >  3:162886/1‑131 (MQ=255)
                                                                                                                                    aaTCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgc               <  6:18376/149‑1 (MQ=255)
                                                                                                                                    aaTCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgc               >  8:296079/1‑149 (MQ=255)
                                                                                                                                     aTCAGCATCAGTGCAGG‑‑TGCACTGTGGCTGGCGACGCTGCCCAGCAACGGTTAAAAACGCGGATCGTGACTGTCTGCCGGCACCGTGTCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgcg              <  1:365484/149‑1 (MQ=255)
                                                                                                                                       cAGCAACAGTGCACC‑‑TGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCt            >  1:403088/1‑149 (MQ=255)
                                                                                                                                        aGCATTAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                                 >  6:226149/1‑111 (MQ=255)
                                                                                                                                        aGCATTAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                                 <  5:226149/111‑1 (MQ=255)
                                                                                                                                         gCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTg           >  7:61887/1‑148 (MQ=255)
                                                                                                                                         gCAACAGTGCACC‑‑TGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTg           >  3:219290/1‑148 (MQ=255)
                                                                                                                                          cATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGGCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCa         <  1:301431/149‑1 (MQ=255)
                                                                                                                                           aacAGTGCACC‑‑TGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAc        >  5:289467/3‑149 (MQ=255)
                                                                                                                                           aTCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGATGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAc        <  6:279071/149‑1 (MQ=255)
                                                                                                                                           aTCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAc        >  4:405061/1‑149 (MQ=255)
                                                                                                                                            tCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAcc       >  1:352400/1‑149 (MQ=255)
                                                                                                                                              aGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCg      >  8:52779/1‑148 (MQ=255)
                                                                                                                                              aGTGCAGG‑‑TGAACTGTTGCTGGCTACGCTGCCCAGCAACGCTGAAAAACGACGATCGTGAACTTCTGCCGGCACCGTGGCGTATTCCTGTTAGAACCACTCCAGCGTCTGCGCTTCAAGCACGTAACCGCTGCCATACCGCTGCAcac      >  4:345317/1‑146 (MQ=255)
                                                                                                                                               gTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCATCAACGGTGAAAAACGCGGATGGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa    <  8:166948/149‑1 (MQ=255)
                                                                                                                                               gTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa    <  7:241749/149‑1 (MQ=255)
                                                                                                                                               gTGCACC‑‑TGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa    <  6:289467/149‑1 (MQ=255)
                                                                                                                                                tGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgc               <  3:66077/137‑1 (MQ=255)
                                                                                                                                                tGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgc               >  4:66077/1‑137 (MQ=255)
                                                                                                                                                tGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg   <  7:115904/149‑1 (MQ=255)
                                                                                                                                                  cAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAgg  <  3:345317/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGG  >  NC_002947/4767909‑4768202

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