Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,079 T→G 9.3% intergenic (+40/+112) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,0790TG 9.3% 424.9 / 3.2 170intergenic (+40/+112)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (57/90);  new base G (19/2);  total (76/94)
Fisher's exact test for biased strand distribution p-value = 6.95e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

TGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCA  >  NC_002947/667930‑668228
                                                                                                                                                     |                                                                                                                                                     
tGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGT                                                                                                                                                       >  2:3934992‑M1/1‑150 (MQ=255)
 gACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTC                                                                                                                                                      <  2:1345284‑M1/150‑1 (MQ=255)
    cccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGGGTCGCG                                                                                                                                                   >  1:3283733‑M1/1‑150 (MQ=255)
         gCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAgg                                                                                                                                              <  1:598553/150‑1 (MQ=255)
             cTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGGGTGCGGCGCGGTGGGGCCg                                                                                                                                          >  1:1015572/1‑150 (MQ=255)
             cTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGCTGGGGCCg                                                                                                                                          >  1:3177926/1‑150 (MQ=255)
                ccGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaa                                                                                                                                       >  1:2013388/1‑150 (MQ=255)
                  gcgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGTGTCGCGATGGGGCCGctttg                                                                                                                                     <  2:3624922/150‑5 (MQ=255)
                   cgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                                    <  1:3960622/150‑1 (MQ=255)
                   cgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                                    <  2:2974804/150‑1 (MQ=255)
                   cgTGCCGCCGCCGCAGGTCGAATCGATGGACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                                    <  2:3444404/150‑1 (MQ=255)
                    gTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCg                                                                                                                                   <  2:257053/150‑1 (MQ=255)
                     tGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGTCAAGGGGAGATGACTGTGAAAGCCCGCTAGATACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCgg                                                                                                                                  <  2:2600261/150‑1 (MQ=255)
                      gccgccgccGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagc                                                                                                                                 <  1:1966961‑M1/150‑18 (MQ=255)
                      gccgccgccGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                                 <  1:1527171/150‑1 (MQ=255)
                          ccgccgcAGGTCGAATCGATGGACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGcccta                                                                                                                             <  1:2864001/150‑3 (MQ=255)
                           cgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGa                                                                                                                            <  2:895410/150‑1 (MQ=255)
                           cgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGa                                                                                                                            <  2:3001383/150‑1 (MQ=255)
                             ccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctatcg                                                                                                                          <  2:1856619/150‑6 (MQ=255)
                              cgcCGGTCGAATCTTTGTAGCGGCCTAACCTGCAGGCCGTGGGCTTTGGGCATGGGGAGATGAGGGTGAACGCCCGCTAGAGCCCCGCTCACGGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATTGGGCCGAAATGCGGCCCCGACaa                                                                                                                         <  2:433604/150‑1 (MQ=255)
                                   gTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGGGCGGGGCGGTAGGGCCGCCAAGCGGCCCCCACaaagaaa                                                                                                                    >  1:2805082/1‑145 (MQ=255)
                                    tCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCCCGAAAGGGCCGCAAACCGGCCCCCACAAATCaac                                                                                                                   >  2:1240438/1‑149 (MQ=255)
                                      gAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATATGGCCGCCAAGCGGCCCCGCCCATTCAAgcg                                                                                                                 >  1:4114361/1‑150 (MQ=255)
                                       aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcgg                                                                                                                <  2:5384063/150‑1 (MQ=255)
                                           gATGTACCGGCCGAACCTGCAGGCGGTGGGCGTTGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                            <  2:2013388/150‑1 (MQ=255)
                                                 ccGGCCGCACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaa                                                                                                      <  2:4952926/150‑1 (MQ=255)
                                                   ggCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAg                                                                                                    >  2:3126717/1‑150 (MQ=255)
                                                       gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                <  2:4009135/150‑1 (MQ=255)
                                                       gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcggc                                                                                                               >  1:5354457/1‑135 (MQ=255)
                                                       gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                <  2:605854/150‑1 (MQ=255)
                                                       gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCGGAAAGCCTg                                                                                                >  2:1309580/1‑150 (MQ=255)
                                                        ccccTGCCGGCCGGGGGCGCTGGGCACGGGGAGAGGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGGGGGCGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                               <  2:1502895/148‑1 (MQ=255)
                                                        acccTGCAGGCCGTGGGCGTTGGGCATGGGGAGATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGTAATGTGTGGGAGCGGCCTTGCGTCGCTATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                               <  1:4858314/148‑1 (MQ=38)
                                                        aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                               <  2:1503142/150‑1 (MQ=255)
                                                         aCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTggg                                                                                              <  2:2736370/150‑1 (MQ=255)
                                                           cTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGCTAGGGCCCCCAAGCCGCCCCGACAAGTCAAGCGGCGAAGCGGGAAGCCGGGGgc                                                                                            >  2:3904337/1‑150 (MQ=255)
                                                            ggCAGGCCGTGGGCGATGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGACAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAGTCAAGCGGCGAAGCTGAAAGCCTGGGgct                                                                                           <  2:5452727/149‑1 (MQ=255)
                                                               aGGGCGTGGGCGATGGGCAAGGGTTGTTGTGGGTGAAAGCCCGCTTGATCCCCGGTATACGGTGAATGTGTGGGTGCGGGCTTGGGTCGCGTTGGGGCCGCAAAGGGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                        <  1:397866/150‑1 (MQ=255)
                                                               aGGCCGTGGGCGATGGGCATGGGTAGATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCTATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                        <  1:2606763/150‑1 (MQ=255)
                                                                 gCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgt                                                                                      <  1:4606702/150‑1 (MQ=255)
                                                                 gCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGCCCATTCAAGCGGCGAAGCTGTAAGCCGGGGGCTGCtgt                                                                                      >  2:4469617/1‑150 (MQ=38)
                                                                  ccGTGGGGTTTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                     <  1:3479100/150‑1 (MQ=255)
                                                                  ccGTGGGCTATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                     <  1:5569447/150‑1 (MQ=255)
                                                                  ccGTGGGCGTTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGGGGAATGTGTGGGAGCGGCCTTGCGCCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                     <  1:3870283/150‑1 (MQ=255)
                                                                    gggggCGATGGGCACGGGGAGATGACGGTGAACGCCCGCTAGACACCCGCTCTACGGGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:5312035/148‑1 (MQ=38)
                                                                    gTGGGCTTTGGGAATGGGGAGATGACGGTTAAAGCCCGCTAGACACCCGCTCTACGCGTAAAGTTTGGGAGCGGCATTGCTTCGCTATAGGCCCCAAAAGCGCCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:5207688/150‑1 (MQ=255)
                                                                    gTGGGCTATGGGCAAGGGTAGATGACGTTGAAAGCCCGCTAGATACCCGCTCTACGCGTAAAGTTTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:2233514/150‑1 (MQ=38)
                                                                    gTGGGCGTTGGGCAAGGGGATATGACGGTGAACGCCCGCTAGAGACCCGCTCTCCGGGGAAAGTGTGGGCGCGGCCTTGGGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:1046523/150‑1 (MQ=255)
                                                                    gTGGGCGTGGGGCATGGGGAGATGCGGGTTAAAGCCCGCTAGAGACCCGCTCTAGGCGGAAAGTGTGGGCGCGGCCTGGCGTCGCTATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:477941/150‑1 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  1:3165011/150‑1 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAtagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcg                                                                                   <  1:1556057‑M1/150‑64 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   <  2:4358450/150‑1 (MQ=255)
                                                                    gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                   >  2:3304674/1‑150 (MQ=255)
                                                                     tGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATTGGGCCGCTAAGCGGCCCCGACAATTCAAGCGGCGATGCTGAAAGCCTGGGGCTGCTGTGCAg                                                                                  <  2:457741/150‑1 (MQ=255)
                                                                         cgttGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                              <  2:4166679/147‑1 (MQ=255)
                                                                         cGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTAGGGCTGCACAGCAGcccc                                                                              <  1:2089691/150‑2 (MQ=38)
                                                                         cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGGTTGGGGGGGGGTTGGGGCGGGTAGGGGGGCCGGCAATTCAAGCGGGGAAGCTGATAGGCTGGGGGTGGTGTGGGGCCCt                                                                              >  1:4671530/1‑150 (MQ=255)
                                                                         cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGCTAGGGCCGCCAAGCGGCCCCGACAAATCAAGCGGCGAAGCGGCAAGCCTGGGGCTGCTGTGCAGcccg                                                                              >  2:4847493/1‑149 (MQ=255)
                                                                          gATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATCAAGCGGCGCAGCTGAAAACCAGGGGCAGCGGTGCAGcccag                                                                             >  1:4617969/1‑148 (MQ=38)
                                                                           ctgGGCACGGGGAGCTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgttttgcggccctat                                                                            <  1:2834835‑M1/149‑71 (MQ=255)
                                                                            tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatc                                                                           <  1:300473‑M1/150‑72 (MQ=255)
                                                                             gggCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcg                                                                          >  1:3677015‑M2/1‑147 (MQ=255)
                                                                               gCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcg                                                                        <  1:1297734‑M1/150‑75 (MQ=255)
                                                                               gCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGGGGGTGGCGATAGGGCCGCAAAGCGGGCCCGACACATCAAGCGGCGCAGCGGAAAGCCTGGGGCTGCTGTGCGGCCCAGtcgcg                                                                        <  1:5040422‑M2/150‑6 (MQ=255)
                                                                               gCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCCTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGCCAATTCAAGCGGCGAAACGGAAAGCCGGGGGGAGCTGTGCCGCCCGGtcgcg                                                                        <  2:475246‑M2/150‑6 (MQ=255)
                                                                                catggggAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGa                                                                       <  1:1609651/147‑1 (MQ=255)
                                                                                 aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgac                                                                      <  1:4408169‑M1/150‑77 (MQ=255)
                                                                                   ggggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgc                                                                    <  2:5531117‑M2/150‑10 (MQ=255)
                                                                                     ggtgaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCaa                                                                  <  1:4522785/147‑1 (MQ=255)
                                                                                       tgatGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCTATCGGGCCGCACAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                >  1:734476‑M2/2‑137 (MQ=255)
                                                                                       agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                >  1:2171185‑M2/1‑137 (MQ=255)
                                                                                       agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                >  1:2171188‑M2/1‑137 (MQ=255)
                                                                                        gaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCAACGCGGCCCCGACAATTCCCGCGGCGACGCGGAAAGCCTGGGGCGGCGGTGCCGCCCTGtcgcgacgcaaggc                                                               <  2:5106950‑M2/150‑15 (MQ=255)
                                                                                              ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGCTAGGGCCGCAAAGCGGGCCCGCCACGTCAAGCGGCGGAGCGGGGAGCCGGGGGCGGCGGTGGAGCCCGGTCGCGGCGCGCGGCCGCGcc                                                         >  2:1640877/1‑150 (MQ=255)
                                                                                                       ccGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCACGGCCGCTCCCACAAAGAc                                                >  2:3621885/1‑150 (MQ=255)
                                                                                                        cGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCGGAAAGCCTGGGGCTGCGGTGCCGCCCTGTCGCGACGCCAGGCCGCTCCCACAAAGAcc                                               >  2:1834852/1‑150 (MQ=255)
                                                                                                                     atctACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGGGCAGCCCTGTCGCGACGCCAGGCCGCTCCCCCAAAGACCGCGCAGATCAGTg                                  >  2:2168910/2‑150 (MQ=255)
                                                                                                                        tACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGAGGCGAAGCAGAAACCCAGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCGCTCCCACAAATACCGCGCAGATCAGTGGAt                               >  2:4206578/1‑150 (MQ=25)
                                                                                                                        tACGCGGAAAGTGTGGGAGCGGCCGTTCGTGGCGGTATGGCCCCCAAGCGGCCCCCACCAATCAACCGGCGAAGCTGTAAGCCTGGGGCTGGGGGGCAGCCCTGGCGCGACGCAAGGCCGCTCCCCCCAAGCACGCGCCGCTAAGTGGAt                               >  2:5172560/1‑150 (MQ=255)
                                                                                                                        tACGCGGAAAGTGTGGGAGCGGCCATGCCTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATAAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGGGCAGCCCAGTCGCGACGCAAGGCCGCTCCCCCCAAGACCGCGCCCATAACTGGAt                               >  2:5172604/1‑150 (MQ=25)
                                                                                                                             ggAAAGTGTGGGCGCGGCGTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCCAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGATGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCAccc                          <  2:5097144/150‑1 (MQ=37)
                                                                                                                                     tgggtgggGCCTTGGGTCGCTActgggctgcacagcagcccctggctttcagcttcgccgcttgatttgtcggggccgcttggcggccatatggcgacgaaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttca                  <  2:1979007‑M1/150‑129 (MQ=255)
                                                                                                                                     tgGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTAAAGCGGCGAAGCTGAAAGCCTTGGGCTGCAGTGCAGCCCAGTCGCGACGCAAAGCCGCGCCCACAAAGACCCCGCACATAAGTGGAACACCCAGcacaac                  >  2:2681591/1‑148 (MQ=25)
                                                                                                                                       ggAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTAAAGCGGCGAACCTGAAGGCCAGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACCAAGACCGCGCAGATCAGTGGATAACCCACCTCAacac                >  2:1228058/1‑150 (MQ=37)
                                                                                                                                             gCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatc          >  2:2781507‑M1/1‑14 (MQ=255)
                                                                                                                                             gCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgcgttgcggcccaatcgcgacgcaaggccgctcccacacagtccgcggagagcgggtctctatcgggcgttcacccgcaac          >  1:2819441‑M1/1‑14 (MQ=255)
                                                                                                                                             gCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg          >  2:4414966/1‑150 (MQ=255)
                                                                                                                                             gCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg          >  2:3841226/1‑150 (MQ=255)
                                                                                                                                             gCCTTGCGTCGCGATACGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGTCCTCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAACCACAg          >  2:3841560/1‑150 (MQ=255)
                                                                                                                                              ccTTGCGTCGCTATTGGCCCGCAACGCGGCCCCGAAAATTCAAGCGGCGAAGCTTAAAGCGTGGGGCTGGTGTGCAGCCCTGTCGCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGc         <  2:4355945/150‑1 (MQ=255)
                                                                                                                                              ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGc         >  2:2310893/1‑150 (MQ=255)
                                                                                                                                              ccTTGAGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacacattccgcgtagagagggtctctagcgggctttcaccgtcatct         >  2:4561160‑M1/1‑13 (MQ=255)
                                                                                                                                               cTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgcttttcggccctatcgcgacgcaaggccgctcccacacgttacccgtagatcgggtctctagcgggggttcaccgtcaaata        >  1:228253‑M1/1‑12 (MQ=255)
                                                                                                                                                ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGTGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc       >  1:8139/1‑150 (MQ=255)
                                                                                                                                                 agcGTCGCGAcagggctccacaccacccccaggatatcagcttcgccgctcgaattgtcggggccgctttgcggcccaaccgcgacgaaagcccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccc      <  2:3617344‑M1/149‑141 (MQ=255)
                                                                                                                                                  gCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcccc     <  2:5531117‑M1/150‑142 (MQ=255)
                                                                                                                                                    gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc   >  1:5440040/1‑150 (MQ=255)
                                                                                                                                                     tCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttg  >  2:2193562‑M1/1‑6 (MQ=255)
                                                                                                                                                     tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCa  >  1:5606731/1‑150 (MQ=255)
                                                                                                                                                     tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCa  >  1:2797237/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCA  >  NC_002947/667930‑668228

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