Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,319,011 T→G 8.5% F79V (TTC→GTC) ‡ PP_1152 → membrane fusion efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,319,0110TG 8.5% 182.1 / 7.9 59F79V (TTC→GTC) ‡PP_1152membrane fusion efflux protein
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base G (3/2);  total (30/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 = 7.76e-01

GCTCACCCCATGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCC  >  NC_002947/1318863‑1319157
                                                                                                                                                    |                                                                                                                                                  
gCTCACCCCATGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCGTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCt                                                                                                                                                    >  4:41886/1‑149 (MQ=255)
gCTCACCCCATGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGCCCTCGAGCGCt                                                                                                                                                    >  8:237629/1‑149 (MQ=255)
    aCCCCATGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTcc                                                                                                                                                 <  6:169395/148‑1 (MQ=255)
                 gTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAgg                                                                                                                                              <  5:148330/138‑1 (MQ=255)
                 gTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAgg                                                                                                                                              >  6:148330/1‑138 (MQ=255)
                        ccGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCa                                                                                                                            <  2:94165/149‑1 (MQ=255)
                                  gcCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGAc                                                                                                                  >  7:44496/1‑149 (MQ=255)
                                            tggtgATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACAcgc                                                                                                        >  8:50647/1‑149 (MQ=255)
                                                         tGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTgc                                                                                           >  7:161976/1‑149 (MQ=255)
                                                             gtgtCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCga                                                                                                               >  2:138383/1‑125 (MQ=255)
                                                             gtgtCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCga                                                                                                               <  1:138383/125‑1 (MQ=255)
                                                             gtgtCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACATGCGCCCAGCAACTgc                                                                                           >  8:33643/1‑145 (MQ=255)
                                                             gtgtCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACATGCGCCCAGCAACTgc                                                                                           <  7:33643/145‑1 (MQ=255)
                                                              tgtCCGGCTGGGTGCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCtt                                                                                       <  2:94135/148‑1 (MQ=255)
                                                                tCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCa                                                                                                                                 >  7:83328/1‑104 (MQ=255)
                                                                tCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCa                                                                                                                                 <  8:83328/104‑1 (MQ=255)
                                                                             gcgAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAag                                                                         >  5:49410/1‑147 (MQ=255)
                                                                                   tCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGcc                                                                    <  2:294025/146‑1 (MQ=255)
                                                                                   tCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGcc                                                                    >  1:294025/1‑146 (MQ=255)
                                                                                     aaGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCg                                                               >  8:142166/1‑149 (MQ=255)
                                                                                     aaGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCg                                                               >  5:265686/1‑149 (MQ=255)
                                                                                          ccATGACAAGCAGGAGTTCAAGGCTGGCGTCTTGTTTATGAGTATTGACCGCGAGCGCTTACAGGCAGCGCTCGATCAGGCCAGCGCGTTGACCGATACAGGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                                           <  3:37164/148‑1 (MQ=255)
                                                                                          ccATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTg                                                          >  8:108548/1‑149 (MQ=255)
                                                                                               aCAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGcttc                                                      <  6:265686/148‑4 (MQ=255)
                                                                                                cAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTgg                                                    <  8:161976/149‑1 (MQ=255)
                                                                                                    cagcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTg                                                          <  2:12415/139‑1 (MQ=255)
                                                                                                    cagcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTg                                                          >  1:12415/1‑139 (MQ=255)
                                                                                                       cagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCga                                                                                                               >  8:105242/1‑83 (MQ=255)
                                                                                                          gTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                          >  4:128391/1‑149 (MQ=255)
                                                                                                          gTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                          >  8:54833/1‑149 (MQ=255)
                                                                                                           tCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCACGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                                          >  5:223305/1‑148 (MQ=255)
                                                                                                             aaGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGa                                        <  6:49410/148‑1 (MQ=255)
                                                                                                                 cTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGAt                                       <  7:167730/145‑1 (MQ=255)
                                                                                                                 cTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGAt                                       >  8:167730/1‑145 (MQ=255)
                                                                                                                     cGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGcc                                >  1:126551/1‑148 (MQ=255)
                                                                                                                         ttgttgATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGtt                           >  5:304369/1‑149 (MQ=255)
                                                                                                                            ttgATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGc                         <  7:237629/148‑1 (MQ=255)
                                                                                                                            ttgATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCg                        <  7:142166/149‑1 (MQ=255)
                                                                                                                              gatgaGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGt                       <  3:153333/148‑1 (MQ=255)
                                                                                                                                tgaGTATTGACCGCGAGCGGTTCCAGGGCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGaa                    <  7:102227/149‑1 (MQ=255)
                                                                                                                                tgaGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGaa                    <  3:41886/149‑1 (MQ=255)
                                                                                                                                tgaGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGa                     >  6:12284/1‑148 (MQ=255)
                                                                                                                                 gaGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTgga                    >  5:303610/1‑146 (MQ=255)
                                                                                                                                 gaGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGaa                    >  8:101269/1‑148 (MQ=255)
                                                                                                                                  aGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGaaaa                  >  2:192586/1‑149 (MQ=255)
                                                                                                                                        gACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGcc              >  8:167714/1‑147 (MQ=255)
                                                                                                                                         aCCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCa            <  6:304369/148‑1 (MQ=255)
                                                                                                                                          gatcgAGCGCGGCCTGGAAGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTggg                                                   <  5:134098/105‑1 (MQ=255)
                                                                                                                                          gatcgAGCGCGGCCTGGAAGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTggg                                                   >  6:134098/4‑108 (MQ=255)
                                                                                                                                          gatcgAGCGCGGCCTGGAAGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGAt                                       <  4:100032/117‑1 (MQ=255)
                                                                                                                                          gatcgAGCGCGGCCTGGAAGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGAt                                       >  3:100032/4‑120 (MQ=255)
                                                                                                                                          ccGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCagag                                            <  3:116645/115‑1 (MQ=255)
                                                                                                                                          ccGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCagag                                            <  3:112007/115‑1 (MQ=255)
                                                                                                                                          ccGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCagag                                            >  4:116645/1‑115 (MQ=255)
                                                                                                                                          ccGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCagag                                            >  4:112007/1‑115 (MQ=255)
                                                                                                                                           cgcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGAt         <  7:54833/149‑1 (MQ=255)
                                                                                                                                            tcgAGCGCGGCCTGGAAGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGAt         >  1:275379/2‑148 (MQ=255)
                                                                                                                                                gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATg    <  7:50647/149‑1 (MQ=255)
                                                                                                                                                 cgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATgc   <  7:101269/149‑1 (MQ=255)
                                                                                                                                                  gcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATgcc  <  5:12284/149‑1 (MQ=255)
                                                                                                                                                  gcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATgcc  <  7:108548/149‑1 (MQ=255)
                                                                                                                                                    ttCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAg                             >  8:124557/1‑120 (MQ=255)
                                                                                                                                                    ttCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAg                             <  7:124557/120‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTCACCCCATGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACGCCCAGATCAATGCC  >  NC_002947/1318863‑1319157

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: