Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,768,058 G→C 5.6% P252A (CCT→GCT) ‡ 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,0580GC 5.6% 288.5 / 6.4 89P252A (CCT→GCT) ‡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 = 7.40e-01

ATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGC  >  NC_002947/4767913‑4768206
                                                                                                                                                 |                                                                                                                                                      
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                                                                      >  3:245360/1‑149 (MQ=255)
                                                                             cAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                                      <  8:394134/149‑1 (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                                                                      <  8:30671/126‑1 (MQ=255)
                                                                                                    ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                                      >  7:30671/1‑126 (MQ=255)
                                                                                                    ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                              <  3:46723/134‑1 (MQ=255)
                                                                                                    ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                              >  4:46723/1‑134 (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                                     >  3:162886/1‑131 (MQ=255)
                                                                                                                                aaTCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCa                                     <  4:162886/131‑1 (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)
                                                                                                                                                 g‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGgccgc  <  4:245360/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
ATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGC  >  NC_002947/4767913‑4768206

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: