Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,303,359:1 +G 6.3% coding (414/3675 nt) recB ← Chi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,303,3591.G 6.3% 245.2 / 11.4 63coding (414/3675 nt)recBChi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base . (32/27);  new base G (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TGCCGCGCACCCAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGTGGATCGGCGCC  >  NC_002947/5303214‑5303502
                                                                                                                                                       |                                                                                                                                               
tGCCGCGCACCCAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCgcag                                                                                                                                               >  1:370543/1‑148 (MQ=255)
  ccGCGCACCCAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCAt                                                                                                                                            >  2:307506/1‑149 (MQ=255)
       cACCCAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTGCGCGAGCATgcctt                                                                                                                                       <  2:232712/149‑4 (MQ=255)
          ccAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTgg                                                                                                                                     >  1:252239/1‑148 (MQ=255)
          ccAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTgg                                                                                                                                     >  5:211646/1‑148 (MQ=255)
          ccAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGc                                                                                                                                    <  8:218330/149‑1 (MQ=255)
              gCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCa                                                                                                                                <  7:334063/149‑1 (MQ=255)
               ccAGTGAGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGCTTTCCAGTGTCTGGGT‑‑GAACAGGCTCCCGCAGTCGAAGGCATGCTCGCGCAGCATGCCCTCGCACCAg                                                                                                                               >  4:269954/1‑149 (MQ=255)
                          ccTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGAGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCAAGGGTCTGGGT‑‑GAACAGGCTGCAGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGt                                                                                                                    >  5:192021/1‑149 (MQ=255)
                               cATCCCATAGCAGAAGGGCCGCCAGTAATCGGGCTTGACCTGGCCCAACTACTCGCTGTGGGCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATa                                                                                                               <  2:168510/149‑1 (MQ=255)
                                aTCCCATTGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCATGGTCTGGGT‑‑GAACAGGCTTCCGCTGTCGTAGGCATGCTCGCGCAGCATGCGCTGGCACCATCCATGGATGGTGGATAc                                                                                                              <  3:269954/149‑1 (MQ=255)
                                aTCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATAc                                                                                                              >  8:325273/1‑149 (MQ=255)
                                aTCCCATAGCAGAAGCGCCGCCAATAAAAGCGCATTACCTGGCCCAACAAATAGATATGGTCGGTTTCAAAGGTATGAGT‑‑GAACAGGCTGCCGCTGACGAAGGCATGCACGCGCAGAATgc                                                                                                                                          >  7:56911/1‑121 (MQ=255)
                                 tCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATAcc                                                                                                             <  5:285418/149‑1 (MQ=255)
                                   ccATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGc                                                                                                           <  4:207495/149‑1 (MQ=255)
                                    cATAGCAGAAGCGCCGCCAGTAATAGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCg                                                                                                            >  5:169466/1‑147 (MQ=255)
                                       aGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGcc                                                                                                       >  5:126976/1‑149 (MQ=255)
                                        gCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACGTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGcc                                                                                                       >  1:72860/1‑149 (MQ=255)
                                         gagtatcgccgcGAGTTATAGCGCATGTCCTTGCCCAAAAACTCTGTGTGGTGGGTTTCTCATGGTTTTGGCTCCTTACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCa                                                                                                          <  1:307506/142‑1 (MQ=255)
                                         cAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTc                                                                                                     >  6:346358/1‑149 (MQ=255)
                                         cAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTc                                                                                                     <  1:277831/149‑1 (MQ=255)
                                           gatccgcctcCAGTAATCGCGCACGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAATAGGCTTCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGt                                                                                                   <  8:347880/145‑1 (MQ=255)
                                           gAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGt                                                                                                   <  2:252239/149‑1 (MQ=255)
                                                    ccAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACt                                                                                          <  7:48086/149‑1 (MQ=255)
                                                    ccAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACt                                                                                          <  6:126976/149‑1 (MQ=255)
                                                     cAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTc                                                                                                                   >  2:156445/1‑123 (MQ=255)
                                                     cAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTc                                                                                                                   <  1:156445/123‑1 (MQ=255)
                                                     cAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTg                                                                                         >  3:238609/1‑149 (MQ=255)
                                                      aGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGc                                                                                        >  5:304443/1‑149 (MQ=255)
                                                                   gACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTcc                                                                            >  8:46036/1‑148 (MQ=255)
                                                                           ccAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAAccggccg                                                                   >  6:43854/1‑149 (MQ=255)
                                                                            caacaaCTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAAccggccgg                                                                  >  2:337551/1‑149 (MQ=255)
                                                                              acaacTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGAt                                                                                >  8:316779/1‑133 (MQ=255)
                                                                              acaacTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGAt                                                                                <  7:316779/133‑1 (MQ=255)
                                                                                  cTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAAc                                                            >  6:120594/1‑149 (MQ=255)
                                                                                      gCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAAcgcg                                                         >  4:248499/1‑149 (MQ=255)
                                                                                          tgGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGc                                                       <  2:72860/147‑1 (MQ=255)
                                                                                          tgGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCa                                                     <  6:304443/149‑1 (MQ=255)
                                                                                            gTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGt                                                   >  1:102772/1‑149 (MQ=255)
                                                                                                gTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTc                                               <  7:325273/149‑1 (MQ=255)
                                                                                                gTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTc                                               >  5:318529/1‑149 (MQ=255)
                                                                                                 tttCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCt                                              <  1:205323/149‑1 (MQ=255)
                                                                                                   tCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGTATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCtt                                             <  6:211646/148‑1 (MQ=255)
                                                                                                      cAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTtgtg                                          >  8:251831/1‑149 (MQ=255)
                                                                                                      cAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTtgt                                           <  4:238609/148‑1 (MQ=255)
                                                                                                      cAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTtgt                                           <  2:370543/148‑1 (MQ=255)
                                                                                                         ggttCTGTGT‑‑GATCAGGCTGCCTCTGTCTAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGACATGGATTTTCTATACCGCAGCCTCGTCCATCCACTGCACTTCGATTTCCAACCGGACGGCGCAACGCGGCCAGTTTTCTTGTggg                                        <  6:192021/146‑1 (MQ=255)
                                                                                                                ggT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAAt                                    <  7:46036/146‑1 (MQ=255)
                                                                                                                ggT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGt                                 >  3:337617/1‑149 (MQ=255)
                                                                                                                            cTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCACCCGCCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTg                       <  6:169466/149‑1 (MQ=255)
                                                                                                                                cGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTg                    >  1:326339/1‑148 (MQ=255)
                                                                                                                                    gTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCaa                 >  6:383783/1‑147 (MQ=255)
                                                                                                                                     tCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAAc                <  5:43854/147‑1 (MQ=255)
                                                                                                                                     tCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAg              <  3:248499/149‑1 (MQ=255)
                                                                                                                                     tCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAg              >  5:111960/1‑149 (MQ=255)
                                                                                                                                      ccagcgCATGCTGCGCGAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGt             >  6:208928/6‑149 (MQ=255)
                                                                                                                                       gAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGTg            >  6:319699/1‑149 (MQ=255)
                                                                                                                                       gAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGTg            >  8:135297/1‑149 (MQ=255)
                                                                                                                                        agcgCATGCTGCGCGAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGt             <  6:256840/144‑1 (MQ=255)
                                                                                                                                        agcgCATGCTGCGCGAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGt             >  5:256840/4‑147 (MQ=255)
                                                                                                                                             aTGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGTGGATCgg      >  1:285291/1‑149 (MQ=255)
                                                                                                                                               gCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTgg                                         >  5:34063/1‑112 (MQ=255)
                                                                                                                                               gCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTgg                                         <  6:34063/112‑1 (MQ=255)
                                                                                                                                               gCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGTGGATcatgg    >  8:327834/1‑145 (MQ=255)
                                                                                                                                                   gcgcAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGTGGATCGGCGcc  <  3:63142/147‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
TGCCGCGCACCCAGGCCAGTGCGTCACCTTGCATCCCATAGCAGAAGCGCCGCCAGTAATCGCGCATGACCTGGCCCAACAACTCGCTGTGGTCGGTTTCCAGGGTCTGGGT‑‑GAACAGGCTGCCGCTGTCGAAGGCATGCTCGCGCAGCATGCGCTGGCACCAGCCATGGATGGTCGATACCGCAGCCTCGTCCATCCACTGCACGGCGATTTCCAACCGGCCGGCGCAACGCGGCCAGTTTTCTTGTGGGTAATCGTCACGCAGTTGGTGCAACAGTGGATCGGCGCC  >  NC_002947/5303214‑5303502

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