Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,272 G→C 6.7% intergenic (+48/+94) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,2720GC 6.7% 230.8 / 5.1 74intergenic (+48/+94)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base C (2/3);  total (35/39)
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.69e-01

TGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGA  >  NC_002947/1707124‑1707421
                                                                                                                                                    |                                                                                                                                                      
tgTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcggg                                                                                                                                                       <  1:1044485/150‑1 (MQ=255)
 gTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggc                                                                                                                                                      >  1:240551/1‑150 (MQ=255)
  ttCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggcg                                                                                                                                                     >  2:1120720/1‑150 (MQ=255)
      acacCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                                 >  2:223214/1‑150 (MQ=255)
        acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               <  2:1373559/150‑1 (MQ=255)
        acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               <  1:93684/150‑1 (MQ=255)
        acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                               <  2:20656/150‑1 (MQ=255)
           cagcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               <  2:771009/147‑1 (MQ=255)
           cagcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               >  1:771009/1‑147 (MQ=255)
             gcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                          >  2:369062/1‑150 (MQ=255)
                gAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                       >  2:834034/1‑150 (MQ=255)
                     gCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcg                                                                                                                                  >  2:10550/1‑150 (MQ=255)
                       cTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAggcggc                                                                                                                                >  1:653817/1‑150 (MQ=255)
                         gACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGcac                                                                                                                              >  1:982672/1‑150 (MQ=255)
                             ggCGTGCTGGTGGACAAGCGCAAGCGCGGGGTGAGGTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCGCGCTCTCACAGGCGGCacaaca                                                                                                                          <  2:208168/150‑1 (MQ=39)
                                  gCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATc                                                                                                                     <  1:355818/150‑1 (MQ=255)
                                   cTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCg                                                                                                                    <  1:73739/150‑1 (MQ=255)
                                     ggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                                  <  2:1436702/150‑1 (MQ=255)
                                            aaCCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaaa                                                                                                           >  2:253670/1‑150 (MQ=255)
                                              ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGa                                                                                                         >  2:892996/1‑150 (MQ=255)
                                                   ggcgcgcGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTAGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGc                                                                                                    <  2:240551/149‑1 (MQ=255)
                                                           gTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑Tc                                                                                           <  1:369062/150‑1 (MQ=255)
                                                           gTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑Tc                                                                                           <  1:223214/150‑1 (MQ=255)
                                                             gcgcTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGa                                                                                         <  2:653817/150‑1 (MQ=255)
                                                                   cATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCgg                                                                                   <  1:1419186/150‑1 (MQ=255)
                                                                      aggCGAGGGAACGCCGCGGCAGGTGCAGTGAGGTTTGGGGGCAGCATTTTTGGCGCCAGCGAGATGGAGGGCCGCGCGGGCGGCGCTCGCTCTCACAGGGGGCACAACATTATCGGCGAAAAAAGAAAAGCCCCGCG‑TCGATCGCtgggg                                                                                >  2:497033/2‑149 (MQ=255)
                                                                        gCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCtt                                                                              <  2:961708/150‑1 (MQ=255)
                                                                          gAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGCtttt                                                                            <  1:892996/150‑3 (MQ=255)
                                                                            ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTcc                                                                          >  1:180238/1‑150 (MQ=255)
                                                                             gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCt                                                                         >  1:801770/1‑150 (MQ=255)
                                                                             gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCt                                                                         <  1:401580/150‑1 (MQ=255)
                                                                                 cacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCa                                                                     <  1:253670/150‑1 (MQ=255)
                                                                                  acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAg                                                                    >  1:564246/1‑150 (MQ=255)
                                                                                           ggTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                            >  1:1357664/1‑70 (MQ=255)
                                                                                           ggTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                            <  2:1357664/70‑1 (MQ=255)
                                                                                                aGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAg                                                      >  1:477483/1‑150 (MQ=255)
                                                                                                       tttCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGc                                                                   <  1:1418435/130‑1 (MQ=255)
                                                                                                       tttCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGc                                                                   >  2:1418435/1‑130 (MQ=255)
                                                                                                       tttCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATAGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGaa                                               >  1:1263731/1‑150 (MQ=255)
                                                                                                           ggCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCgg                                                                                   <  1:534192/110‑1 (MQ=255)
                                                                                                           ggCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCgg                                                                                   >  2:534192/1‑110 (MQ=255)
                                                                                                            gcgcAGCATTTTTGGCGCCCGCGCGATCGAGGGGCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                          <  2:146911/150‑1 (MQ=255)
                                                                                                            gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                          <  2:180238/150‑1 (MQ=255)
                                                                                                            gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                          <  2:982672/150‑1 (MQ=255)
                                                                                                              gcagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTa                                                                 >  2:1174599/1‑125 (MQ=255)
                                                                                                              gcagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTa                                                                 <  1:1174599/125‑1 (MQ=255)
                                                                                                              gcagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCt                                        >  1:419036/1‑150 (MQ=255)
                                                                                                               cagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                        >  1:1503141/1‑54 (MQ=255)
                                                                                                               cagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                        <  2:1503141/54‑1 (MQ=255)
                                                                                                               cagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTc                                       >  2:1333287/1‑150 (MQ=255)
                                                                                                                agcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAAc                                                                                                              >  1:334164/1‑79 (MQ=255)
                                                                                                                agcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAAc                                                                                                              <  2:334164/79‑1 (MQ=255)
                                                                                                                   aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCg                                   >  1:48500/1‑150 (MQ=255)
                                                                                                                   aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCg                                   >  1:495751/1‑150 (MQ=255)
                                                                                                                         ggCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtg                             <  2:564246/150‑1 (MQ=255)
                                                                                                                         ggCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtg                             <  2:801770/150‑1 (MQ=255)
                                                                                                                          gcgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑t                                                                                            >  1:1237524/1‑86 (MQ=255)
                                                                                                                          gcgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑t                                                                                            <  2:1237524/86‑1 (MQ=255)
                                                                                                                           cgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTagccagcc                                                          <  1:1330467/119‑1 (MQ=255)
                                                                                                                           cgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTagccagcc                                                          >  2:1330467/1‑119 (MQ=255)
                                                                                                                                    gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGa                                  <  2:871709/134‑1 (MQ=255)
                                                                                                                                    gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGa                                  >  1:871709/1‑134 (MQ=255)
                                                                                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg                <  2:109783/150‑1 (MQ=255)
                                                                                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg                >  1:109783/1‑150 (MQ=255)
                                                                                                                                       cgggcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTgg               <  1:46554/147‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCgggg                                                                                 >  2:100590/1‑83 (MQ=255)
                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCgggg                                                                                 <  1:100590/83‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt              <  2:1263731/150‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt              <  2:477483/150‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt              >  2:1258849/1‑150 (MQ=255)
                                                                                                                                         aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt             <  1:802329/150‑1 (MQ=255)
                                                                                                                                         aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt             >  1:1030958/1‑150 (MQ=255)
                                                                                                                                          gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGa                                  >  2:811964/1‑128 (MQ=255)
                                                                                                                                          gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGa                                  <  1:811964/128‑1 (MQ=255)
                                                                                                                                           cgccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGt           >  1:411547/1‑150 (MQ=255)
                                                                                                                                              cgcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGt        >  1:924933/1‑150 (MQ=255)
                                                                                                                                               gcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt       <  1:834034/150‑1 (MQ=255)
                                                                                                                                                    ggcggcGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGgaga  <  1:1258849/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGA  >  NC_002947/1707124‑1707421

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