Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 122,222 A→C 6.4% intergenic (+32/+29) katE → / ← PP_0116 hydroperoxidase/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947122,2220AC 6.4% 224.1 / 2.0 72intergenic (+32/+29)katE/PP_0116hydroperoxidase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base C (2/3);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCA  >  NC_002947/122074‑122369
                                                                                                                                                      |                                                                                                                                                    
gaggaCAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGa                                                                                                                                                      >  5:224044/1‑149 (MQ=255)
   gaCAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATg                                                                                                                                                    >  6:327705/1‑148 (MQ=255)
      aaGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAg                                                                                                                                                 >  7:13041/1‑148 (MQ=255)
           gCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCgg                                                                                                                                           >  6:235379/1‑149 (MQ=255)
           gCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCgg                                                                                                                                           <  3:291480/149‑1 (MQ=255)
            cTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCggg                                                                                                                                          <  8:13041/149‑1 (MQ=255)
                        gATGAGCAGAAGGCAGTGGATGCCTTTGTCAGAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCcgcg                                                                                                                                >  2:103359/1‑147 (MQ=255)
                         aTGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAg                                                                                                                             >  5:244220/1‑149 (MQ=255)
                          tGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGc                                                                                                                            <  7:278964/149‑1 (MQ=255)
                            aGCAGAAGGGAGGGGAGGCCTTTGTCAAAGCGGTAGAAGGCCTCCTGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCGGTTCATGCCTGACAAG‑‑GCCTTTTCGCGGGTAAAGCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCg                                                                                                                           <  5:100356/148‑1 (MQ=255)
                             gCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGtt                                                                                                                         >  7:132569/1‑149 (MQ=255)
                              cAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGaa                                                                                                                              >  5:141487/1‑143 (MQ=255)
                              cAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGaa                                                                                                                              <  6:141487/143‑1 (MQ=255)
                                gAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACGCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTAc                                                                                                                      <  2:250586/149‑1 (MQ=255)
                                gAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTAc                                                                                                                      >  8:217307/1‑149 (MQ=255)
                                  aGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTActg                                                                                                                    <  8:132569/149‑1 (MQ=255)
                                   ggCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTActgc                                                                                                                  <  2:91346/149‑1 (MQ=255)
                                    gCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCGGATGTAGGAGCCGGTTTACCCGCGCAGCGTTTActgct                                                                                                                 >  1:378012/1‑149 (MQ=255)
                                     cAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTActgct                                                                                                                 >  6:334599/1‑148 (MQ=255)
                                      aGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGCCTCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTActgct                                                                                                                 <  1:113340/149‑1 (MQ=255)
                                        tGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTTAAACCCCGACGGATGTAGGGGCGGGTTTAACCCCGAAGCGTTTActgggttgt                                                                                                              >  8:54137/1‑143 (MQ=255)
                                              ccTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                                                                       <  7:217307/149‑1 (MQ=255)
                                               cTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCCGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAcc                                                                                                      >  2:350465/1‑149 (MQ=255)
                                                tttGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                                                                       >  6:118965/1‑147 (MQ=255)
                                                tttGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                                                                       >  5:359817/1‑147 (MQ=255)
                                                    tCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGca                                                                                                  <  5:235379/148‑1 (MQ=255)
                                                     cAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCTGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcac                                                                                                 <  6:244220/148‑1 (MQ=255)
                                                      aaaGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTTGGGGCGGGGTTTCCCGCCAAGCGTTTAATGCTGAATCGGGACCCGcacca                                                                                               >  1:378164/1‑149 (MQ=255)
                                                        aGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCaccacc                                                                                             >  7:13596/1‑149 (MQ=255)
                                                         gCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCt                                                                                            <  6:323542/149‑1 (MQ=255)
                                                         gCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCt                                                                                            <  5:334599/149‑1 (MQ=255)
                                                           ggTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCAGATCCGGGACCAGcacaaacggt                                                                                          >  6:265944/1‑142 (MQ=255)
                                                                  ggCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCgg                                                                                     >  7:108746/1‑147 (MQ=255)
                                                                      aCCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCagcga                                                                               <  8:108746/149‑1 (MQ=255)
                                                                        cGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCa                                                                                   <  5:36083/143‑1 (MQ=255)
                                                                        cGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCa                                                                                   >  6:36083/1‑143 (MQ=255)
                                                                            tgtgGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTg                                                                         <  6:224044/149‑1 (MQ=255)
                                                                               gggAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCtgt                                                                                          >  1:340716/1‑129 (MQ=255)
                                                                               gggAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCtgt                                                                                          <  2:340716/129‑1 (MQ=255)
                                                                               gggAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGAt                                                                      >  6:369802/1‑149 (MQ=255)
                                                                                ggAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATc                                                                     >  7:11744/1‑149 (MQ=255)
                                                                                    cGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAt                                                                 <  5:369802/149‑1 (MQ=255)
                                                                                         ggcggcggcTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGGGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgcttc                                                            <  5:327705/149‑1 (MQ=255)
                                                                                           cggcggcTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgcttctg                                                          >  7:243767/1‑149 (MQ=255)
                                                                                              cggcTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgcttctgctt                                                       <  1:350465/149‑1 (MQ=255)
                                                                                                  tGAAGCTGTTCCTGCCTGACAAT‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGc                                                   <  6:358886/149‑1 (MQ=255)
                                                                                                   gAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCt                                                                                            <  4:376500/107‑1 (MQ=255)
                                                                                                   gAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCt                                                                                            >  3:376500/1‑107 (MQ=255)
                                                                                                    aaGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGc                                                   <  6:359817/147‑1 (MQ=255)
                                                                                                    aaGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTc                                                 <  7:54137/149‑1 (MQ=255)
                                                                                                             ctgcctgACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACt                                        <  3:375641/149‑1 (MQ=255)
                                                                                                                    aCAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCTCCACCTGTGCCGGCAGCGAGCTCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAAt                                  <  5:265944/148‑1 (MQ=255)
                                                                                                                     cAAG‑‑GCCTTTTCGGGGGTAAAGCCCCTCCTGCTGTCGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCt                                <  1:376535/149‑1 (MQ=255)
                                                                                                                     cAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCt                                                                          <  4:101316/107‑1 (MQ=255)
                                                                                                                     cAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCt                                                                          >  3:101316/1‑107 (MQ=255)
                                                                                                                      aaG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCtt                               >  6:384179/1‑149 (MQ=255)
                                                                                                                        acgcccTTTCCGCGGGAAACCCCCGCCCTCTGGATGAGGCGGTTTAACCGGGTATGGTTTACTGCCTCTTGCGGCACCCCAACAACCGGGCCTGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCtt                               <  1:103359/145‑1 (MQ=255)
                                                                                                                         taaacgcttCGCGGGTAAACCCGCTCCTACATCAGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTg                              >  2:23268/8‑149 (MQ=25)
                                                                                                                           gCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTgcg                            >  4:96945/1‑149 (MQ=255)
                                                                                                                            acgcttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTgcgc                           <  5:1396/145‑1 (MQ=255)
                                                                                                                              ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcaccac                                                                                              <  8:351061/80‑1 (MQ=255)
                                                                                                                              ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcaccac                                                                                              >  7:351061/1‑80 (MQ=255)
                                                                                                                              ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGa                                                                       >  2:131653/1‑103 (MQ=255)
                                                                                                                              ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGa                                                                       <  1:131653/103‑1 (MQ=255)
                                                                                                                              ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCAc                         >  8:117961/1‑149 (MQ=255)
                                                                                                                                   gcgGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACGTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCGTCAGCTCGATCCCACTGTCAATCTTGCGCACgcgcg                    >  8:88144/1‑148 (MQ=255)
                                                                                                                                   gcgGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACgcgc                     <  7:117961/148‑1 (MQ=255)
                                                                                                                                      ggTAAACCCGCTCCTACATCAGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCtttt                 >  1:106845/1‑149 (MQ=38)
                                                                                                                                       gTAAACCCGCTCCTACATCAGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTc                <  1:23268/149‑1 (MQ=38)
                                                                                                                                          accccccTCCTGATGTAGGAGCGGGTTAACCCTCGAAGCGTTTACTGCTGATGTGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCa              <  8:122206/146‑1 (MQ=255)
                                                                                                                                          aaCCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCg                                                <  4:136054/114‑1 (MQ=255)
                                                                                                                                          aaCCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCg                                                >  3:136054/1‑114 (MQ=255)
                                                                                                                                            cccccTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGc            >  4:6037/1‑148 (MQ=255)
                                                                                                                                                      aTGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCa  <  5:118965/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAG‑‑GCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCA  >  NC_002947/122074‑122369

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