Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,768,066 T→A 5.5% N249I (AAC→ATC)  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,0660TA 5.5% 270.9 / 3.1 89N249I (AAC→ATC) PP_4218lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base T (46/38);  new base A (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.92e-01

AACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACT  >  NC_002947/4767920‑4768214
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                >  8:371867/1‑127 (MQ=255)
                                          cAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                                                <  7:371867/127‑1 (MQ=255)
                                             gCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACgg                                                                                                       >  8:115904/1‑149 (MQ=255)
                                               tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTc                                                                                                     >  8:325240/1‑149 (MQ=255)
                                                 cACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACg                                                                                                                              >  6:429814/1‑122 (MQ=255)
                                                 cACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACg                                                                                                                              <  5:429814/122‑1 (MQ=255)
                                                    cGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCg                                                                                                <  8:357644/149‑1 (MQ=255)
                                                        aCCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCAc                                                                                            >  4:55002/1‑149 (MQ=255)
                                                          cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGc                                                                                                                                               >  5:268644/1‑96 (MQ=255)
                                                          cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGc                                                                                                                                               <  6:268644/96‑1 (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                                                                            <  2:9691/149‑1 (MQ=255)
                                                                        aGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                            <  3:55002/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                                             >  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)
                                                                                                                                             tGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGgccgcc         >  4:22086/1‑149 (MQ=255)
                                                                                                                                               cACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGccgccg        <  4:265607/148‑1 (MQ=255)
                                                                                                                                               cACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCgg       <  3:22086/149‑1 (MQ=255)
                                                                                                                                                    ttGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACt  >  7:79834/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGG‑‑TGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCCGGCAACT  >  NC_002947/4767920‑4768214

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: