Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,820,199 A→T 6.6% L3664Q (CTG→CAG)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,820,1990AT 6.6% 176.6 / 3.1 61L3664Q (CTG→CAG) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base T (2/2);  total (28/33)
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.21e-01

CGCTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGA  >  NC_002947/4820057‑4820347
                                                                                                                                              |                                                                                                                                                    
cgcTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGACTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCAGCTCCCCTTCAAAACCTTGTACCAGCGCcagcag                                                                                                                                                  <  5:121972/147‑1 (MQ=255)
 gcTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCcagcagcag                                                                                                                                               >  1:59556/1‑149 (MQ=255)
 gcTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCcagcagcag                                                                                                                                               <  4:223347/149‑1 (MQ=255)
 gcTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCcagcagcag                                                                                                                                               <  5:325695/149‑1 (MQ=255)
   ttCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGgc                                                                                                                                             >  4:362597/1‑149 (MQ=255)
      gAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGgcgc                                                                                                                                           <  7:157545/148‑1 (MQ=255)
      gAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGgcgc                                                                                                                                           <  2:193524/148‑1 (MQ=255)
       aaCTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCTGCTGCTGGCGCtg                                                                                                                                         <  3:75349/149‑2 (MQ=255)
       aaCTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCtt                                                                                                                                         >  3:331567/1‑149 (MQ=255)
         cTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGa                                                                                                                                       <  5:62328/149‑1 (MQ=255)
            cgcgATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGaa                                                                                                                                      <  3:136863/147‑1 (MQ=255)
            cgcgATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAAtt                                                                                                                                    >  4:206889/1‑149 (MQ=255)
              cgATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCg                                                                                                                                  <  2:184404/149‑1 (MQ=255)
               gATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCgg                                                                                                                                 >  5:157734/1‑149 (MQ=255)
                     gCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAACCCTTGTACCAGCGCCTGCTGCTGGCGCtg                                                                                                                                         >  3:229984/1‑134 (MQ=255)
                     gCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCTGCTGCTGGCGCtg                                                                                                                                         <  4:229984/135‑2 (MQ=255)
                       tctGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAgg                                                                                                                         >  5:242594/1‑149 (MQ=255)
                         tGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGACcagcgctac                                                                                                                       >  4:356407/1‑144 (MQ=255)
                            tCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTc                                                                                                                    >  8:45843/1‑149 (MQ=255)
                              cAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGa                                                                                                                  >  5:133659/1‑149 (MQ=255)
                              cAGCGGATAAACCCGCTCGGTCCGGGTGCAGGCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCGCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGa                                                                                                                  <  8:75015/149‑1 (MQ=255)
                                   gATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTc                                                                                                             <  6:242594/149‑1 (MQ=255)
                                      aaaCCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCggg                                                                                                          >  1:322578/1‑149 (MQ=255)
                                        aCCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCggg                                                                                                          >  5:152712/1‑147 (MQ=255)
                                            gCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCg                                                                                                    <  3:206889/149‑1 (MQ=255)
                                               cGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAg                                                                                                 <  8:222688/149‑1 (MQ=255)
                                                      gTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTGACCTCCCCTTCAATACCTTGTAACAGCGCCAGCAGCAGGCGGTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTGGGAGTAGCGCt                                                                                          >  8:140204/1‑149 (MQ=255)
                                                       ttGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTTGGGTAGTCGAAGTAGCGCTg                                                                                         <  8:360711/149‑1 (MQ=255)
                                                             tcctccAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGaa                                                                                   <  4:371464/149‑1 (MQ=255)
                                                               ctccAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAAc                                                                                  <  8:113153/148‑1 (MQ=255)
                                                                  cagcagGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACa                                                                              <  7:45843/149‑1 (MQ=255)
                                                                     cagGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGt                                                                           >  7:299283/1‑149 (MQ=255)
                                                                        gAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGt                                                                           <  2:58544/146‑1 (MQ=255)
                                                                        gAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGt                                                                           >  1:58544/1‑146 (MQ=255)
                                                                             cTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGcc                                                                    >  2:370495/1‑148 (MQ=255)
                                                                             cTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTCGTCGAACGACAGGTGCAGGcc                                                                    >  5:90364/1‑148 (MQ=255)
                                                                              ttGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGcc                                                                    >  3:120514/1‑147 (MQ=255)
                                                                                gCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAg                                                                 <  2:59556/148‑1 (MQ=255)
                                                                                 cTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAgg                                                                <  3:231590/148‑1 (MQ=255)
                                                                                  tCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGt                                                               >  6:9640/1‑148 (MQ=255)
                                                                                      gggCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGtcatca                                                          <  8:299283/149‑1 (MQ=255)
                                                                                       ggCGCTCAGCAGCGGTAGCTCGCTCACCTGCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGtcatca                                                          <  3:87054/148‑1 (MQ=255)
                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCgg                                                     <  1:370495/149‑1 (MQ=255)
                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCgg                                                     <  2:322578/149‑1 (MQ=255)
                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCg                                                      <  2:88763/148‑1 (MQ=255)
                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCg                                                      >  1:88763/1‑148 (MQ=255)
                                                                                             cagcagCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGCCGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCggg                                                    >  6:59913/1‑148 (MQ=255)
                                                                                             cagcagCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCggg                                                    >  5:80309/1‑148 (MQ=255)
                                                                                               gcagcGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTaa                                                 <  5:268544/149‑1 (MQ=255)
                                                                                                     gTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACAcc                                           <  7:140204/149‑1 (MQ=255)
                                                                                                     gTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACAcc                                           >  6:169927/1‑149 (MQ=255)
                                                                                                          tggctcACCTCCCCTTCAACACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAgg                                                                      <  8:166422/115‑1 (MQ=255)
                                                                                                              tcACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGc                                  >  5:307800/1‑149 (MQ=255)
                                                                                                                   tCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGATTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCggg                               >  8:360394/1‑147 (MQ=255)
                                                                                                                       cTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTg                         >  7:297137/1‑149 (MQ=255)
                                                                                                                       cTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTg                         <  4:286898/149‑1 (MQ=255)
                                                                                                                       cTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTg                         <  6:80309/149‑1 (MQ=255)
                                                                                                                          cAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCTGTCAGCGGGAAGTTGgt                       <  5:9640/148‑1 (MQ=255)
                                                                                                                          cAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTgg                        <  6:152712/147‑1 (MQ=255)
                                                                                                                          cAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGgt                       <  7:360394/148‑1 (MQ=255)
                                                                                                                            aaaCCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGtcat                                                            <  6:330065/109‑1 (MQ=255)
                                                                                                                             aaCCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTgg                   <  6:90364/149‑1 (MQ=255)
                                                                                                                                cTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTa                >  3:196694/1‑149 (MQ=255)
                                                                                                                                 ttGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTAc               >  7:2053/1‑149 (MQ=255)
                                                                                                                                       aaGCGCCTGCTGCTGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGccc          >  7:63026/2‑148 (MQ=255)
                                                                                                                                              agcagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGa  <  4:196694/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CGCTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGA  >  NC_002947/4820057‑4820347

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