Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,218,661:1 +C 8.8% intergenic (+44/+39) oprJ → / ← nfxB outer membrane protein OprJ/HTH‑type transcriptional regulator NfxB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,218,6611.C 8.8% 256.5 / 17.7 68intergenic (+44/+39)oprJ/nfxBouter membrane protein OprJ/HTH‑type transcriptional regulator NfxB
Reads supporting (aligned to +/- strand):  ref base . (29/33);  new base C (4/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

TACGTGAACGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACACC  >  NC_002947/3218516‑3218808
                                                                                                                                                  |                                                                                                                                                     
tACGTGAACGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATc                                                                                                                                                    >  6:292514/1‑149 (MQ=255)
     gAACGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGc                                                                                                                                               >  5:339000/1‑149 (MQ=255)
       aCGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGccc                                                                                                                                             >  6:78125/1‑149 (MQ=255)
          aaGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTgg                                                                                                                                          <  8:258817/149‑1 (MQ=255)
          aaGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTgg                                                                                                                                          <  1:300175/149‑1 (MQ=255)
            gcagcCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTgggg                                                                                                                                        >  2:115702/1‑149 (MQ=255)
               gcCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACATCCTCCTTTGGTTTGCTTGTCCAGCCACTTGGGCGGTGGGATCGATGCCCTGGGGttt                                                                                                                                     <  6:166857/149‑1 (MQ=255)
               gcCTTCATCGAAGCCAGAACCCAGCGACAACTGGCGTAGGTCAACCTGTTCCGCGCACTCGGCGGCGCCTCCCCGAGCAACGGCTGAAACAGCCTCCATTGGTTGGCTTGTCCAGCCACTTGCCCCCAGGGATCGAAGCCCTgggctat                                                                                                                                     >  1:178628/1‑145 (MQ=255)
                 cTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGttttt                                                                                                                                   <  5:351298/149‑1 (MQ=255)
                   tCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGCATCGATCCCTggg                                                                                                                                         <  5:313750/141‑1 (MQ=255)
                   tCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGCATCGATCCCTggg                                                                                                                                         >  6:313750/1‑141 (MQ=255)
                     aTCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTgg                                                                                                                               <  5:368793/149‑1 (MQ=255)
                     aTCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTgg                                                                                                                               >  2:239127/1‑149 (MQ=255)
                      tCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGCATCGATCCCTggg                                                                                                                                         <  1:357573/138‑1 (MQ=255)
                      tCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGCATCGATCCCTggg                                                                                                                                         >  2:357573/1‑138 (MQ=255)
                        gAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGTTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTg                                                                                                                            <  1:115702/149‑1 (MQ=255)
                        gAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTg                                                                                                                            >  2:138905/1‑149 (MQ=255)
                            ccAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTg                                                                                                                         >  4:225849/1‑148 (MQ=255)
                             cAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGaa                                                                                                                       <  5:292514/149‑1 (MQ=255)
                                   ccAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAgg                                                                                                                >  3:141904/1‑149 (MQ=255)
                                        gACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCagc                                                                                                           <  4:141904/149‑1 (MQ=255)
                                          cAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCagcga                                                                                                         <  4:92254/149‑1 (MQ=255)
                                           aaCTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCagcgag                                                                                                        >  7:91702/1‑149 (MQ=255)
                                                gCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCgg                                                                                                     <  6:304029/147‑1 (MQ=255)
                                                  gTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGc                                                                                                 >  2:111499/1‑149 (MQ=255)
                                                     ggTCAACCGGTTCCGCGCACTCGGCTGCGGCTGGCCGAGCACGGGCTGAAACAGCCTCCGTTGGTTGGCTTGTCCAGCCACGTGGCCCCAGGGATCGATGCCCTGGGGATTTTTTGGGTGGTGAAGGCAGGTCAGCGACCGGGTGgcgc                                                                                              >  6:202011/1‑149 (MQ=255)
                                                         aaCCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCgcagca                                                                                          >  1:42634/1‑149 (MQ=255)
                                                            cTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTgg                                                                                                                                          >  2:126909/1‑99 (MQ=255)
                                                            cTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTgg                                                                                                                                          <  1:126909/99‑1 (MQ=255)
                                                            cTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCg                                                                                       >  1:364621/1‑149 (MQ=255)
                                                                 ccGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGa                                                                                   <  6:351827/148‑1 (MQ=255)
                                                                  cgcgCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAg                                                                                  <  1:239127/149‑1 (MQ=255)
                                                                     gcACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAg                                                                                  <  5:143082/145‑1 (MQ=255)
                                                                     gcACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAg                                                                                  >  6:143082/1‑145 (MQ=255)
                                                                        cTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAagc                                                                          <  2:364621/149‑1 (MQ=255)
                                                                           ggcggcggcTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAagca                                                                         <  3:225849/147‑1 (MQ=255)
                                                                               gcggcTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCa                                                                   <  5:137941/149‑1 (MQ=255)
                                                                                     ggCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGttt                                                             <  2:176017/149‑1 (MQ=255)
                                                                                      gCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTg                                                            >  1:364780/1‑149 (MQ=255)
                                                                                      gCCGAGCAAGGGCTGAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTg                                                                                                  <  6:199050/111‑1 (MQ=255)
                                                                                      gCCGAGCAAGGGCTGAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTg                                                                                                  >  5:199050/1‑111 (MQ=255)
                                                                                        cGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCa                                                          >  5:153032/1‑149 (MQ=255)
                                                                                             aaGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCATGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACAGAc                                                     <  6:339000/149‑1 (MQ=255)
                                                                                             aaGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGAc                                                     <  8:91702/149‑1 (MQ=255)
                                                                                             aaGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGAc                                                     <  1:138905/149‑1 (MQ=255)
                                                                                                ggCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTc                                                  <  5:78220/149‑1 (MQ=255)
                                                                                                   tGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTcgcg                                               >  7:331166/1‑149 (MQ=255)
                                                                                                     aaaCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCg                                                                                               >  1:139675/1‑100 (MQ=255)
                                                                                                     aaaCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCg                                                                                               <  2:139675/100‑1 (MQ=255)
                                                                                                     aaaCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTcgcg                                               <  7:107285/147‑1 (MQ=255)
                                                                                                      aaCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTg                                            >  6:362046/1‑149 (MQ=255)
                                                                                                      aaCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTg                                            >  2:308565/1‑149 (MQ=255)
                                                                                                         aGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGTCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTc                                                  >  5:104397/1‑140 (MQ=255)
                                                                                                         aGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGTCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTc                                                  <  6:104397/140‑1 (MQ=255)
                                                                                                                        gCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGcgtcgt                          >  5:176651/1‑149 (MQ=255)
                                                                                                                          ttGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTc                        >  3:126057/1‑149 (MQ=255)
                                                                                                                               cAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCAt                   >  2:300358/1‑149 (MQ=255)
                                                                                                                                   cACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATcgac               <  5:251209/149‑1 (MQ=255)
                                                                                                                                        ggCCCCAGGGATCTTTGCCCTGGGGTTTTTTTGGCTGGTTAAGGCAGGGCAGCTAGCGGGTGCCGCAGCACCGTGGAGTATCTAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATcgacga             <  1:111499/147‑1 (MQ=255)
                                                                                                                                        ggCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGAt            <  3:196949/147‑1 (MQ=255)
                                                                                                                                        ggCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGAt            <  1:203555/147‑1 (MQ=255)
                                                                                                                                         gCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGc          <  6:153032/148‑1 (MQ=255)
                                                                                                                                          ccccAGGGCATCGATCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCg        >  6:237121/1‑149 (MQ=255)
                                                                                                                                           cccAGGGCATCGATCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCg        >  2:51950/1‑148 (MQ=255)
                                                                                                                                           cccAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATg           >  6:384287/1‑145 (MQ=255)
                                                                                                                                           cccAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATg           <  5:384287/145‑1 (MQ=255)
                                                                                                                                              aGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTaca    >  4:38425/1‑149 (MQ=255)
                                                                                                                                              aGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTaca    <  5:78125/149‑1 (MQ=255)
                                                                                                                                                ggATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTGAGCGAGCGGGTGCCGCAGCACGGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGGATCGACGATGCCGTACAcc  <  2:178628/149‑1 (MQ=255)
                                                                                                                                                ggATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACAcc  <  5:362046/149‑1 (MQ=255)
                                                                                                                                                ggATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGACCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACAcc  <  1:308565/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TACGTGAACGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACACC  >  NC_002947/3218516‑3218808

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: