Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,066,013 G→T 7.9% H1017Q (CAC→CAA)  mdtC ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,066,0130GT 7.9% 192.3 / 5.6 63H1017Q (CAC→CAA) mdtCmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base T (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

GGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAA  >  NC_002947/4065867‑4066156
                                                                                                                                                  |                                                                                                                                               
ggCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTg                                                                                                                                               <  7:29295/149‑1 (MQ=255)
  ccACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgc                                                                                                                                             <  5:154703/149‑1 (MQ=255)
  ccACGCACAACCCGACACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgc                                                                                                                                             >  5:331841/1‑149 (MQ=255)
   cACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgcg                                                                                                                                            >  8:79403/1‑149 (MQ=255)
          aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGTGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                     <  7:234993/149‑1 (MQ=255)
          aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                     <  4:82472/149‑1 (MQ=255)
          aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                     >  7:324045/1‑149 (MQ=255)
          aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                     >  8:39704/1‑149 (MQ=255)
             ccGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                                   >  7:42574/1‑148 (MQ=255)
                gCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcagg                                                                                                                               <  8:115437/149‑1 (MQ=255)
                       gCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggt                                                                                                                        <  3:10956/149‑1 (MQ=255)
                         aGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTaa                                                                                                                      <  8:42574/149‑1 (MQ=255)
                                aGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGGCGGTACGCACGCCGCGACAACGGATGAAACGGTGGCGCAGGCGGTAGAGGTACAGGTAAACCACCg                                                                                                               >  6:13380/1‑149 (MQ=255)
                                  ccacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCACCGTTTCAACCGTTGgcg                                                                                                                                            >  7:334401/1‑118 (MQ=255)
                                  ccacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCACCGTTTCAACCGTTGgcg                                                                                                                                            <  8:334401/118‑1 (MQ=255)
                                  ccacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGc                                                                                                             >  8:259443/1‑149 (MQ=255)
                                        ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAaccac                                                                                                                 <  1:181196/139‑1 (MQ=255)
                                        ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAaccac                                                                                                                 >  2:181196/1‑139 (MQ=255)
                                        ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGta                                                                                                       <  6:331841/149‑1 (MQ=255)
                                          acggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGtat                                                                                                      >  6:9784/1‑148 (MQ=255)
                                            gggTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                             <  1:312125/123‑1 (MQ=255)
                                            gggTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                             >  2:312125/1‑123 (MQ=255)
                                             ggTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAg                                                                                                  <  7:259443/149‑1 (MQ=255)
                                                 aGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGTTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTg                                                                                              <  5:114776/149‑1 (MQ=255)
                                                  gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                     >  7:252461/1‑109 (MQ=255)
                                                  gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                                     <  8:252461/109‑1 (MQ=255)
                                                  gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGa                                                                                             <  7:306224/149‑1 (MQ=255)
                                                          aGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTTGTATAAAGGGTGAGGATCTgg                                                                                     >  4:323658/1‑149 (MQ=255)
                                                           gggCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTTGTATAAAGGGTGAGGATCTGGc                                                                                    <  3:323658/149‑1 (MQ=255)
                                                               aCGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGaa                                                                                >  2:112210/1‑149 (MQ=255)
                                                                 ggTGAAGTGGGGTCTGGGCAAAAGTCATAGGGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAaca                                                                              >  1:29706/1‑149 (MQ=255)
                                                                 ggTGAAGTGGGGTCTGGGCAAAAATCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtcgt                                                                                                        <  5:39114/123‑1 (MQ=255)
                                                                 ggTGAAGTGGGGTCTGGGCAAAAATCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtcgt                                                                                                        >  6:39114/1‑123 (MQ=255)
                                                                         ggggTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCa                                                                      >  2:270769/1‑149 (MQ=255)
                                                                          gggTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCa                                                                      <  7:79403/148‑1 (MQ=255)
                                                                                 ggCAAAAGTCATTGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgatgac                                                              <  2:29706/149‑2 (MQ=255)
                                                                                 ggCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgat                                                                 <  7:39704/146‑1 (MQ=255)
                                                                                  gCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGatgatg                                                             <  8:324045/149‑1 (MQ=255)
                                                                                   cAAAAGTCATAGCGGGCTTTCCAGAGCGGCGTCGCTACGCACTCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATgg                                                            <  4:58306/149‑1 (MQ=255)
                                                                                     aaaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTc                                                          <  5:301027/149‑1 (MQ=255)
                                                                                         tCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCTGCGCCACCGTTTCAACCGTTgg                                                                                                                                              >  8:306854/1‑61 (MQ=21)
                                                                                         tCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCTGCGCCACCGTTTCAACCGTTgg                                                                                                                                              <  7:306854/61‑1 (MQ=21)
                                                                                           aTAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGa                                                                                             >  3:202954/1‑108 (MQ=255)
                                                                                           aTAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGa                                                                                             <  4:202954/108‑1 (MQ=255)
                                                                                             aGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAg                                                  >  8:16315/1‑149 (MQ=255)
                                                                                              gCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGc                                                 <  5:9784/149‑1 (MQ=255)
                                                                                               cGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCg                                                <  7:16315/149‑1 (MQ=255)
                                                                                                          gagCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCa                                       >  1:128278/1‑147 (MQ=255)
                                                                                                                         gcacgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGcc                      >  7:226196/1‑149 (MQ=255)
                                                                                                                         gcacgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGcc                      >  7:201082/1‑149 (MQ=255)
                                                                                                                          cacgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCt                     <  4:314419/149‑1 (MQ=255)
                                                                                                                            cgccggGCAAACGGTGTAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGCGTATCTGGCTCAACACCAGGCCTCCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTgg                   <  1:270769/149‑1 (MQ=255)
                                                                                                                             cctgcgcCACCGTTTCAACCGTTGGCGCAGGCAGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGc                  >  1:43468/4‑149 (MQ=255)
                                                                                                                               cgcgGCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCATGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCt                 <  8:226196/148‑1 (MQ=255)
                                                                                                                                gcgcCAACGGTTGAAACGGTGGCGGAGGCGGTCCAGGTAGAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTca               <  5:13380/149‑1 (MQ=255)
                                                                                                                                    cAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagcag           >  3:250842/1‑149 (MQ=255)
                                                                                                                                     aaCGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagcag           >  7:143063/1‑148 (MQ=255)
                                                                                                                                           tGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTgc                                >  3:271626/1‑121 (MQ=255)
                                                                                                                                           tGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTgc                                <  4:271626/121‑1 (MQ=255)
                                                                                                                                            gAAACGGTGGCGCAGGCCGAACAGGTACACGTAAAACACCGGCGTCTAATAAAGGGTGACGATCTGGCTCAACACCAGGCCACCGATGATGGTCACGCCCACCGGCAGACCCAGTTCTGCGCCTGCCACCTCCCTCAGCAGCAGCGGc    >  8:268635/1‑148 (MQ=255)
                                                                                                                                             aaaCGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCaa  >  1:19001/1‑149 (MQ=255)
                                                                                                                                                cGGTGGCGCAGGCGGTCGAGGTACAGGTAAGCCACCGGCGTGGGATAAAGGGTGAGGATCTGGCTGAACACCAGGCCCCGGATGATGGTCGGGCCCAGCGGCTGACGCCTTTCGGCGCCTTCGGCCTGGGTccgc             >  8:33344/1‑135 (MQ=255)
                                                                                                                                                cGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagc             <  7:33344/135‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAA  >  NC_002947/4065867‑4066156

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