Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,743,758 G→A 8.8% V41V (GTG→GTA) ‡ PP_2398 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,743,7580GA 8.8% 154.3 / 4.0 57V41V (GTG→GTA) ‡PP_2398hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base A (3/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GCGTTGTGCCTGTTGGAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTTCTCGACCGATACCAG  >  NC_002947/2743611‑2743902
                                                                                                                                                   |                                                                                                                                                
gCGTTGTGCCTGTTGGAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTgg                                                                                                                                                 >  1:210210/1‑149 (MQ=255)
 cGTTGTGCCTGTTGGAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTggg                                                                                                                                                <  4:226624/149‑1 (MQ=255)
         cTGTTGGAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGcc                                                                                                                                         >  5:192096/1‑148 (MQ=255)
           gTTGGAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTAt                                                                                                                                      >  6:40266/1‑149 (MQ=255)
            ttGGAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTAt                                                                                                                                      <  5:158581/148‑1 (MQ=255)
              ggAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATgg                                                                                                                                    >  1:181012/1‑148 (MQ=255)
               gagaACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGc                                                                                                                                   >  8:363407/1‑148 (MQ=255)
                ctaaCTTTTTTGTTTCCTGCGATGTTGAACGTTGGCCTAGGCGTGCGTCGCGGCCTGTTGCCGGAACTGCGGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGGGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCaa                                                                                                                                 <  6:39934/147‑1 (MQ=255)
                 gaACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGTGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGTCAAg                                                                                                                                <  5:288310/149‑1 (MQ=255)
                      tttATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCtt                                                                                                                           >  4:284488/1‑149 (MQ=255)
                       ttATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCtt                                                                                                                           <  7:183286/148‑1 (MQ=255)
                       ttATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCATAGGCGCCGCCCACGGCaatccag                                                                                                                            >  6:416537/1‑142 (MQ=255)
                       ttATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCATAGGCGCCGCCCACGGCaatccag                                                                                                                            <  5:416537/147‑6 (MQ=255)
                        tATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTg                                                                                                                          <  5:40266/148‑1 (MQ=255)
                            tttCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGcc                                                                                                                       >  4:368791/1‑147 (MQ=255)
                              tCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCTGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCtt                                                                                                                           >  1:212686/1‑141 (MQ=255)
                              tCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCTGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCtt                                                                                                                           <  2:212686/141‑1 (MQ=255)
                                     aTGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGa                                                                                                            <  3:368791/149‑1 (MQ=255)
                                          gAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAg                                                                                                           >  5:273208/1‑145 (MQ=255)
                                            aCGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCttt                                                                                                     <  3:189391/149‑1 (MQ=255)
                                             cGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTg                                                                                                    <  3:284488/149‑1 (MQ=255)
                                             cGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTg                                                                                                    <  6:192096/149‑1 (MQ=255)
                                               ttGGCTTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCt                                                                                                  <  8:154969/149‑1 (MQ=255)
                                                  gCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTc                                                                                               <  5:26387/149‑1 (MQ=255)
                                                              gTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATggg                                                                                   <  3:356682/149‑1 (MQ=255)
                                                                      cTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCttt                                                                           <  2:76395/149‑1 (MQ=255)
                                                                         ttGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTa                                                                          <  4:295415/147‑1 (MQ=255)
                                                                             cGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTg                                                                    >  4:329650/1‑149 (MQ=255)
                                                                              ggAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCCGAAAACTGGATTGACGAGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTgg                                                                   >  5:288922/1‑149 (MQ=255)
                                                                               gAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTggg                                                                  >  3:103969/1‑149 (MQ=255)
                                                                                      tgGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTCCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGt                                                           >  5:358856/1‑149 (MQ=255)
                                                                                          cATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGccc                                                       >  4:326868/1‑149 (MQ=255)
                                                                                           aTGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCCTTGCTGTCACCTGTCATgg                                                                                    <  6:22240/119‑1 (MQ=255)
                                                                                           aTGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCCTTGCTGTCACCTGTCATgg                                                                                    >  5:22240/1‑119 (MQ=255)
                                                                                              gAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGc                                                    <  7:363407/148‑1 (MQ=255)
                                                                                                 gCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGa                                                <  1:388864/149‑1 (MQ=255)
                                                                                                          gTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCATTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCg                                                     <  8:100381/135‑1 (MQ=255)
                                                                                                          gTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCg                                                     >  7:100381/1‑135 (MQ=255)
                                                                                                            ggACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTcc                                     >  4:11796/1‑149 (MQ=255)
                                                                                                              aCTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCtt                                   <  3:326868/149‑1 (MQ=255)
                                                                                                                tttCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTCCTGTCCCCTGTAATGGGTTGTCTTTATCGCTGGGTGGCAGAGCCCCGCTg                                                  >  6:409006/1‑132 (MQ=255)
                                                                                                                  tCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGa                                 <  3:11796/147‑1 (MQ=255)
                                                                                                                     ggAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGa                            >  4:369022/1‑149 (MQ=255)
                                                                                                                      gAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCtt                                   >  5:123597/1‑141 (MQ=255)
                                                                                                                      gAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCtt                                   <  6:123597/141‑1 (MQ=255)
                                                                                                                      gAGTGTGCGCCGGAAAACTGGATTGCCGTAGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGAc                           >  8:61002/1‑149 (MQ=255)
                                                                                                                           tgCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCgg                       <  4:103969/148‑1 (MQ=255)
                                                                                                                             cgcCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAg                    >  2:279924/1‑149 (MQ=255)
                                                                                                                              gcCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCt                                                   <  7:362197/117‑1 (MQ=255)
                                                                                                                              gcCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATAGCTGGGTGGCAGTGCCCCGCt                                                   >  8:362197/1‑117 (MQ=255)
                                                                                                                                   aaaaCTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTtctc              <  6:358856/149‑1 (MQ=255)
                                                                                                                                      aCTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTTCTCGAc           <  3:329650/149‑1 (MQ=255)
                                                                                                                                          gATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGAccc                         <  7:22033/131‑1 (MQ=255)
                                                                                                                                          gATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGAccc                         >  8:22033/1‑131 (MQ=255)
                                                                                                                                          gATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTTCTCGACCGAt       <  2:210210/149‑1 (MQ=255)
                                                                                                                                          ccttGCCATAGGCGCCGCCCACGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTtctc              <  1:348533/140‑1 (MQ=255)
                                                                                                                                          ccttGCCATAGGCGCCGCCCACGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTtctc              >  2:348533/3‑142 (MQ=255)
                                                                                                                                               tcGTTGGCGGCGACTATGGCAAGGTCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTTCTCGACCGATACCAg  <  6:288922/148‑1 (MQ=255)
                                                                                                                                               ccGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTTCTCGACCGATACCAg  >  1:408089/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGTTGTGCCTGTTGGAGAACTTTTATGTTTCCTGCGATGTTGAACGTTGGCCTAGGCTTGCGTCGCGGCCTGTTGCCGGAACTGCTGGCCATGGAGGCGGGTGCAGTGGACTTTCTGGAGTGTGCGCCGGAAAACTGGATTGCCGTGGGCGGCGCCTATGGCAAGGGCTTGGCCCAGTTGGCGGAGCGCTTTGCTGTCACCTGTCATGGGTTGTCTTTATCGCTGGGTGGCAGTGCCCCGCTGGACCGACATTTCCTTGAGCAGACCCGGCAGTTTCTCGACCGATACCAG  >  NC_002947/2743611‑2743902

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: