Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,066,004 G→T 5.2% N1020K (AAC→AAA mdtC ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,066,0040GT 5.2% 283.9 / 3.3 88N1020K (AAC→AAAmdtCmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/45);  new base T (2/3);  total (40/48)
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.10e-01

AGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACG  >  NC_002947/4065862‑4066152
                                                                                                                                               |                                                                                                                                                    
aGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGaaa                                                                                                                                                <  5:12660/149‑1 (MQ=255)
 gCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAAc                                                                                                                                               >  7:398201/1‑149 (MQ=255)
  cAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACg                                                                                                                                              >  7:276408/1‑149 (MQ=255)
  cAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCAGCGCCAACGGTTGAAACg                                                                                                                                              >  6:391999/1‑149 (MQ=255)
  cAAGGCCACGCACAACACGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACg                                                                                                                                              >  8:253548/1‑149 (MQ=255)
   aaGGCCACGCACAACCCGGCACCGCGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCCCGGTGAAGTGGGGTCTGGGCAAAACTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGGTGAAACgg                                                                                                                                             >  3:198807/1‑149 (MQ=255)
     ggCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACgg                                                                                                                                                      >  4:94965/1‑138 (MQ=255)
     ggCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACgg                                                                                                                                                      <  3:94965/138‑1 (MQ=255)
     ggCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTg                                                                                                                                           <  7:327035/149‑1 (MQ=255)
      gCCACGAACAACCCGGCACCGAGCAGGCGCGAGCCACGGCGACGGGTCAGTGCTTGCAGGGCAGGGTGAAGTGGGGTGTGGGCAAAAGTAATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTgg                                                                                                                                          <  4:198807/149‑1 (MQ=255)
       ccACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgc                                                                                                                                         <  3:79222/149‑1 (MQ=255)
         aCGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgcg                                                                                                                                        >  2:40718/1‑148 (MQ=255)
           gCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAg                                                                                                                                     <  5:84591/149‑1 (MQ=255)
            cacaACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCTGCGCCACCGTTTCAACCGTTGgcgcgg                                                                                                                                     <  6:236875/149‑3 (MQ=255)
            cacaACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAgg                                                                                                                                    >  6:41034/1‑149 (MQ=255)
             acaACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGTGCAGGc                                                                                                                                   <  2:292453/149‑1 (MQ=255)
             acaACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGGGAAGAGGGGTCTGGGCAAAAGTCATAGCGGGGGTTCCAGAGCGGCGTCGGTACGCACGCCGCCCCAACGGTTGAAACGGTGGCGGAGGc                                                                                                                                   >  8:308437/1‑149 (MQ=255)
               aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                 >  3:148276/1‑149 (MQ=255)
                 cccGGCACCGAGCAGGCGCGAGCCACGGCTACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                               <  4:310026/149‑1 (MQ=255)
                 cccGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                               <  5:391973/149‑1 (MQ=255)
                    ggCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcag                                                                                                                            <  6:82508/149‑1 (MQ=255)
                        cggAGCAGGCGCGAGCCACGGCCACGTGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCTGCGTCGGTACGAACGGCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtac                                                                                                                        <  7:59771/147‑1 (MQ=255)
                        ccGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtac                                                                                                                        <  3:198503/149‑1 (MQ=255)
                          gAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacag                                                                                                                      <  8:21551/149‑1 (MQ=255)
                          gAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacag                                                                                                                      >  3:249768/1‑149 (MQ=255)
                            gCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggt                                                                                                                    <  1:334495/149‑1 (MQ=255)
                             cAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggt                                                                                                                    >  8:275143/1‑148 (MQ=255)
                             cAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggt                                                                                                                    >  8:189361/1‑148 (MQ=255)
                              aGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTaa                                                                                                                  <  7:253548/149‑1 (MQ=255)
                              aGGCGCGAGCCACGGCCACGGGGCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggtat                                                                                                                  <  3:244309/149‑2 (MQ=255)
                                gcgcGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGt                                                                                                                                            <  4:169433/121‑1 (MQ=255)
                                gcgcGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGt                                                                                                                                            >  3:169433/1‑121 (MQ=255)
                                    gAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAaccacc                                                                                                            <  8:315235/149‑1 (MQ=255)
                                       ccacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGc                                                                                                         <  5:379257/149‑1 (MQ=255)
                                        cacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcg                                                                                                        >  3:52204/1‑149 (MQ=255)
                                          cggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtc                                                                                                      <  7:189361/149‑1 (MQ=255)
                                             ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacagg                                                                                                                     >  5:288937/1‑131 (MQ=255)
                                             ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacagg                                                                                                                     <  6:288937/131‑1 (MQ=255)
                                                   gTCAGTGCTTGCAGGGCACGGTTAAGTGGGGTCTGGGCAAAAGTCATTGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAgg                                                                                             <  5:427400/149‑1 (MQ=255)
                                                    tCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAggg                                                                                            <  1:40718/149‑1 (MQ=255)
                                                      aGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGCTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGt                                                                                           <  5:13448/148‑1 (MQ=255)
                                                       gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                               >  7:123232/1‑111 (MQ=255)
                                                       gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                               <  8:123232/111‑1 (MQ=255)
                                                          cTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                 <  5:417606/106‑1 (MQ=255)
                                                          cTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                 >  6:417606/1‑106 (MQ=255)
                                                           ttGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGAt                                                                                     >  6:325666/1‑149 (MQ=255)
                                                           ttGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGGCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGAt                                                                                     <  5:296246/149‑1 (MQ=255)
                                                             gCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCt                                                                                   >  4:103800/1‑149 (MQ=255)
                                                             gCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTGCCAGAGCGGCGTCGGTACGCAGGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCt                                                                                   <  4:7935/149‑1 (MQ=255)
                                                               aGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGCTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTgg                                                                                 >  1:11188/1‑149 (MQ=255)
                                                                      ggTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggt                                                                                                                    <  5:161416/107‑1 (MQ=255)
                                                                      ggTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggt                                                                                                                    >  6:161416/1‑107 (MQ=255)
                                                                          aaGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAg                                                                      >  3:96679/1‑149 (MQ=255)
                                                                           aggggggGCTTGGCATAAGTCATAGCGGGGTTTCCAGATGCGCGTGGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAg                                                                      <  8:178468/145‑1 (MQ=255)
                                                                           aGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAg                                                                      <  3:103800/148‑1 (MQ=255)
                                                                            gTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                         <  5:208671/96‑1 (MQ=255)
                                                                            gTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                         >  6:208671/1‑96 (MQ=255)
                                                                             tGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGaa                                                                            <  3:365486/140‑1 (MQ=255)
                                                                             tGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTCAGGATCTGGCTGaa                                                                            >  4:365486/1‑140 (MQ=255)
                                                                               gggTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCAc                                                                 <  8:350417/149‑1 (MQ=255)
                                                                                 gTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCAcc                                                                >  5:35707/1‑148 (MQ=255)
                                                                                 gTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCg                                                               <  5:391999/149‑1 (MQ=255)
                                                                                  tCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCga                                                              >  7:229261/1‑149 (MQ=255)
                                                                                  tCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCga                                                              >  2:143115/1‑149 (MQ=255)
                                                                                    tGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGGGGATCTGGCTGAACACCAGGCCACCgat                                                             >  4:327183/1‑148 (MQ=255)
                                                                                       gCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGatgatg                                                         >  2:83812/1‑149 (MQ=255)
                                                                                       gCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGatgatg                                                         <  4:249768/149‑1 (MQ=255)
                                                                                         aaaaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGt                                                       <  7:288878/149‑1 (MQ=255)
                                                                                         aaaaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGt                                                       <  4:52204/149‑1 (MQ=255)
                                                                                            aGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAg                                                    <  4:148276/149‑1 (MQ=255)
                                                                                                aTAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGGGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGcc                                                 <  3:327183/148‑1 (MQ=255)
                                                                                                      gggTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGc                                          >  5:282105/1‑149 (MQ=255)
                                                                                                        gTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGTCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTg                                        <  8:398201/149‑1 (MQ=255)
                                                                                                         tttCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGTTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCATGCGACCGATGATGGGCAGGCCCAGCGGCTGa                                       <  2:11188/149‑1 (MQ=255)
                                                                                                                agCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCAttt                                <  8:351525/149‑1 (MQ=255)
                                                                                                                        cGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCt                        >  7:145359/1‑149 (MQ=255)
                                                                                                                           tACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCa                     <  5:41034/149‑1 (MQ=255)
                                                                                                                            acgcacgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAg                    >  7:141201/1‑149 (MQ=255)
                                                                                                                                acgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTg                <  8:141201/149‑1 (MQ=255)
                                                                                                                                acgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTg                >  3:421920/1‑149 (MQ=255)
                                                                                                                                  gccgcGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGc              >  1:21983/1‑149 (MQ=255)
                                                                                                                                  gccgcGCCAACGGTTGAAACGGTGGCGCAGGCGGGCCAGGGACAGGTTAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGACCCCGAGGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGc              >  2:337212/1‑149 (MQ=255)
                                                                                                                                    ctgcgcCACCGTTTCAACCGTTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCt             >  8:386844/3‑149 (MQ=255)
                                                                                                                                    cgcgCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCt             >  4:74927/1‑148 (MQ=255)
                                                                                                                                        ggcttcGGTTGAAAGGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagc         <  7:308437/144‑1 (MQ=255)
                                                                                                                                        gcCACCGTTTCAACCGTTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGTTGATGGTCAGGcc                                                 >  1:356851/1‑109 (MQ=255)
                                                                                                                                        gcCACCGTTTCAACCGTTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGcc                                                 <  2:356851/109‑1 (MQ=255)
                                                                                                                                        gcCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcag          >  4:2675/1‑148 (MQ=255)
                                                                                                                                          cAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagcag       <  7:163763/149‑1 (MQ=255)
                                                                                                                                            aCGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGAAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCagcagc      >  3:14432/1‑148 (MQ=255)
                                                                                                                                               gTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACg  >  7:285655/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACG  >  NC_002947/4065862‑4066152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: