Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 687,721 T→A 7.2% L196M (TTG→ATG)  ydgK → Bcr/CflA family multidrug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947687,7210TA 7.2% 193.0 / 3.3 64L196M (TTG→ATG) ydgKBcr/CflA family multidrug resistance transporter
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base A (2/3);  total (31/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 = 1.00e+00

GGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGAGCACTACGG  >  NC_002947/687580‑687869
                                                                                                                                             |                                                                                                                                                    
ggggtgCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCa                                                                                                                                               <  2:329815/149‑1 (MQ=255)
    tgCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCAAGGCGCCAgaca                                                                                                                                           <  2:335582/149‑4 (MQ=255)
      cTGGTGAACCTGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCAAGGCGCCAgacag                                                                                                                                          <  2:211009/148‑5 (MQ=255)
              ccTGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGTCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGgctgc                                                                                                                                 >  6:9872/1‑149 (MQ=255)
               cTGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGgctgct                                                                                                                                <  7:272558/149‑1 (MQ=255)
                tGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGctgctg                                                                                                                               >  6:285369/1‑149 (MQ=255)
                 ggCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTgctgg                                                                                                                              >  2:324685/1‑149 (MQ=255)
                     ggCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGa                                                                                                                          <  1:199567/149‑1 (MQ=255)
                      gCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTgc                                                                                                                                                   <  8:77658/123‑1 (MQ=255)
                      gCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTgc                                                                                                                                                   >  7:77658/1‑123 (MQ=255)
                      gCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGAc                                                                                                                         <  4:205745/149‑1 (MQ=255)
                          tcattCGATCTTCCTTTCGCTGAGCCTGTTCAGTGCCGGCTGCCTTCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCggg                                                                                                                     <  3:218024/145‑1 (MQ=255)
                          gCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCg                                                                                                                       <  6:109650/147‑1 (MQ=255)
                           cAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGt                                                                                                                    <  6:290852/149‑1 (MQ=255)
                              tCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGttt                                                                                                                  <  2:258240/148‑1 (MQ=255)
                              tCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGttt                                                                                                                  >  1:258240/1‑148 (MQ=255)
                                      ccTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCg                                                                                                            <  2:180087/146‑1 (MQ=255)
                                      ccTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGc                                                                                                          >  1:164600/1‑148 (MQ=255)
                                            gcTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGttt                                                                                                                  >  8:65846/1‑134 (MQ=255)
                                            gcTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGGCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGttt                                                                                                                  <  7:65846/134‑1 (MQ=255)
                                              tGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCt                                                                                                  >  2:272690/1‑148 (MQ=255)
                                                  ccTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCAAAAACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACcg                                                                                             >  4:61839/1‑149 (MQ=255)
                                                    tGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACcggc                                                                                           >  4:193482/1‑149 (MQ=255)
                                                      ttCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTgc                                                                                                                                                   <  1:74697/91‑1 (MQ=255)
                                                      ttCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTgc                                                                                                                                                   >  2:74697/1‑91 (MQ=255)
                                                         aGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGt                                                                                        <  5:340396/147‑1 (MQ=255)
                                                         aGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGt                                                                                        >  6:340396/1‑147 (MQ=255)
                                                          gTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGTCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATc                                                                                     <  5:9872/149‑1 (MQ=255)
                                                          gTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATc                                                                                     >  7:23893/1‑149 (MQ=255)
                                                             ccGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGc                                                                                   >  2:56992/1‑148 (MQ=255)
                                                             ccGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATAAACCCCGCCAACCCCAGTAATGCGCCATGCGCCAGTACCTTCGGCTGCTGGCATCCCGGGGTTTTCAACGCCAAGCCCTTCCCGGCGGTAtacc                                                                                   >  1:263529/1‑145 (MQ=255)
                                                                gCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCAtt                                                                               >  5:222168/1‑149 (MQ=255)
                                                                             ggTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCAt                                                                                                       <  6:88036/112‑1 (MQ=255)
                                                                             ggTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCAt                                                                                                       >  5:88036/1‑112 (MQ=255)
                                                                               tAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGcc                                                                                                    >  5:41476/1‑113 (MQ=255)
                                                                               tAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGcc                                                                                                    <  6:41476/113‑1 (MQ=255)
                                                                                 aGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCTGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCt                                                               <  4:126746/148‑1 (MQ=255)
                                                                                          cTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCgg                                                     <  1:272690/149‑1 (MQ=255)
                                                                                            gCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGtt                                                   >  8:184933/1‑149 (MQ=255)
                                                                                             ccGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGtt                                                   >  3:119063/1‑148 (MQ=255)
                                                                                             ccGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGtt                                                   <  1:324685/148‑1 (MQ=255)
                                                                                             ccGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGtt                                                   >  5:91517/1‑148 (MQ=255)
                                                                                               ggAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCg                                                                           <  2:155638/122‑1 (MQ=255)
                                                                                               ggAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCg                                                                           >  1:155638/1‑122 (MQ=255)
                                                                                                     ccTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCg                                          >  3:74189/1‑149 (MQ=255)
                                                                                                        gCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCgtgt                                       >  1:106034/1‑149 (MQ=255)
                                                                                                          ccGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTtc                                     <  5:285369/149‑1 (MQ=255)
                                                                                                                atatACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCaaa                               <  2:263529/149‑1 (MQ=255)
                                                                                                                 tataCCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAAc                              >  2:62766/1‑149 (MQ=255)
                                                                                                                  ataCCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGTGGTATCGCCATTGCCGGCATGTTCGCCtat                                                             <  8:290927/117‑1 (MQ=255)
                                                                                                                  ataCCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGTGGTATCGCCATTGCCGGCATGTTCGCCtat                                                             >  7:290927/1‑117 (MQ=255)
                                                                                                                  ataCCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACt                             >  7:194197/1‑149 (MQ=255)
                                                                                                                     ccccGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTa                          <  1:62766/149‑1 (MQ=255)
                                                                                                                     ccccGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGCGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTAGTGTTAATCAAACTGTa                          >  7:107474/1‑149 (MQ=255)
                                                                                                                            aaCCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTGTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTg                   <  4:23299/149‑1 (MQ=255)
                                                                                                                               cGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGccc                >  6:28389/1‑149 (MQ=255)
                                                                                                                                  gtacTGGCGCAAGGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATg                        <  7:232428/135‑1 (MQ=255)
                                                                                                                                  gtacTGGCGCAAGGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATg                        >  8:232428/4‑138 (MQ=255)
                                                                                                                                   gTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTg            <  3:193482/149‑1 (MQ=255)
                                                                                                                                    tCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGa           <  2:106034/149‑1 (MQ=255)
                                                                                                                                    tCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGa           <  2:164600/149‑1 (MQ=255)
                                                                                                                                     cTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCGTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGa           <  3:61839/148‑1 (MQ=255)
                                                                                                                                          gcCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGAGCACTa     >  7:142943/1‑149 (MQ=255)
                                                                                                                                            cTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGAGCACTACg   <  7:184933/149‑1 (MQ=255)
                                                                                                                                             ttGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGAGCACTACgg  >  3:221040/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GGGGTGCTGGTGAACCTGGCCGGCTGGCAGTCGATCTTCCTTGCGCTGAGCCTGTTCAGTGCCGGCTGCCTGCTTGCGGTAAGCCTGGGCCTGCCGGAAAGCCTGCCCGAGCATATACCCCGCCAACCGCTGTCTGGCGCCTTGCGCCAGTACCTGCGGCTGCTGGCTGACCGGGTTTTTCTCGGCCATGCCCTGACCGGCGGTATCGCCATTGCCGGCATGTTCGCCTATATCGCCGGTTCACCTTTCGTGTTCATCAAACTGTATGGCGTGCCCGCTGAGCACTACGG  >  NC_002947/687580‑687869

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: