Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 949,138 C→A 5.9% intergenic (+433/‑24) hmp → / → dsbB flavohemoprotein/disulfide bond formation protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947949,1380CA 5.9% 216.4 / 2.3 67intergenic (+433/‑24)hmp/dsbBflavohemoprotein/disulfide bond formation protein B
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base A (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TTTGCCGCCCGGCGTTGAACCGACGTCGTTGCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAGGCACCCTGCCCG  >  NC_002947/948990‑949287
                                                                                                                                                     |                                                                                                                                                     
tttGCCGCCCGGCGTTGAACCGACGTCGTTGCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCa                                                                                                                                                      <  1:125970/150‑1 (MQ=255)
      gccCGGCGTTGAACCGACGGCGTTGCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCAGCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCTGACCACGCGGCCACTATcc                                                                                                                                                <  2:1734808/150‑1 (MQ=255)
        ccGGCGTTGAACCGACGTCGTTGCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCAt                                                                                                                                              <  2:41145/150‑1 (MQ=255)
         cGGCGTTGAACCGACGTCGTTGCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATc                                                                                                                                             >  2:1264918/1‑150 (MQ=255)
             gTTGAACCGACGTCGTTTCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCAACTATCCATCGAAg                                                                                                                                         >  2:223052/1‑150 (MQ=255)
             gTTGAACCGACGTCGTTGCGGCCGGTCTCAGACACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCAACTATCCATCGAAg                                                                                                                                         >  2:223011/1‑150 (MQ=255)
               tGAACCGACGCCGTTGCGGCCGGTCTCAGCCCCACTGGCTAGGCCTCCGCGAGCATCCCCCGACGCCCCCTCCGCCTGCTCACAGGCCATGCCCAGCGCGTCGACTTACGCCTAGGCAGCCGGCCCACGCAGCCACTACCCATAAAAggg                                                                                                                                       >  2:270604/1‑149 (MQ=255)
                    cGACGTCGTTGCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCgg                                                                                                                                  >  1:1300901/1‑150 (MQ=255)
                              gCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATgaacga                                                                                                                        <  2:1240044/150‑1 (MQ=255)
                                ggCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCCTGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATgaacgaac                                                                                                                      <  1:434794/150‑1 (MQ=255)
                                   cGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGaacaaa                                                                                                                   <  1:1186967/150‑1 (MQ=255)
                                       ctcAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATc                                                                                                               <  2:492157/150‑1 (MQ=255)
                                        tcAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATcg                                                                                                              >  1:879685/1‑150 (MQ=255)
                                           gCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCgcgc                                                                                                           <  2:311532/150‑1 (MQ=255)
                                           gCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCgcgc                                                                                                           >  1:672212/1‑150 (MQ=255)
                                            ccACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGTGTGTAGACTTGCGCCCCGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGcc                                                                                                          >  2:948973/1‑150 (MQ=255)
                                             cacaCTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCt                                                                                                         <  2:650142/150‑1 (MQ=255)
                                                 cTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAAt                                                                                                     <  2:599466/150‑1 (MQ=255)
                                                   ggATAGCCCTACCTCATGCTTCCGCGCTCATCCCTCGCAATGCTAGCCGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAACCATCGCGCCTGAATCg                                                                                                   <  2:803467/150‑1 (MQ=255)
                                                       aGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGaa                                                                                               <  2:1489265/150‑1 (MQ=255)
                                                        gCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGAAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAAc                                                                                              <  2:768567/150‑1 (MQ=255)
                                                         cccTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACg                                                                                             >  1:504587/1‑150 (MQ=255)
                                                         cccTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACg                                                                                             <  2:1159713/150‑1 (MQ=255)
                                                                  gTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCt                                                                                    <  1:1411674/150‑1 (MQ=255)
                                                                  gTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCt                                                                                    <  1:847964/150‑1 (MQ=255)
                                                                   tGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCTAAGTAACCGGCAATGAACGAACAACCATCGCGCCTGAATCGGGAACGGCGCTTTCTg                                                                                   <  1:948973/150‑1 (MQ=255)
                                                                   tGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTg                                                                                   <  2:541343/150‑1 (MQ=255)
                                                                   tGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTg                                                                                   <  2:1293931/150‑1 (MQ=255)
                                                                       ttCCGCGACGCTCCCTCTGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCAACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTCAATCGGGAACGGCGCTTTCTGGtg                                                                                <  1:223011/150‑1 (MQ=255)
                                                                       ttCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCAACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGtg                                                                                <  1:223052/150‑1 (MQ=255)
                                                                              tcGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAAGCGGGAACGGCGCTTTCTGGTGCTGCTgg                                                                         <  1:270604/149‑1 (MQ=255)
                                                                              aCGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTgg                                                                         <  2:531995/150‑1 (MQ=255)
                                                                              aCGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTgg                                                                         <  2:504587/150‑1 (MQ=255)
                                                                                    ccTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATc                                                                                                    >  1:1713702/1‑117 (MQ=255)
                                                                                    ccTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATc                                                                                                    <  2:1713702/117‑1 (MQ=255)
                                                                                             gCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCttt                                                          <  1:390942/150‑1 (MQ=255)
                                                                                               tCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCg                                                                                                   >  2:634418/1‑107 (MQ=255)
                                                                                               tCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCg                                                                                                   <  1:634418/107‑1 (MQ=255)
                                                                                               tCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCa                                                        <  2:1300901/150‑1 (MQ=255)
                                                                                                 acaGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACt                                                      >  1:375351/1‑150 (MQ=255)
                                                                                                     gCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCggg                                                                                                 >  2:1142208/1‑103 (MQ=255)
                                                                                                     gCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCggg                                                                                                 <  1:1142208/103‑1 (MQ=255)
                                                                                                        aTGCCCAGCGTGTAGACTTGCGCCCAAGCAGCCCGCCCACCATGCCACTACCGATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATcggc                                               >  1:243654/1‑150 (MQ=255)
                                                                                                          gCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATcggcgg                                             >  2:413583/1‑150 (MQ=255)
                                                                                                             cAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGgcgc                                          >  1:1314469/1‑150 (MQ=255)
                                                                                                             cAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCACGCCTGAATCGGGAACGGCACTTTCTGGTGCTGCTGGGGCTGATCTGCCTTGCACTGATCGGCGgcgc                                          >  2:409297/1‑150 (MQ=255)
                                                                                                              aGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGcc                                         >  1:739804/1‑150 (MQ=255)
                                                                                                                  tgtAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCctct                                     >  2:1077865/1‑150 (MQ=255)
                                                                                                                   gtAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCg                                                                                                   <  1:930676/87‑1 (MQ=255)
                                                                                                                   gtAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCg                                                                                                   >  2:930676/1‑87 (MQ=255)
                                                                                                                   gtAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTa                                    <  1:1077865/150‑1 (MQ=255)
                                                                                                                         ttGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCa                              >  1:586797/1‑150 (MQ=255)
                                                                                                                         tggcgcCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCa                              >  2:1485737/3‑150 (MQ=255)
                                                                                                                             gcCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAggtg                          <  2:243654/150‑1 (MQ=255)
                                                                                                                              cccAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAggtgg                         >  2:1428454/1‑150 (MQ=255)
                                                                                                                               ccAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAggtggt                        <  1:1485737/150‑1 (MQ=255)
                                                                                                                               ccAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAggtggt                        >  2:833806/1‑150 (MQ=255)
                                                                                                                                   gcaggcGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTc                    >  1:715334/1‑150 (MQ=255)
                                                                                                                                        cGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGa               <  1:1264918/150‑1 (MQ=255)
                                                                                                                                         ggACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAg              <  2:1344548/150‑1 (MQ=255)
                                                                                                                                          gACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAgg             <  2:879685/150‑1 (MQ=255)
                                                                                                                                              aCGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCAc                                                       <  2:479490/104‑1 (MQ=255)
                                                                                                                                              aCGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCAc                                                       >  1:479490/1‑104 (MQ=255)
                                                                                                                                              aCGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAGGCAcc         >  2:434246/1‑150 (MQ=255)
                                                                                                                                                  tGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTgg                                                                         >  1:837900/1‑82 (MQ=255)
                                                                                                                                                  tGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTgg                                                                         <  2:837900/82‑1 (MQ=255)
                                                                                                                                                  tGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAGGCACCCTGc     >  2:1095890/1‑150 (MQ=255)
                                                                                                                                                    ccACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAGGCCCGCTGccc   <  2:672212/150‑1 (MQ=255)
                                                                                                                                                     cACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAGGCACCCTGCCCg  >  2:68966/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGCCGCCCGGCGTTGAACCGACGTCGTTGCGGCCGGTCTCAGCCACACTGGATAGCCCTGCGCAGTGCTTCCGCGACGCTCCCTCGGCATGCTCACAGGCCATGCCCAGCGTGTAGACTTGCGCCCAGGCAGGCGGACCACGCTGCCACTATCCATCGAAGGAACCGGCAATGAACGAACAAACATCGCGCCTGAATCGGGAACGGCGCTTTCTGGTGCTGCTGGGGCTGATCTGCCTTTCACTGATCGGCGGCGCCCTCTACATGCAGGTGGTGCTCGGCGAGGCACCCTGCCCG  >  NC_002947/948990‑949287

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