Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 593,434 A→C 9.1% E391A (GAG→GCG) ‡ PP_0503 → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947593,4340AC 9.1% 178.5 / 7.3 64E391A (GAG→GCG) ‡PP_0503MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base C (4/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCG  >  NC_002947/593300‑593573
                                                                                                                                       |                                                                                                                                           
cAGCGGTTTGGGGGTGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAAt                                                                                                                               >  2:248669/1‑149 (MQ=255)
     gTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCa                                                                                                                          <  5:138150/149‑1 (MQ=255)
      tttGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                                                         >  7:15296/1‑149 (MQ=255)
      tttGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                                                         <  5:284256/149‑1 (MQ=255)
          ggggCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTc                                                                                                                      <  3:94878/148‑1 (MQ=255)
            ggCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAg                                                                                                                    >  6:143412/1‑148 (MQ=255)
             gCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGTGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCGGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGGCGGCCTAGGCCGCCCTGCAATGTTCAg                                                                                                                         <  1:205876/143‑1 (MQ=255)
             gCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTCGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGCGGCCGCCCTGCAATGTTCAg                                                                                                                         >  2:205876/1‑142 (MQ=255)
                   ttGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTg                                                                                                              >  4:235221/1‑147 (MQ=255)
                   ttGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                                                            <  3:112174/149‑1 (MQ=255)
                           cTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgac                                                                                                    <  1:248669/149‑1 (MQ=255)
                              ggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccga                                                                                                 <  5:143412/149‑1 (MQ=255)
                                 gtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCa                                                                                              >  8:346910/1‑149 (MQ=255)
                                     gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGtt                                                                                                                            <  2:316951/115‑1 (MQ=255)
                                     gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGtt                                                                                                                            >  1:316951/1‑115 (MQ=255)
                                     gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                           <  8:249876/132‑1 (MQ=255)
                                     gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                           >  7:249876/1‑132 (MQ=255)
                                     gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                           <  6:30893/132‑1 (MQ=255)
                                     gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGAGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                           >  5:30893/1‑132 (MQ=255)
                                       acacTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGtata                                                                                        >  3:268570/1‑149 (MQ=255)
                                                        cGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTg                                                                                                                                   >  6:272300/1‑89 (MQ=255)
                                                        cGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTg                                                                                                                                   <  5:272300/89‑1 (MQ=255)
                                                          gtgtgtGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgac                                                                     >  3:149549/1‑149 (MQ=255)
                                                             tgtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgt                                                                   <  2:214005/148‑1 (MQ=255)
                                                             tgtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgt                                                                   >  1:214005/1‑148 (MQ=255)
                                                              gtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtc                                                                  >  7:301583/1‑148 (MQ=255)
                                                                      ggCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTg                                                           <  3:235221/147‑1 (MQ=255)
                                                                      ggCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGt                                                         >  7:321563/1‑149 (MQ=255)
                                                                              gcacAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGt                                                 >  5:323580/4‑149 (MQ=255)
                                                                              cATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTg                                                  >  2:44530/1‑148 (MQ=255)
                                                                                  aGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                                                            <  7:190632/86‑1 (MQ=255)
                                                                                  aGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                                                            >  8:190632/1‑86 (MQ=255)
                                                                                   gTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgac                                                                           >  3:216413/1‑118 (MQ=255)
                                                                                   gTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgac                                                                           <  4:216413/118‑1 (MQ=255)
                                                                                    tGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCaa                                           <  7:346910/149‑1 (MQ=255)
                                                                                     gccgccGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAg                                          >  3:260970/1‑149 (MQ=255)
                                                                                     gccgccGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAg                                          >  5:139394/1‑149 (MQ=255)
                                                                                            cACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTc                                   <  4:268570/149‑1 (MQ=255)
                                                                                             aCTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCg                                  >  2:64832/1‑149 (MQ=255)
                                                                                                 tgCTACTGGATTTGGCTGCAAGGCGCGAGGGCGGCCGAGGCCTCCCGGCAATGTTCAGGTCAGATCTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCggg                                <  8:15296/147‑1 (MQ=255)
                                                                                                   cTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCg                             <  8:321563/148‑1 (MQ=255)
                                                                                                    tACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGtt                           >  3:220932/1‑149 (MQ=255)
                                                                                                         gCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGc                      >  2:8342/1‑149 (MQ=255)
                                                                                                          cTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCt                     <  6:323580/149‑1 (MQ=255)
                                                                                                          cTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCt                     >  3:145271/1‑149 (MQ=255)
                                                                                                           tATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                              >  1:177340/1‑91 (MQ=255)
                                                                                                           tATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                              <  2:177340/91‑1 (MQ=255)
                                                                                                           tATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTg                    <  4:149549/149‑1 (MQ=255)
                                                                                                           tATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTg                    >  3:76200/1‑149 (MQ=255)
                                                                                                               gCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCt                >  2:57320/1‑149 (MQ=255)
                                                                                                                cTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCGTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCtt               >  1:193313/1‑149 (MQ=255)
                                                                                                                 tGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCAGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGcatag              >  4:82282/1‑145 (MQ=255)
                                                                                                                  gCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                           <  6:28111/55‑1 (MQ=255)
                                                                                                                  gCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                           >  5:28111/1‑55 (MQ=255)
                                                                                                                      cgcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtc         <  8:301583/149‑1 (MQ=255)
                                                                                                                       gcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcg        <  4:260970/149‑1 (MQ=255)
                                                                                                                        cgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcga       >  5:322309/1‑149 (MQ=255)
                                                                                                                         tgcaGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGgt            <  3:136655/140‑1 (MQ=255)
                                                                                                                         tgcaGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGgt            >  4:136655/4‑143 (MQ=255)
                                                                                                                            aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGc                  >  7:140420/1‑134 (MQ=255)
                                                                                                                            aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGc                  <  8:140420/134‑1 (MQ=255)
                                                                                                                            aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgtc   >  5:339039/1‑149 (MQ=255)
                                                                                                                              gggCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgtcg  <  1:64832/149‑1 (MQ=255)
                                                                                                                              gggCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCGTGCTGTAGCCAAGGGCGGGCGGGCCGTTGACGCGGGCCTGGGtggcgatggcg  >  8:181133/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCG  >  NC_002947/593300‑593573

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