Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,518,908 C→T 6.8% M1V (GTG→GTA) † PP_4851 ← phosphate starvation‑inducible protein PsiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,518,9080CT 6.8% 173.0 / 2.8 59M1V (GTG→GTA) †PP_4851phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base T (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGCCAAGGTCGAAGGCAATGGCC  >  NC_002947/5518771‑5519055
                                                                                                                                           |                                                                                                                                                   
tCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGa                                                                                                                                          <  4:130668/149‑1 (MQ=255)
  ggCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGc                                                                                                                                        <  1:133161/149‑1 (MQ=255)
  ggCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGc                                                                                                                                        <  7:213084/149‑1 (MQ=255)
   gCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCt                                                                                                                                       >  4:14155/1‑149 (MQ=255)
     gTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGa                                                                                                                                     <  5:81125/149‑1 (MQ=255)
       tGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAgt                                                                                                                                   >  4:31577/1‑149 (MQ=255)
        gCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAgtg                                                                                                                                  >  3:301814/1‑149 (MQ=255)
        gCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAgt                                                                                                                                   <  7:208/148‑1 (MQ=255)
          aGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAgtg                                                                                                                                  >  4:239624/1‑147 (MQ=255)
          aGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAgtg                                                                                                                                  >  3:229293/1‑147 (MQ=255)
                 cATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGTACGGTGATCGAGCTGAGTGTAGCCAAg                                                                                                                          >  4:36209/1‑148 (MQ=255)
                 cATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAg                                                                                                                          >  3:132673/1‑148 (MQ=255)
                       ttCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAAt                                                                                                                       <  4:132673/145‑1 (MQ=255)
                       ttCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGtt                                                                                                                    <  2:54658/148‑1 (MQ=255)
                        tCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTc                                                                                                                   <  4:153188/148‑1 (MQ=255)
                          ctgctgCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATACACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCa                                                                                                                  <  8:77712/147‑1 (MQ=255)
                          ctgctgCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATACACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCa                                                                                                                  >  7:77712/1‑147 (MQ=255)
                                          gCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCAGCGCCCGATGTACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGAGCCCTGTCCGCGCCTg                                                                                                >  7:127054/1‑149 (MQ=255)
                                           cAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGc                                                                                               <  3:14155/149‑1 (MQ=255)
                                              cGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCcgg                                                                                            <  4:8006/149‑1 (MQ=255)
                                               gAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCcggc                                                                                           <  4:301814/149‑1 (MQ=255)
                                                aaaCCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCcggcg                                                                                          <  6:44835/149‑1 (MQ=255)
                                                aaaCCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCcggcg                                                                                          <  3:31577/149‑1 (MQ=255)
                                                  aCCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTCTGAGTGACTACCCGCCCGATACACGGTGATAGAGCTGAGTGTATCCAAGAATGTTCAGCGCCCCGTCAGCGCCAGCcggcg                                                                                          >  1:264760/1‑149 (MQ=255)
                                                   ccTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTa                                                                                       >  6:328470/1‑149 (MQ=255)
                                                        gtctcTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGAACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGTACGGTGATCGAGCTGAGTGTAGGCAAGAATGTTCAGCGCCCTGTCCGCGCCTGCCGGCGGTACTGc                                                                                   <  8:127054/144‑1 (MQ=255)
                                                           aCTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCgg                                                                               >  3:56327/1‑149 (MQ=255)
                                                                acaacaCCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTa                                                                          <  7:286656/149‑1 (MQ=255)
                                                                  aacaCCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTaca                                                                        >  2:283642/1‑149 (MQ=255)
                                                                  aacaCCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTaca                                                                        >  5:217974/1‑149 (MQ=255)
                                                                     acCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTc                                                                                 >  3:72897/1‑137 (MQ=255)
                                                                     acCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTc                                                                                 <  4:72897/137‑1 (MQ=255)
                                                                      ccGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGcc                                                                                              <  2:313989/123‑1 (MQ=255)
                                                                      ccGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGcc                                                                                              >  1:313989/1‑123 (MQ=255)
                                                                          cagcgccagTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCt                                                                >  5:218643/1‑149 (MQ=255)
                                                                             cgccagTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGAACGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGgc                                                             <  1:278043/149‑1 (MQ=255)
                                                                                cagTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGTACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCt                                                          <  3:36209/149‑1 (MQ=255)
                                                                                     aCGGGATTTGCAGCACTTTCTTGGAGTGGTCTC‑‑CTTGTCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTATCCAAGAATGTTCAGGGCCCCTTCCGCGCCTGCCGTCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGaa                                                      <  2:264760/148‑1 (MQ=255)
                                                                                     aCGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGaa                                                      <  1:283642/148‑1 (MQ=255)
                                                                                     aCGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGaa                                                      >  3:108080/1‑148 (MQ=255)
                                                                                     aCGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGaa                                                      >  4:281359/1‑148 (MQ=255)
                                                                                      cGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACa                                                    >  2:229720/1‑149 (MQ=255)
                                                                                              gcagcaCTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGAt                                             <  1:115658/148‑1 (MQ=255)
                                                                                                     tttCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGcc                                      >  6:308114/1‑148 (MQ=255)
                                                                                                        cATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGaa                                   <  6:217974/148‑1 (MQ=255)
                                                                                                          tGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTg                                                         >  1:297635/1‑124 (MQ=255)
                                                                                                          tGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTg                                                         <  2:297635/124‑1 (MQ=255)
                                                                                                          tGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAgcg                                >  3:117878/1‑149 (MQ=255)
                                                                                                             aGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCgctg                             <  4:56327/149‑1 (MQ=255)
                                                                                                                 gTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATg                                                  <  4:195568/124‑1 (MQ=255)
                                                                                                                    tc‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCGAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGc                      >  7:172760/1‑149 (MQ=255)
                                                                                                                        cTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGcc                     <  6:218643/148‑1 (MQ=255)
                                                                                                                            gCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGCCAAGGt                >  5:5899/1‑149 (MQ=255)
                                                                                                                             cTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCgctg                             >  3:123535/1‑135 (MQ=255)
                                                                                                                             cTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCgctg                             <  4:123535/135‑1 (MQ=255)
                                                                                                                                 gcgcCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGCCAAGGTCGAAg           <  4:108080/149‑1 (MQ=255)
                                                                                                                                        aTGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGCCGATGAACGCCGAAGCGCTGCTGTAGCCAAGGTCGAAGGCAATgg    <  1:229720/149‑1 (MQ=255)
                                                                                                                                         tGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGCCAAGGTCGAAGGCAATGGc   <  4:117878/149‑1 (MQ=255)
                                                                                                                                          gCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGCCAAGGTCGAAGGCAATGGcc  >  4:116630/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
TCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTGGTCTC‑‑CTTGGCTCCGCGCCCGATGCACGGTGATCGAGCTGAGTGTAGCCAAGAATGTTCAGCGCCCCGTCCGCGCCTGCCGGCGGTACTGCTCGGGCGTACACCGGGCCTGGCGCTGGAACATGGCGATGAACGCCGAAGCGCTGCTGTAGCCAAGGTCGAAGGCAATGGCC  >  NC_002947/5518771‑5519055

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