Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,944,533 T→G 6.3% intergenic (+34/+228) actP‑I → / ← PP_1744 acetate permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,944,5330TG 6.3% 191.6 / 2.5 62intergenic (+34/+228)actP‑I/PP_1744acetate permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base G (2/2);  total (33/29)
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.65e-01

GGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCA  >  NC_002947/1944394‑1944673
                                                                                                                                           |                                                                                                                                            
ggTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCAGCTTGTCGTGGCGtt                                                                                                                                     >  7:224782/1‑149 (MQ=255)
 gTCACCGACATGTGCAAGTCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGAGTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCATTCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGttt                                                                                                                                    <  8:224782/149‑1 (MQ=255)
  tCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGACAAGCGGGGCGtttt                                                                                                                                   <  5:312588/149‑1 (MQ=255)
      cGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTAt                                                                                                                               >  5:112761/1‑149 (MQ=255)
         cAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTc                                                                                                                            <  6:112761/149‑1 (MQ=255)
             tCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATg                                                                                                                         >  6:226535/1‑148 (MQ=255)
               cAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCAt                                                                                                                          >  8:84192/1‑145 (MQ=255)
               cAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCAt                                                                                                                          <  7:84192/145‑1 (MQ=255)
                aaGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGt                                                                                                                     >  3:312715/1‑149 (MQ=255)
                      gcTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGa                                                                                                               >  6:88107/1‑149 (MQ=255)
                        tGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                             >  1:12893/1‑149 (MQ=255)
                          aGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                             <  4:74220/147‑1 (MQ=255)
                               gAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTc                                                                                                        >  7:302262/1‑147 (MQ=255)
                                   gcgcgcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGa                                                                                                    >  2:239524/1‑147 (MQ=255)
                                       gcTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGtt                                                                                              <  5:226535/149‑1 (MQ=255)
                                        cTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTg                                                                                             <  1:239524/149‑1 (MQ=255)
                                        cTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTg                                                                                             <  4:312715/149‑1 (MQ=255)
                                           ttcttcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTg                                                                                          <  2:287008/149‑1 (MQ=255)
                                           ttcttcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTg                                                                                          <  5:88107/149‑1 (MQ=255)
                                               tcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGa                                                                                                    >  5:172106/1‑135 (MQ=255)
                                               tcCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGa                                                                                                    <  6:172106/135‑1 (MQ=255)
                                                   gCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGctct                                                                                      >  5:178642/1‑145 (MQ=255)
                                                   gCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGctct                                                                                      <  6:178642/145‑1 (MQ=255)
                                                          gTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                             >  7:249688/1‑115 (MQ=255)
                                                          gTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATg                                                                                                             <  8:249688/115‑1 (MQ=255)
                                                                 cGCAGACTGGCCTGGGGGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCAt                                                                     <  8:277591/148‑1 (MQ=255)
                                                                        tGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAg                                                             >  6:304590/1‑149 (MQ=255)
                                                                                   ccAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATgg                                                  >  5:109677/1‑149 (MQ=255)
                                                                                     aGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCa                                                                                                                       <  8:73096/78‑1 (MQ=255)
                                                                                     aGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCa                                                                                                                       >  7:73096/1‑78 (MQ=255)
                                                                                     aGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATggag                                                >  7:224728/1‑149 (MQ=255)
                                                                                          ggCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACgg                                           >  4:183870/1‑149 (MQ=255)
                                                                                             ggTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGACGACGGca                                         >  6:211091/1‑148 (MQ=255)
                                                                                                   gCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTgg                                  >  5:222937/1‑149 (MQ=255)
                                                                                                     aCTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGc                                                                                                           <  8:72246/74‑1 (MQ=255)
                                                                                                     aCTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGc                                                                                                           >  7:72246/1‑74 (MQ=255)
                                                                                                       tAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGt                                                                <  4:257004/115‑1 (MQ=255)
                                                                                                       tAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGt                                                                >  3:257004/1‑115 (MQ=255)
                                                                                                         aCATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATg                                                   <  7:83255/126‑1 (MQ=255)
                                                                                                         aCATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATg                                                   >  8:83255/1‑126 (MQ=255)
                                                                                                         aCATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTg                            >  5:337028/1‑149 (MQ=255)
                                                                                                           aTCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGt                          >  4:75546/1‑149 (MQ=255)
                                                                                                            tCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGt                          >  6:69626/1‑148 (MQ=255)
                                                                                                            tCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTc                         <  6:109677/149‑1 (MQ=255)
                                                                                                              aaTCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCCCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCaa                       >  8:292462/1‑149 (MQ=255)
                                                                                                               aTCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAAt                      <  3:183870/149‑1 (MQ=255)
                                                                                                               aTCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAAt                      >  1:44738/1‑149 (MQ=255)
                                                                                                                tCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAAt                      >  8:35224/1‑148 (MQ=255)
                                                                                                                  cccTGCATTGAAAACGCCCGGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGaa                   <  5:304590/149‑1 (MQ=255)
                                                                                                                  cccTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGaa                   <  5:126978/149‑1 (MQ=255)
                                                                                                                    cTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGc                 >  1:275944/1‑149 (MQ=255)
                                                                                                                       cATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCCCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCa               <  7:292462/148‑1 (MQ=255)
                                                                                                                        ataaaaaaCGCCCCGACAAGCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAAt                      >  7:169820/5‑140 (MQ=255)
                                                                                                                        ataaaaaaCGCCCCGACAAGCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAAt                      <  8:169820/136‑1 (MQ=255)
                                                                                                                            aaaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGt         <  5:69626/149‑1 (MQ=255)
                                                                                                                            aaaaCGCCCCGACAAGCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGt         >  2:189944/1‑149 (MQ=255)
                                                                                                                             aaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGtt        <  7:35224/149‑1 (MQ=255)
                                                                                                                             aaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGtt        <  3:75546/149‑1 (MQ=255)
                                                                                                                             aaaCGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGt         <  2:275944/148‑1 (MQ=255)
                                                                                                                                 gCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAAt    <  6:222937/149‑1 (MQ=255)
                                                                                                                                   cccGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCa  >  1:133576/1‑149 (MQ=255)
                                                                                                                                   cccGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCa  >  8:231259/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GGTCACCGACAAGTCCAAGGCGGCTGAGGATGAGCGCGCGCTGTTCTTCCCGCAGTTCGTGCGTTCGCAGACTGGCCTGGGTGCCAGCGGGGCGGTCTCGCACTAACATCAATCCCCTGCATTGAAAACGCCCCGCTTGTCGGGGCGTTTTTTATTTCATGCAGTCGGCGATGGCCTCCAGACAGGTTGGTGCTCTGCCCATGCTCCCCTCATGTCGTGAGCGCCTTCATGGAGGACGGCACATTTGGGGGCTGGTCAATGAAGCCATACTGTTTAATCA  >  NC_002947/1944394‑1944673

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