Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,123,837 T→C 6.0% intergenic (‑37/+184) PP_mr13 ← / ← gcvT‑I cspA regulatory region/aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,123,8370TC 6.0% 192.3 / 2.1 66intergenic (‑37/+184)PP_mr13/gcvT‑IcspA regulatory region/aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base C (2/2);  total (31/35)
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

TGATGACAGTTACTTTCGCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAA  >  NC_002947/1123692‑1123984
                                                                                                                                                 |                                                                                                                                                   
tgatgaCAGTTACTTTCGCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGa                                                                                                                                                   <  4:99378/148‑1 (MQ=255)
   tgaCAGTTACTTTCGCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTATAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGActc                                                                                                                                                <  3:156070/148‑1 (MQ=255)
       acttACTTTCGCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGc                                                                                                                                           <  5:102577/147‑1 (MQ=255)
       aGTTACTTTCGCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGc                                                                                                                                           >  5:39846/1‑149 (MQ=255)
         ttACTTTCGCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCAGGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGcc                                                                                                                                          >  1:224558/1‑148 (MQ=255)
                 gCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGa                                                                                                                                 <  5:74585/149‑1 (MQ=255)
                 gCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTGGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCg                                                                                                                                  <  5:40972/148‑1 (MQ=255)
                     ttttttGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCg                                                                                                                              >  1:280511/1‑148 (MQ=255)
                            tAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAAc                                                                                                                        >  8:121422/1‑147 (MQ=255)
                             aaGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAg                                                                                                                     >  2:89328/1‑149 (MQ=255)
                              aGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAg                                                                                                                     <  2:224558/148‑1 (MQ=255)
                               gCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGcc                                                                                                                   >  2:14951/1‑149 (MQ=255)
                                cGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGa                                                                                                                                 <  2:104057/134‑1 (MQ=255)
                                cGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGa                                                                                                                                 >  1:104057/1‑134 (MQ=255)
                                cGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGccc                                                                                                                  >  6:59644/1‑149 (MQ=255)
                                  aaaGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCtt                                                                                                                >  2:71526/1‑149 (MQ=255)
                                     gATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCACAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCATGGt                                                                                                             <  5:105534/149‑1 (MQ=255)
                                      aTTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAg                                                                                                                                <  4:211630/129‑1 (MQ=255)
                                      aTTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAg                                                                                                                                >  3:211630/1‑129 (MQ=255)
                                          ctgggctgggTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTAt                                                                                                        <  7:271252/149‑1 (MQ=255)
                                          ctgggctgggTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTAt                                                                                                        <  3:146287/149‑1 (MQ=255)
                                           tgggctgggTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATg                                                                                                       >  2:10791/1‑149 (MQ=255)
                                             ggctgggTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATgg                                                                                                      >  2:251389/1‑148 (MQ=255)
                                             ggctgggTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATg                                                                                                       <  2:280511/147‑1 (MQ=255)
                                             ggctgggTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGa                                                                                                     <  1:10791/149‑1 (MQ=255)
                                                tgggTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAATACAAGCCCTTGGTCGTATGGaaaa                                                                                                  <  8:52684/149‑1 (MQ=255)
                                                    ttGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTgg                                                                                              >  5:158281/1‑149 (MQ=255)
                                                       cAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCg                                                                                            <  1:14951/148‑1 (MQ=255)
                                                        aGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGt                                                                                          >  8:261969/1‑149 (MQ=255)
                                                        aGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGt                                                                                          >  7:46871/1‑149 (MQ=255)
                                                                cgagaCGTCCATCTGGTGGTCGCAAATAGCGGCTACTGGCCCAGGGCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGAt                                                                                  <  8:103444/148‑1 (MQ=255)
                                                                 gagaCGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGAt                                                                                  >  3:227183/1‑148 (MQ=255)
                                                                 gagaCGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGAt                                                                                  >  8:160147/1‑148 (MQ=255)
                                                                  agaCGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATaa                                                                                <  5:59644/149‑1 (MQ=255)
                                                                         gAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAaccgacc                                                                          >  2:305063/1‑148 (MQ=255)
                                                                          aaCGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAATACTGGCGGTTTGCGGATAAccgaccga                                                                        >  2:248835/1‑149 (MQ=255)
                                                                            cGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATaa                                                                                <  8:60584/139‑1 (MQ=255)
                                                                            cGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATaa                                                                                >  7:60584/1‑139 (MQ=255)
                                                                               gTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCTAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGt                                                                   <  1:251389/149‑1 (MQ=255)
                                                                                ttGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTc                                                                  <  7:121422/149‑1 (MQ=255)
                                                                                    aGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGt                                                                                                          <  2:153316/105‑1 (MQ=255)
                                                                                    aGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGt                                                                                                          >  1:153316/1‑105 (MQ=255)
                                                                                     gCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGaaa                                                             <  7:261969/149‑1 (MQ=255)
                                                                                       aaaTCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAg                                                                                                                     >  8:231542/1‑91 (MQ=255)
                                                                                       aaaTCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAg                                                                                                                     <  7:231542/91‑1 (MQ=255)
                                                                                       aaaTCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGAt                                                                       >  5:225121/1‑137 (MQ=255)
                                                                                       aaaTCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGAt                                                                       <  6:225121/137‑1 (MQ=255)
                                                                                       aaaTCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACg                                                           <  1:71526/149‑1 (MQ=255)
                                                                                       aaaTCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACg                                                           <  6:158281/149‑1 (MQ=255)
                                                                                         aTCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGc                                                         <  4:227183/149‑1 (MQ=255)
                                                                                              ggCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCAt                                                       <  1:305063/146‑1 (MQ=255)
                                                                                                    tGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGCGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctgg                                              >  6:89774/1‑149 (MQ=255)
                                                                                                    tGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGCGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCgcctgg                                              >  6:85473/1‑149 (MQ=255)
                                                                                                                     tGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCa                              >  8:176977/1‑148 (MQ=255)
                                                                                                                       ttGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaa                             >  4:46846/1‑147 (MQ=255)
                                                                                                                       ttGCGACCCATGCAAACACAGTGCAGCGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaa                            <  5:85473/148‑1 (MQ=255)
                                                                                                                       ttGCGACCCATGCAAACACAGTGCAGCGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaa                            <  5:89774/148‑1 (MQ=255)
                                                                                                                         gCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCctggcctggc                                        >  3:239418/1‑134 (MQ=255)
                                                                                                                         gCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCctggcctggc                                        <  4:239418/134‑1 (MQ=255)
                                                                                                                         gCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                           >  7:149099/1‑147 (MQ=255)
                                                                                                                         gCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCaaaa                           <  8:149099/147‑1 (MQ=255)
                                                                                                                         gCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAg                          <  7:176977/148‑1 (MQ=255)
                                                                                                                           gACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCGGGCGCCATGTTCAAAAGTTa                       <  8:46871/149‑1 (MQ=255)
                                                                                                                            aCCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGATTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAt                      <  3:246458/149‑1 (MQ=255)
                                                                                                                               cATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg                   >  2:47397/1‑149 (MQ=255)
                                                                                                                                        acaGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTc           >  2:24724/1‑148 (MQ=255)
                                                                                                                                                gtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGaa  <  7:160147/149‑1 (MQ=255)
                                                                                                                                                gtgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTACAAAAGTTATTGGGCAGGTTCGAAAACGaa  <  1:24724/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGATGACAGTTACTTTCGCACTTTTTTGTAAGCGAAAGATTACTGGGCTGGGTTGCAGGAAGTAAGAGACGTCGAACGGGTTGTAGCAAATCTCGGCTACTGGCCCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAA  >  NC_002947/1123692‑1123984

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: ‑