Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,820,202 A→T 6.6% L3663Q (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,2020AT 6.6% 180.2 / 3.3 61L3663Q (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 = 8.63e-01

CGCTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGAGC  >  NC_002947/4820057‑4820349
                                                                                                                                                 |                                                                                                                                                   
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                                                       <  2:322578/149‑1 (MQ=255)
                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCgg                                                       <  1:370495/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                           <  4:286898/149‑1 (MQ=255)
                                                                                                                       cTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTg                           >  7:297137/1‑149 (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)
                                                                                                                                                cagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGccc            <  7:78960/139‑1 (MQ=255)
                                                                                                                                                cagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGccc            >  8:78960/1‑139 (MQ=255)
                                                                                                                                                cagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGAGc  >  8:219875/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGCTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGAGC  >  NC_002947/4820057‑4820349

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