Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 593,892 A→G 8.0% L83P (CTG→CCG) ‡ oprG ← outer membrane protein OprG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947593,8920AG 8.0% 173.5 / 5.0 62L83P (CTG→CCG) ‡oprGouter membrane protein OprG
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base G (2/3);  total (28/34)
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.95e-01

GCGTACGGCTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCGGTGGTGATG  >  NC_002947/593745‑594038
                                                                                                                                                   |                                                                                                                                                   
gCGTACGGCTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCcag                                                                                                                                                    <  5:106378/149‑1 (MQ=255)
 cGTACGGCTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCcag‑c                                                                                                                                                  <  6:29889/149‑1 (MQ=255)
      ggCTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTc                                                                                                                                             >  1:330073/1‑149 (MQ=255)
        cTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCTGG‑CAGCtacc                                                                                                                                           <  4:330316/149‑4 (MQ=255)
        cTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCTGG‑CAGCtacc                                                                                                                                           <  4:335679/149‑4 (MQ=255)
        cTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCa                                                                                                                                           <  7:225247/149‑1 (MQ=255)
        cTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGGCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCTGG‑CAGCtacc                                                                                                                                           <  2:64258/149‑4 (MQ=255)
           gAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAgg                                                                                                                                         <  3:292186/148‑1 (MQ=255)
           gAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAg                                                                                                                                          <  5:123497/147‑1 (MQ=255)
           gAAGCGCGAATTTGCCTCCATTTCGTACTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACTGTGTTCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGc                                                                                                                                        <  5:202304/149‑1 (MQ=255)
               cgcgAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGAt                                                                                                                                    <  6:366597/149‑1 (MQ=255)
               cgcgAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGAt                                                                                                                                    <  3:234943/149‑1 (MQ=255)
                  gAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACTGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGAtg                                                                                                                                   <  8:146608/147‑1 (MQ=255)
                        ggCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCa                                                                                                                                           >  4:322182/1‑133 (MQ=255)
                        ggCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCa                                                                                                                                           <  3:322182/133‑1 (MQ=255)
                         gCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTg                                                                                                                          <  6:52013/149‑1 (MQ=255)
                          cTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGa                                                                                                                         <  3:288464/149‑1 (MQ=255)
                             cATTGGGTAGTACTGCAGCGACAGCGTTGGTGGAAGTTGCTTGATGTGGGCCAGCTTGCCGTCCAGCCCGCAGCACAGAGCTTTTACGACGATGGTGTGATTGATCGGGTTAGATTCCAG‑CAGCTCCAGGAAGATTTGGTGGGTGAGc                                                                                                                       <  6:152017/148‑1 (MQ=255)
                               ttGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAGCCAGCTCCAGGCCGATGTGGTCGGTGAGCAt                                                                                                                     <  1:34729/149‑1 (MQ=255)
                                 gggTAGTACTGCAGCGTCTGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGGTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGt                                                                                                                   <  8:54605/148‑1 (MQ=255)
                                  ggTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACC‑TTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAg                                                                                                                 <  6:300494/148‑1 (MQ=255)
                                   gtagtaCTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAgg                                                                                                                <  4:111858/149‑1 (MQ=255)
                                    tagtaCTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAg                                                                                                                 >  4:223231/1‑147 (MQ=255)
                                    tagtaCTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAg                                                                                                                 <  3:223231/147‑1 (MQ=255)
                                         cTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAg                                                                                                           <  2:172922/148‑1 (MQ=255)
                                          tGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGAtg                                                                                                                                   <  1:13135/123‑1 (MQ=255)
                                          tGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGAtg                                                                                                                                   >  2:13135/1‑123 (MQ=255)
                                          tGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAg                                                                                                           >  2:259222/1‑147 (MQ=255)
                                          tGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAg                                                                                                           <  1:259222/147‑1 (MQ=255)
                                            cAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGAt                                                                                                                                    <  6:249486/120‑1 (MQ=255)
                                            cAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGAt                                                                                                                                    >  5:249486/1‑120 (MQ=255)
                                              gCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAgg                                                                                                     >  5:283861/1‑149 (MQ=255)
                                                 aCAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGcc                                                                                                   >  5:326955/1‑148 (MQ=255)
                                                          gcggcAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTa                                                                                         >  3:165544/1‑149 (MQ=255)
                                                                            ccGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGGGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGcc                                                                       >  5:241112/1‑149 (MQ=255)
                                                                              gCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCg                                                                                      <  6:280103/132‑1 (MQ=255)
                                                                              gCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCg                                                                                      >  5:280103/1‑132 (MQ=255)
                                                                                 aGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGt                                                                  <  4:890/149‑1 (MQ=255)
                                                                                 aGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGt                                                                  <  8:28053/149‑1 (MQ=255)
                                                                                    ttGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGcc                                                               >  5:24850/1‑149 (MQ=255)
                                                                                            cAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCt                                                                                    <  3:196411/120‑1 (MQ=255)
                                                                                            cAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCt                                                                                    >  4:196411/1‑120 (MQ=255)
                                                                                                ccgccgCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCttgtt                                                   <  7:225950/149‑1 (MQ=255)
                                                                                                ccgccgCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCttgtt                                                   <  6:283861/149‑1 (MQ=255)
                                                                                                    cgccCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCg                                               <  6:214919/149‑1 (MQ=255)
                                                                                                          gACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAAc                                         <  6:326955/149‑1 (MQ=255)
                                                                                                             cTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACtt                                       >  3:198075/1‑148 (MQ=255)
                                                                                                              ttttACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACtt                                       >  6:160256/1‑147 (MQ=255)
                                                                                                               tttACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAg                                    >  7:51407/1‑149 (MQ=255)
                                                                                                                ttACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAg                                    <  2:330073/148‑1 (MQ=255)
                                                                                                                    gCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTcgc                               >  1:258278/1‑149 (MQ=255)
                                                                                                                       gACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTcgccgc                            <  4:165544/149‑1 (MQ=255)
                                                                                                                             gtgCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGt                                              >  4:316627/1‑125 (MQ=255)
                                                                                                                             gtgCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGACGTCGAACTTGAGGTCGCCGCTGctct                      >  3:260196/1‑149 (MQ=255)
                                                                                                                                cTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATAGCTGTCCAGGTTCGCCTTGGTGCCCGCCACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGt                   >  2:17902/1‑149 (MQ=255)
                                                                                                                                 ttGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGAATCGATGTCCAGGTTCGCCTTCGTGCCCGACACCTTGTTGACGTCGAACTTGAGGTCGCCGCTGCTCTCGtt                  >  7:45914/1‑149 (MQ=255)
                                                                                                                                       cctggaGCTGCTGG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGc            >  6:186297/6‑149 (MQ=255)
                                                                                                                                       cctggaGCTGCTGG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGc            >  4:53268/6‑149 (MQ=255)
                                                                                                                                           tAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCggtg        >  6:11307/1‑149 (MQ=255)
                                                                                                                                             gCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCggtggt      >  2:326540/1‑149 (MQ=255)
                                                                                                                                             gCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCggtggt      >  5:140442/1‑149 (MQ=255)
                                                                                                                                             gCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCggtggt      >  8:255785/1‑149 (MQ=255)
                                                                                                                                                 ccAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCGGTGGTGATg  <  6:140442/149‑1 (MQ=255)
                                                                                                                                                 ccAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCGGTGGTGATg  >  4:121245/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGTACGGCTGGAAGCGCGAATTTGGCTCCATTGGGTAGTACTGCAGCGACAGCGTTGGCGGCAGTTGCTTGATGTCCGCCAGCTTGCCGTCCAGCCCGCCGCCCAGACCTTTTACGCCGACGGTGTGCTTGAACGGGGTAGCTGCCAG‑CAGCTCCAGGCCGATGTGGTCGGTGAGCATGTAGGCAAAGGTCAGGCCCAACTGGGTATCGCTGTCCAGGGTCGCCTTGGTGCCCGACACCTTGTTGCCGTCGAACTTGAGGTCGCCGCTGCTCTCGTTCGGGGCGGTGGTGATG  >  NC_002947/593745‑594038

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