Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 827,244 C→T 5.9% N144N (AAC→AAT ycaC‑I → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947827,2440CT 5.9% 1008.5 / 30.7 302N144N (AAC→AATycaC‑Iputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (130/154);  new base T (9/9);  total (139/163)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

GATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTA  >  NC_002947/827095‑827393
                                                                                                                                                     |                                                                                                                                                     
gATTTCGTCTAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTAGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                                                                                                                                       <  2:3148783/150‑1 (MQ=255)
gATTTCGTCAATGCGGTAATAGCCACCGTCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                                                                                                                                       <  1:6822885/150‑1 (MQ=255)
gATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                                                                                                                                       <  1:4580590/150‑1 (MQ=255)
gATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                                                                                                                                       >  2:5221314/1‑150 (MQ=255)
gATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                                                                                                                                       <  2:7769423/150‑1 (MQ=255)
  tttCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGa                                                                                                                                                     >  2:4669026/1‑150 (MQ=255)
   ttCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGaa                                                                                                                                                    <  1:2418246/150‑1 (MQ=255)
   ttCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGaa                                                                                                                                                    <  2:7741100/150‑1 (MQ=255)
    tCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAAc                                                                                                                                                   <  2:2550946/150‑1 (MQ=255)
         aaaGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGAcc                                                                                                                                              >  2:2362911/1‑150 (MQ=255)
         aaaGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGAcc                                                                                                                                              >  1:7596644/1‑150 (MQ=255)
          aaGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGAccc                                                                                                                                             <  2:7460937/150‑1 (MQ=255)
          aaGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGAccc                                                                                                                                             <  1:4860148/150‑1 (MQ=255)
           aGCGGTAAAAGCCACCGGCAAGAAGTAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCg                                                                                                                                            <  1:2510247/150‑1 (MQ=255)
           aGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCg                                                                                                                                            >  1:3973041/1‑150 (MQ=255)
            gCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGt                                                                                                                                           <  1:6242717/150‑1 (MQ=255)
            gCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGt                                                                                                                                           <  1:1384249/150‑1 (MQ=255)
            gCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGt                                                                                                                                           <  2:5898322/150‑1 (MQ=255)
             cGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTg                                                                                                                                          >  1:2235565/1‑150 (MQ=255)
                        aCCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTgg                                                                                                                               <  2:6273364/150‑1 (MQ=255)
                        aCCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTgg                                                                                                                               <  1:6315248/150‑1 (MQ=255)
                         ccGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc                                                                                                                              <  1:908474/150‑1 (MQ=255)
                         ccGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc                                                                                                                              <  1:7060200/150‑1 (MQ=255)
                            gCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGAc                                                                                                                           >  1:1069821/1‑150 (MQ=255)
                             cAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACg                                                                                                                          >  2:5399984/1‑150 (MQ=255)
                              aagaagCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGc                                                                                                                         >  1:8762027/1‑150 (MQ=255)
                              aagaagCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGc                                                                                                                         >  2:4999775/1‑150 (MQ=255)
                                   gcagcTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGa                                                                                                                    <  2:3701398/150‑1 (MQ=255)
                                   gcagcTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCACTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGa                                                                                                                    <  2:8238540/150‑1 (MQ=255)
                                      gctgctGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGc                                                                                                                 >  2:7796608/1‑150 (MQ=255)
                                       ctgctgATTGCTGGTGTGGTCACTGCAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCTGCCTGGCGACGCATGGAGGCt                                                                                                                <  2:7157907/150‑1 (MQ=255)
                                       ctgctgATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCt                                                                                                                >  1:6934063/1‑150 (MQ=255)
                                         gctgATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGc                                                                                                              >  2:8561190/1‑150 (MQ=255)
                                         gctgATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCAGTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGc                                                                                                              >  2:6541264/1‑150 (MQ=255)
                                           tgATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCg                                                                                                            <  2:5561972/150‑1 (MQ=255)
                                            gATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCgg                                                                                                           <  1:5031904/150‑1 (MQ=255)
                                             aTTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCggg                                                                                                          <  1:5988497/150‑1 (MQ=255)
                                             aTTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCggg                                                                                                          <  2:177493/150‑1 (MQ=255)
                                              ttGCTGGTGTGGTCACTGAAGTATGCGTGGCCGTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCgggg                                                                                                         <  1:5967203/150‑1 (MQ=255)
                                               tGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGgc                                                                                                        <  2:8953670/150‑1 (MQ=255)
                                                       tgGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtga                                                                                                <  1:828683/150‑1 (MQ=255)
                                                        ggTCCCTGCAGTATGCGTGGCCTGCCGGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgat                                                                                               <  2:7394680/150‑1 (MQ=255)
                                                          tCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgatga                                                                                             <  1:3057691/150‑1 (MQ=255)
                                                            aCTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGAcc                                                                                           >  1:2719785/1‑150 (MQ=255)
                                                              tGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGAccag                                                                                         >  2:4346383/1‑148 (MQ=255)
                                                              tGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTg                                                                                         >  1:6170379/1‑150 (MQ=255)
                                                                  gTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTc                                                                                     <  2:3072066/150‑1 (MQ=255)
                                                                     tGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGc                                                                                  <  1:8526479/150‑1 (MQ=255)
                                                                     tGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGc                                                                                  <  1:7728748/150‑1 (MQ=255)
                                                                     tGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGc                                                                                  <  2:2342030/150‑1 (MQ=255)
                                                                        gTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTg                                                                               <  1:8592907/150‑1 (MQ=255)
                                                                        gTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTg                                                                               >  2:1379419/1‑150 (MQ=255)
                                                                        gTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTg                                                                               >  2:2731776/1‑150 (MQ=255)
                                                                         tGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTgg                                                                              >  1:8045943/1‑150 (MQ=255)
                                                                             cTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTg                                                                          >  2:7540000/1‑150 (MQ=255)
                                                                              ttCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGc                                                                         >  1:580407/1‑150 (MQ=255)
                                                                                  cGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAAc                                                                     <  2:8872146/150‑1 (MQ=255)
                                                                                    gCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTg                                                                   <  1:2077696/150‑1 (MQ=255)
                                                                                    gCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTg                                                                   <  1:5178735/150‑1 (MQ=255)
                                                                                        tGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc                                                                                                                              >  2:7650394/1‑87 (MQ=255)
                                                                                        tGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCAcc                                                               <  2:1089949/150‑1 (MQ=255)
                                                                                        ggTCGGCATGGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCAcc                                                               <  2:7273680/149‑1 (MQ=255)
                                                                                             gCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGAc                                                          <  2:4747334/150‑1 (MQ=255)
                                                                                                    aggaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCa                                                   <  1:7796608/150‑1 (MQ=255)
                                                                                                      gaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAAc                                                 >  2:2818904/1‑150 (MQ=255)
                                                                                                         ggtttCGATGTCTTTGTCGTCACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGat                                              >  2:879724/4‑149 (MQ=2)
                                                                                                            ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                           >  2:5091676/1‑150 (MQ=255)
                                                                                                            ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                           <  1:426111/150‑1 (MQ=255)
                                                                                                            ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                           >  2:548276/1‑150 (MQ=255)
                                                                                                            ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                           >  1:6460969/1‑150 (MQ=255)
                                                                                                               gAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGaa                                        >  1:4496451/1‑150 (MQ=255)
                                                                                                                 ggTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAgg                                      >  1:3923459/1‑150 (MQ=255)
                                                                                                                  gtctttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAggt                                     <  2:6220660/144‑2 (MQ=2)
                                                                                                                     tttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTg                                  >  2:3289974/4‑150 (MQ=2)
                                                                                                                      ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                 <  2:6730346/148‑1 (MQ=2)
                                                                                                                      ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                 <  1:5675025/148‑1 (MQ=2)
                                                                                                                         tcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTa                              >  2:6037330/1‑150 (MQ=255)
                                                                                                                          cgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATACGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTAc                             >  1:741701/1‑150 (MQ=2)
                                                                                                                           gtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAcc                            <  1:2473166/150‑1 (MQ=255)
                                                                                                                           gtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAcc                            <  1:1720886/150‑1 (MQ=255)
                                                                                                                           gtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTAcg                            <  2:8713582/150‑2 (MQ=2)
                                                                                                                              aCCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTg                         >  1:4069318/1‑150 (MQ=255)
                                                                                                                                cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                       >  2:3928323/1‑150 (MQ=255)
                                                                                                                                 gACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTtc                      <  1:3534740/150‑1 (MQ=255)
                                                                                                                                    gCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTcc                   <  2:2100582/150‑1 (MQ=255)
                                                                                                                                     ccTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCa                  <  1:5421054/150‑1 (MQ=255)
                                                                                                                                      ctctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCaa                 <  2:4029641/150‑1 (MQ=255)
                                                                                                                                      ctctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCaa                 >  1:841306/1‑150 (MQ=255)
                                                                                                                                        ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc               <  1:6151342/150‑1 (MQ=255)
                                                                                                                                          ggCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgat                                                                                               >  1:3772210/1‑68 (MQ=255)
                                                                                                                                          ggCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgat                                                                                               <  2:3772210/68‑1 (MQ=255)
                                                                                                                                             atgTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCt                          >  2:8434193/3‑134 (MQ=255)
                                                                                                                                             aCGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATc          >  2:6911451/1‑150 (MQ=255)
                                                                                                                                               gTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATccc        <  1:2573212/150‑1 (MQ=255)
                                                                                                                                                 tCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCgg      >  2:4895325/1‑150 (MQ=2)
                                                                                                                                                 tCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCgg      >  1:4679173/1‑150 (MQ=2)
                                                                                                                                                  cAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGa     <  1:8915581/150‑1 (MQ=2)
                                                                                                                                                  cAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGa     <  1:5529975/150‑1 (MQ=255)
                                                                                                                                                  cAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGAGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGa     <  2:7013115/150‑1 (MQ=255)
                                                                                                                                                    aCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACt   >  2:4251018/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTA  >  NC_002947/827095‑827393

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