Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 926,157 G→A 8.5% R157W (CGG→TGG)  PP_0805 ← outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947926,1570GA 8.5% 164.5 / 4.4 59R157W (CGG→TGG) PP_0805outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (2/3);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

GCGGGCCATGCCGGCAATGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAAT  >  NC_002947/926011‑926301
                                                                                                                                                  |                                                                                                                                                
gCGGGCCATGCCGGCAATGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCGAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCggg                                                                                                                                                <  2:111506/149‑1 (MQ=255)
     ccATGCCGGCAATGTCGCCCACGCCGCAAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGggcggcg                                                                                                                                           >  1:15095/1‑149 (MQ=255)
         gCCGGCAATGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCAggctgccgccg                                                                                                                                       >  2:225228/1‑147 (MQ=255)
          ccGGCAATGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGc                                                                                                                                      >  6:306711/1‑149 (MQ=255)
          ccGGCAATGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGc                                                                                                                                      <  5:306711/149‑1 (MQ=255)
             gCAATGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTg                                                                                                                                   >  3:273087/1‑149 (MQ=255)
                 tGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCTTAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCAGGCTGCCGCCGCCCGGGcc                                                                                                                               <  4:234880/149‑1 (MQ=255)
                 tGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGcc                                                                                                                               >  3:214922/1‑149 (MQ=255)
                  gTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCcg                                                                                                                              >  6:324394/1‑149 (MQ=255)
                     gCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCAGGCTGCCGCCGCCCGGGcccag                                                                                                                            <  1:225228/148‑4 (MQ=255)
                     gCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCAGGCTGCCGCCGCCCGGGcccag                                                                                                                            <  4:167800/148‑4 (MQ=255)
                             cgcgAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGTCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGgcg                                                                                                                   >  3:112406/1‑149 (MQ=255)
                             cgcgAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGgcg                                                                                                                   <  1:1786/149‑1 (MQ=255)
                                       gTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGAc                                                                                                           <  1:270482/147‑1 (MQ=255)
                                             ggCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCg                                                                                                    >  6:153520/1‑148 (MQ=255)
                                              gCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGaa                                                                                                  <  7:155490/149‑1 (MQ=255)
                                                aaTCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGt                                                                                                <  4:273087/149‑1 (MQ=255)
                                                  tCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCa                                                                                              >  3:237268/1‑149 (MQ=255)
                                                    gTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCAt                                                                                             <  2:15095/148‑1 (MQ=255)
                                                    gTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATa                                                                                            <  5:219475/149‑1 (MQ=255)
                                                       atcaGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGGCGACCACCAGCAGGATATTGGCCTGGGCCAGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATaga                                                                                          <  7:287135/148‑1 (MQ=255)
                                                         caGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCa                                                                                       >  7:193814/1‑149 (MQ=255)
                                                            ttGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCAttt                                                                                    >  5:282566/1‑149 (MQ=255)
                                                               cTCGTAGCGGCGGGTCTCGATCCAGTCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGc                                                                                  <  4:112406/148‑1 (MQ=255)
                                                               cTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGc                                                                                  <  4:214922/148‑1 (MQ=255)
                                                               cTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGc                                                                                  <  1:61041/148‑1 (MQ=255)
                                                                tCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCtt                                                                                >  1:72652/1‑149 (MQ=255)
                                                                tCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCACGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCtt                                                                                >  2:60448/1‑149 (MQ=255)
                                                                  gTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGa                                                                              <  2:191816/149‑1 (MQ=255)
                                                                      cggcggGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGa                                                                          <  1:66001/149‑1 (MQ=255)
                                                                      cggcggGTCTCGATCCAGGCATAGGCTGTGTCACGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGa                                                                          <  1:60448/149‑1 (MQ=255)
                                                                          gggTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAgg                                                                      >  4:260756/1‑149 (MQ=255)
                                                                               tcGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCAACAGGTTCAg                                                                 >  3:295176/1‑149 (MQ=255)
                                                                                cGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAgg                                                                >  3:39841/1‑149 (MQ=255)
                                                                                 gATCCAGGCATAGGCTGTGTCGCGGGCCAGGTGGTCGACGACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAgg                                                                <  5:153520/148‑1 (MQ=255)
                                                                                 gATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAggg                                                               >  3:47762/1‑149 (MQ=255)
                                                                                 gATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAgg                                                                <  4:237268/148‑1 (MQ=255)
                                                                                  aTCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCa                                                                                              >  6:2801/1‑117 (MQ=255)
                                                                                  aTCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCa                                                                                              <  5:2801/117‑1 (MQ=255)
                                                                                     cAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCt                                                            <  2:13396/148‑1 (MQ=255)
                                                                                      aGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTgg                                                          <  6:220188/149‑1 (MQ=255)
                                                                                         cATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTa                                                       <  5:324394/149‑1 (MQ=255)
                                                                                             ggCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTg                                                    <  8:193814/148‑1 (MQ=255)
                                                                                              gCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCg                                                  >  2:48083/1‑149 (MQ=255)
                                                                                                      gcgGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTcgccgc                                          >  4:199725/1‑149 (MQ=255)
                                                                                                            cAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGc                                    >  2:302760/1‑149 (MQ=255)
                                                                                                              ggTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCg                                  >  7:281175/1‑149 (MQ=255)
                                                                                                              ggTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCg                                  <  4:39841/149‑1 (MQ=255)
                                                                                                               gtcgtcGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATGGCGGCGACCGCAATAGccag                                 <  1:129058/149‑3 (MQ=255)
                                                                                                                   tcgACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTc                             <  5:294127/149‑1 (MQ=255)
                                                                                                                     gACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCCAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGGCa                            >  7:120633/1‑148 (MQ=255)
                                                                                                                        caccaGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAAc                        <  1:48083/149‑1 (MQ=255)
                                                                                                                         accaGCAGGATATTGGCCTGCGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAAcc                       <  8:137400/149‑1 (MQ=255)
                                                                                                                             gcagGATATTGGCCTGGGACCTGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTTAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACAggag                   >  6:268999/1‑147 (MQ=255)
                                                                                                                                   tatTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCTAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCCGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATa             <  8:297387/149‑1 (MQ=255)
                                                                                                                                       ggCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATAccgc         >  2:212484/1‑149 (MQ=255)
                                                                                                                                        gCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATAccgcc        >  5:208811/1‑149 (MQ=255)
                                                                                                                                        aacgcggCCCAGGCTGCCGCCGCCCGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATAccgcc        >  7:283470/6‑149 (MQ=255)
                                                                                                                                           tGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCgg      >  4:255040/1‑148 (MQ=255)
                                                                                                                                              gCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAAt  <  4:47762/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCGGGCCATGCCGGCAATGTCGCCCACGCCGCGAATCTGGTCGCGGGCAATCGTCATCAGTTGCTCGTAGCGGCGGGTCTCGATCCAGGCATAGGCTGTGTCGCGGGCCAGGTCGTCGACCACCAGCAGGATATTGGCCTGGGCCCGGGCGGCGGCAGCCTGGGCCGCGTTGACGGCGCTGTCGACCTTGCCGAAGTCATAGAGCATTTGCTTGAGCGACAGGTTCAGGGCCTGGCTATTGCGGTCGCCGCCATAGCCGGTGTCATAACCGGTGCGGATACCGCCGGAAAT  >  NC_002947/926011‑926301

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