Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 47,141 G→C 7.6% A745G (GCG→GGG) ‡ czcA‑I ← cation efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294747,1410GC 7.6% 195.9 / 5.1 64A745G (GCG→GGG) ‡czcA‑Ication efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (25/32);  new base C (2/3);  total (28/36)
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.71e-01

CCGGGATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGA  >  NC_002947/46995‑47285
                                                                                                                                                   |                                                                                                                                                
ccGGGATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCcgc                                                                                                                                                 <  6:211606/148‑1 (MQ=255)
    gATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGt                                                                                                                                            >  5:102043/1‑149 (MQ=255)
        aGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTTATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCGCGGATGTTCa                                                                                                                                         <  4:18794/148‑1 (MQ=255)
        aGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCGCGGATGTTCAg                                                                                                                                        <  8:340747/149‑1 (MQ=255)
         gCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGAt                                                                                                                                              <  4:211963/142‑1 (MQ=255)
         gCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGAt                                                                                                                                              >  3:211963/1‑142 (MQ=255)
          cAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGc                                                                                                                                      >  6:127477/1‑149 (MQ=255)
              cTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTa                                                                                                                                  <  4:31055/149‑1 (MQ=255)
               tGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAg                                                                                                                                 >  1:16496/1‑149 (MQ=255)
                 gACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCTCCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCg                                                                                                                               >  1:302079/1‑149 (MQ=255)
                      aCCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGtgctg                                                                                                                          >  8:249049/1‑149 (MQ=255)
                          gCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCttt                                                                                                                      >  4:41045/1‑149 (MQ=255)
                              cGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCGGACGGCCACCCACGGCGATAGCGATCGAGTGCTGAACATCCGCGATGTTCAGGCGGTAGCGTGCCGCTTTTTc                                                                                                                   <  5:275457/148‑1 (MQ=255)
                                tCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTcgcg                                                                                                                <  5:127477/149‑1 (MQ=255)
                                      cGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAAt                                                                                                          <  7:335265/149‑1 (MQ=255)
                                       gAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATg                                                                                                         >  7:316068/1‑149 (MQ=255)
                                           gggctctGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACACACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGa                                                                                                     <  1:306731/146‑1 (MQ=255)
                                            gTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGAt                                                                                                    >  8:1359/1‑149 (MQ=255)
                                            gTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGAt                                                                                                    <  7:270439/149‑1 (MQ=255)
                                             tctcTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATg                                                                                                   <  8:316068/149‑1 (MQ=255)
                                                 tGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGt                                                                                                        <  3:359633/140‑1 (MQ=255)
                                                 tGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGt                                                                                                        >  4:359633/1‑140 (MQ=255)
                                                 tGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATgg                                                                                                  <  1:171746/146‑1 (MQ=255)
                                                   gTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAg                                                                                              >  5:192850/1‑148 (MQ=255)
                                                   gTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAg                                                                                              >  4:197065/1‑148 (MQ=255)
                                                   gTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAg                                                                                              >  3:319034/1‑148 (MQ=255)
                                                         gTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCgg                                                                                        >  1:318103/1‑148 (MQ=255)
                                                          tactacCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCa                                                                                      <  7:310989/149‑1 (MQ=255)
                                                             tacCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGc                                                                                   >  8:260692/1‑149 (MQ=255)
                                                               ccATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTTCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGc                                                                                   <  8:178134/147‑1 (MQ=255)
                                                               ccATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGc                                                                                   <  2:318103/147‑1 (MQ=255)
                                                               ccATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCg                                                                                 <  7:1359/149‑1 (MQ=255)
                                                                cATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCgg                                                                                <  2:331967/149‑1 (MQ=255)
                                                                cATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCgg                                                                                <  6:192850/149‑1 (MQ=255)
                                                                    tCGCACCGACGGGACACCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGAt                                                                              <  1:311344/148‑1 (MQ=255)
                                                                      gAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGCACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCt                                                                          >  5:264499/1‑149 (MQ=255)
                                                                         ccGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCt                                                                       <  7:249049/149‑1 (MQ=255)
                                                                          cGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAg                                                                                              >  3:336663/1‑125 (MQ=255)
                                                                          cGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAg                                                                                              <  4:336663/125‑1 (MQ=255)
                                                                                    cccTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCAc                                                             <  3:41045/149‑1 (MQ=255)
                                                                                     ccTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACt                                                            >  6:118168/1‑149 (MQ=255)
                                                                                     ccTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACt                                                            >  7:279847/1‑149 (MQ=255)
                                                                                         aTACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGC‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGa                                                       >  2:160463/1‑149 (MQ=255)
                                                                                           aCAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGAt                                                      <  4:220920/149‑1 (MQ=255)
                                                                                             aGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGAt                                                      >  2:51161/1‑147 (MQ=255)
                                                                                             aGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGAt                                                      <  1:51161/147‑1 (MQ=255)
                                                                                                gTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGcc                                                 <  8:279847/149‑1 (MQ=255)
                                                                                                 tGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGAt                                                      >  8:183400/1‑143 (MQ=255)
                                                                                                 tGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGAt                                                      <  7:183400/143‑1 (MQ=255)
                                                                                                     ggCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGGCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGa                                            <  6:269233/149‑1 (MQ=255)
                                                                                                      gCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCACGGGAc                                           >  7:299413/1‑149 (MQ=255)
                                                                                                            cGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCttt                                     <  5:118168/149‑1 (MQ=255)
                                                                                                             ggCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCttt                                     >  2:39579/1‑148 (MQ=255)
                                                                                                              gCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGC‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAg                                  <  1:160463/149‑1 (MQ=255)
                                                                                                               ccacccacGGCGATAGGGATCGAGTTCTGAACATCCTGGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCTGTCAATGTTGATGGTCAGCACCGGCAGGCAGGATGTCTTCTCAACTTTCACTTCCGATGAGCCCGGGACCGGTTTGa                                   <  6:264499/148‑1 (MQ=255)
                                                                                                                cacccacGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAgc                                 <  3:197065/149‑1 (MQ=255)
                                                                                                                   ccacGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCt                              <  7:260692/149‑1 (MQ=255)
                                                                                                                            tAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATc                     <  4:319034/149‑1 (MQ=255)
                                                                                                                                   tacggccTGAACATCGCGGATGTTCAGACCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTgg              >  7:181662/7‑149 (MQ=255)
                                                                                                                                         cTGAACATCGCGGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTc                                                          >  7:73804/1‑99 (MQ=255)
                                                                                                                                         cTGAACATCGCGGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTc                                                          <  8:73804/99‑1 (MQ=255)
                                                                                                                                         cTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGgtg         <  2:16496/148‑1 (MQ=255)
                                                                                                                                         cTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTAAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGgtgt        >  6:273929/1‑149 (MQ=255)
                                                                                                                                          tGAACATCCTCGAGGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACGTCCGATGAGCCAGGGACCGCTTTGAGCGCTGCCGCGCGCTTGTTGGCGgtggt       >  8:16279/1‑147 (MQ=255)
                                                                                                                                               aTCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGa  <  7:17638/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGGGATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGAC‑GGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCTTTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGA  >  NC_002947/46995‑47285

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