Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,293 Δ1 bp 9.0% intergenic (+111/+31) 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,2930T. 9.0% 232.0 / 17.4 78intergenic (+111/+31)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base . (3/4);  total (38/40)
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

CTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCG  >  NC_002947/1707146‑1707439
                                                                                                                                                    |                                                                                                                                                    
cTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATc                                                                                                                                                     >  6:294750/1‑149 (MQ=255)
     tGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCggg                                                                                                                                                <  4:36648/149‑1 (MQ=255)
      gggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCgggg                                                                                                                                               >  4:307166/1‑149 (MQ=255)
       ggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGc                                                                                                                                              <  5:242093/149‑1 (MQ=255)
         cGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCt‑‑t                                                                                                                                          <  7:218935/149‑1 (MQ=255)
            ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGCGCGCCCCGCGGGGGGCGCTCGGACTCACAGGGGGCCCCACATTTTCGGCGAAAAAAAAAAAGCCCCCCG‑GGGAGCCCGGGGGT‑‑TGTc                                                                                                                                       >  2:192869/1‑149 (MQ=255)
            ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGCT‑‑TTTc                                                                                                                                       >  1:164234/1‑149 (MQ=255)
            ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTc                                                                                                                                       <  2:49434/149‑1 (MQ=255)
             gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTcc                                                                                                                                      <  1:243573/149‑1 (MQ=255)
             gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTcc                                                                                                                                      >  1:64950/1‑149 (MQ=255)
               aaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGCT‑‑TTTCttt                                                                                                                                    <  5:390519/149‑1 (MQ=38)
                acacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGCT‑‑TTTCtttg                                                                                                                                   <  2:164234/149‑2 (MQ=38)
                  acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCttt                                                                                                                                   >  1:55800/1‑147 (MQ=255)
                    cgcgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAg                                                                                                                                >  1:334519/1‑148 (MQ=255)
                        gCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTag                                                                                                                            >  8:98986/1‑148 (MQ=255)
                        gCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTag                                                                                                                            <  7:98986/148‑1 (MQ=255)
                            gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCt‑‑t                                                                                                                                          >  7:284571/1‑130 (MQ=255)
                            gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCt‑‑t                                                                                                                                          <  8:284571/130‑1 (MQ=255)
                            gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTagccagc                                                                                                                       >  8:211491/1‑149 (MQ=255)
                            gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTagccagc                                                                                                                       >  4:261228/1‑149 (MQ=255)
                                 gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAg                                                                                                                  >  4:157740/1‑149 (MQ=255)
                                  tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAgc                                                                                                                 <  1:159612/149‑1 (MQ=255)
                                  tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑GCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAgc                                                                                                                 >  1:136323/1‑149 (MQ=255)
                                    aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAgcgc                                                                                                               <  6:31107/149‑1 (MQ=255)
                                       tttCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGa                                                                                                            <  4:301158/149‑1 (MQ=255)
                                          cGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCGTTTCAGCTAGCCAGCCTTAGCGCTGAAc                                                                                                          <  5:294750/148‑1 (MQ=255)
                                            tcgcATCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGTCGGCACAACATTATCGGCGAACAAAGAAAAGGCCCTCG‑TCTATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                                                                                       <  8:54112/148‑1 (MQ=39)
                                            gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                                                                                       <  6:317139/149‑1 (MQ=255)
                                             cgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                      <  2:334519/149‑1 (MQ=255)
                                                   aTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAcc                                                                                                <  7:155558/149‑1 (MQ=255)
                                                       ttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGc                                                                                            >  3:50247/1‑149 (MQ=255)
                                                          gcgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtg                                                                                         <  3:261228/149‑1 (MQ=255)
                                                            gcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                                                                       <  7:360840/149‑1 (MQ=255)
                                                               cgcgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGc                                                                                     >  2:9212/1‑148 (MQ=255)
                                                                gcgAGATCGAGCGCCGCGCGGGCGGCGCTCGAGCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCGTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCa                                                                                   >  3:108529/1‑149 (MQ=255)
                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCt                                                                             >  4:243168/1‑149 (MQ=255)
                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGCTCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTgg                                                                           >  7:27044/1‑149 (MQ=255)
                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTgg                                                                           >  7:50460/1‑149 (MQ=255)
                                                                         cgcgcCCCCCGCGCGGCGCTCGATCTCGGGGGGGGCAAAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTCTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg                                                                            <  1:116832/148‑1 (MQ=255)
                                                                         aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt                                                                          >  8:247934/1‑149 (MQ=255)
                                                                         aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt                                                                          >  8:17238/1‑149 (MQ=255)
                                                                          gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                                                                         >  5:22544/1‑149 (MQ=255)
                                                                          gcgcCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                                                                         <  8:50460/149‑1 (MQ=255)
                                                                             ccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTc                                                                      <  5:48953/149‑1 (MQ=255)
                                                                             ccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTc                                                                      >  3:123768/1‑149 (MQ=255)
                                                                             ccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTc                                                                      <  3:307166/149‑1 (MQ=255)
                                                                                  cgggcggCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTgg                                                                 <  3:157740/149‑1 (MQ=255)
                                                                                  cgggcggCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTgg                                                                 <  4:50247/149‑1 (MQ=255)
                                                                                      cggcgCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGAt                                                             <  1:192869/149‑1 (MQ=255)
                                                                                       ggcgCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                            <  6:22544/149‑1 (MQ=255)
                                                                                       ggcgCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                            <  4:123768/149‑1 (MQ=255)
                                                                                          gcTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCttt                                                                                                                                   <  8:292958/75‑1 (MQ=255)
                                                                                          gcTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCttt                                                                                                                                   >  7:292958/1‑75 (MQ=255)
                                                                                           cTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCa                                                        >  4:134955/1‑149 (MQ=255)
                                                                                              gATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCAc                                                       <  4:149809/147‑1 (MQ=255)
                                                                                              gATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCt                                                     >  3:148662/1‑149 (MQ=255)
                                                                                                   cacaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCAcc                                                >  5:4841/1‑149 (MQ=255)
                                                                                                    acgggcggcACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTc                                                    <  4:361012/141‑1 (MQ=255)
                                                                                                    acaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCAccc                                               >  6:394023/1‑149 (MQ=255)
                                                                                                         cggcACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGGGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACgg                                          <  2:136323/149‑1 (MQ=255)
                                                                                                             acaacaTTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCt                                      <  6:4841/149‑1 (MQ=255)
                                                                                                                        ggCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTccc                           >  8:48179/1‑149 (MQ=255)
                                                                                                                           gAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAg                        >  1:170395/1‑149 (MQ=255)
                                                                                                                             aCAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTGCCGAGtt                      >  1:372939/1‑149 (MQ=255)
                                                                                                                             aCAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGtt                      <  7:247934/149‑1 (MQ=255)
                                                                                                                             aCAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGtt                      <  8:27044/149‑1 (MQ=255)
                                                                                                                              cAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTg                     <  3:243168/149‑1 (MQ=255)
                                                                                                                                 aGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTc                  <  5:394023/149‑1 (MQ=255)
                                                                                                                                 aGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTc                  >  3:76042/1‑149 (MQ=255)
                                                                                                                                   acaaGCCCCGCGATCGA‑CGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTc                                                         <  5:262090/106‑1 (MQ=39)
                                                                                                                                   acaaGCCCCGCGATCGA‑CGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTc                                                         >  6:262090/3‑108 (MQ=39)
                                                                                                                                   aaaaGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCTCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCAGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTATGCGCACGTCCCGACTTGcac                  >  3:280417/1‑145 (MQ=255)
                                                                                                                                      aGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGtt             <  2:55800/149‑1 (MQ=255)
                                                                                                                                       gCCCCGCGATCGA‑CGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTg            >  1:261404/1‑149 (MQ=255)
                                                                                                                                        ccccGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGt           >  6:342967/1‑149 (MQ=255)
                                                                                                                                                  gATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCt                                      >  7:216560/1‑113 (MQ=255)
                                                                                                                                                  gATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCt                                      <  8:216560/113‑1 (MQ=255)
                                                                                                                                                  gATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCg  >  1:370129/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCT‑‑TGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCG  >  NC_002947/1707146‑1707439

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