Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,547,134 T→A 6.0% Y71F (TAC→TTC)  pydP ← NCS1 family transporter PydP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,547,1340TA 6.0% 210.2 / 3.1 67Y71F (TAC→TTC) pydPNCS1 family transporter PydP
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCGAACGACGAACGCAGCAGCACCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGTCGGTGCTATGTCGGGGTTGTAGCGCGGGCTGTCGAGG  >  NC_002947/4546996‑4547277
                                                                                                                                          |                                                                                                                                                
ccGTACGACGCACGCAGCAGCACCGGAACCGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGa                                                                                                                                         <  8:135629/148‑1 (MQ=255)
ccGAACGACTAACGCAGCAGCACCTGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGa                                                                                                                                         <  8:140956/148‑1 (MQ=255)
ccGAACGACGAACGCAGCAGCACCGGAAACGGGATGCCGTACTTGCTGCCGGGAAAGGCATTGAGGCTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGa                                                                                                                                         <  2:181855/148‑1 (MQ=255)
  gAACGACGAACGCAGCAGCACCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCa                                                                                                                                      <  5:186853/149‑1 (MQ=255)
    acgacgAACGCAGCAGCACCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCAGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCACTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCAc                                                                                                                                     <  2:63290/148‑1 (MQ=255)
         gaACGCAGCAGCACCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACgccgcc                                                                                                                               >  6:217584/1‑149 (MQ=255)
           aCGCAGCAGCACCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGa                                                                                                                             <  6:375246/149‑1 (MQ=255)
           aCGCAGCAGCACCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGCGTGCTCaccgccta                                                                                                                             <  7:243084/149‑3 (MQ=255)
             gcagcagcaCCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAgg                                                                                                                           <  2:324856/149‑1 (MQ=255)
             gcagcagcaCCCGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAgg                                                                                                                           <  8:316937/149‑1 (MQ=255)
               agcagcaCCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGgt                                                                                                                         >  6:219800/1‑149 (MQ=255)
                 cagcaCCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGgtgt                                                                                                                       <  4:376288/149‑1 (MQ=255)
                  agcaCCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTa                                                                                                                      <  5:23026/149‑1 (MQ=255)
                     aCCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGt                                                                                                                   <  8:166339/149‑1 (MQ=255)
                      ccGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑g                                                                                                                 >  2:159071/1‑149 (MQ=255)
                               ggATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCg                                                                                                                                             >  6:367935/1‑113 (MQ=255)
                               ggATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCg                                                                                                                                             <  5:367935/113‑1 (MQ=255)
                                 aTGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGAt                                                                                                      <  3:115143/149‑1 (MQ=255)
                                 aTGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGAt                                                                                                      <  3:121792/149‑1 (MQ=255)
                                          ttGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGcc                                                                                             >  7:26587/1‑149 (MQ=255)
                                            ggTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGccca                                                                                           >  3:172212/1‑149 (MQ=255)
                                              tGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGcccacc                                                                                         >  6:146204/1‑149 (MQ=255)
                                                  gggAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGcccaccca                                                                                       >  1:326637/1‑147 (MQ=255)
                                                  gggAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGAAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGcccaccca                                                                                       >  8:91795/1‑147 (MQ=255)
                                                  gggAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGAAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCcaca                                                                                     >  8:86411/1‑149 (MQ=255)
                                                     aaaGGCATTGCGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTGGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAggg                                                                                  <  8:71640/149‑1 (MQ=255)
                                                       aGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCg                                                                                >  7:279614/1‑149 (MQ=255)
                                                             tGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGAt                                                                           <  8:26587/148‑1 (MQ=255)
                                                               aGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCg                                                                                >  7:289042/1‑141 (MQ=255)
                                                               aGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCg                                                                                <  8:289042/141‑1 (MQ=255)
                                                                  gTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCa                                                                     >  5:193997/1‑149 (MQ=255)
                                                                   tAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCAt                                                                    >  3:162747/1‑149 (MQ=255)
                                                                     aGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCAttt                                                                  >  6:160922/1‑149 (MQ=255)
                                                                      gCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTg                                                                 >  4:35784/1‑149 (MQ=255)
                                                                        ggTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGtt                                                               <  4:154391/149‑1 (MQ=255)
                                                                         gTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCTGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTc                                                              <  4:172212/149‑1 (MQ=255)
                                                                         gTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTc                                                              <  1:64846/149‑1 (MQ=255)
                                                                          tATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTcc                                                             <  5:219800/149‑1 (MQ=255)
                                                                            tCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAg                                                           <  5:217584/149‑1 (MQ=255)
                                                                              aGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGt                                                         >  2:198933/1‑149 (MQ=255)
                                                                               gCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTg                                                        >  1:44777/1‑149 (MQ=255)
                                                                               gCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTCATGTGCCATTTGTTCCAGGTg                                                        >  8:198240/1‑149 (MQ=255)
                                                                                cACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTgc                                                       >  8:36956/1‑149 (MQ=255)
                                                                                cACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTgc                                                       >  2:302191/1‑149 (MQ=255)
                                                                                 aCGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTgcg                                                      <  8:279614/149‑1 (MQ=255)
                                                                                    aTTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCAAGGTgcacca                                                   <  2:326637/149‑5 (MQ=255)
                                                                                       aCGTTGGCCAGCAAAATCGCCAGCAGAGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGt                                                 <  1:159071/148‑1 (MQ=255)
                                                                                         gTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGAAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCa                                              <  7:86411/149‑1 (MQ=255)
                                                                                         gTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGAAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCa                                              <  7:91795/149‑1 (MQ=255)
                                                                                               cAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGAtgt                                                                         >  7:309144/1‑116 (MQ=255)
                                                                                               cAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGAtgt                                                                         <  8:309144/116‑1 (MQ=255)
                                                                                                aGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCACCAGCGCACCtggaa                                       <  4:300731/149‑5 (MQ=255)
                                                                                                  cAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGTCg                                     <  4:162747/149‑1 (MQ=255)
                                                                                                    aaaTCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGTCGGt                                   <  5:146204/149‑1 (MQ=255)
                                                                                                       tCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGt                                                         >  7:52265/1‑124 (MQ=255)
                                                                                                       tCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGt                                                         <  8:52265/124‑1 (MQ=255)
                                                                                                                 cgcTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTg                                                   <  6:333964/120‑1 (MQ=255)
                                                                                                                 cgcTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTGCCAGGTGCGCTg                                                   >  5:333964/1‑120 (MQ=255)
                                                                                                                  gcTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGTCGGTGCTATGTCGGGGtt                     >  8:68933/1‑149 (MQ=255)
                                                                                                                   cttagacAACGCTGAGGACGAAGTAGGCGGTGAGCACGCAGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGCGCGGTGAAGTGCAAATTGATCCAGGTGCGCTGCTGCACCTTGGTCgg                                    >  2:349487/7‑133 (MQ=255)
                                                                                                                   cTTCGCCATCGCTGATGCCGAAGTAGGCGGTTAGGAGGGCGGCGTTTGTGTAGTT‑GTGTACGCATATGGACATGCACACCCACAGGGCGGTTTTGTGCCATTTGTTCCAGGTGCGGTGGTGCACCTTGGTCgg                                    <  1:349487/133‑1 (MQ=255)
                                                                                                                         cAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGt                                       <  7:84860/124‑1 (MQ=255)
                                                                                                                         cAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGt                                       >  8:84860/1‑124 (MQ=255)
                                                                                                                           aCGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCAGAGGGTGTAGGT‑GGGTACTCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTACAGGTGCGCTGGTGCACCTTGGTCGTTTCTATGTCGGGGTTGTATCGCgg            >  6:387921/1‑149 (MQ=255)
                                                                                                                             gCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGTCGGTGCTATGTCGGGGTTGTAGCGCGGGc          >  3:115696/1‑149 (MQ=255)
                                                                                                                                    ccGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGTCGGTGCTATGTCGGGGTTGTAGCGCGGGCTGTCGAg   <  1:302191/149‑1 (MQ=255)
                                                                                                                                       aaGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGTGGGGTACGAAGATGGACATGCCAACCAACAGGGCGGTGATGTGCAATTTGTTCAAGGTCGGATGTTGAACATTGTTCGTTCGTATGACGGGGTTGAAGCGCGGGATGACAAgg  >  1:64248/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                
CCGAACGACGAACGCAGCAGCACCGGAAACGGGATGCCGTACTTGGTGCCGGGAAAGGCATTGAGGGTAAGCGGTATCAGCACGATTACGTTGGCCAGCAAAATCGCCAGCAGCGCTTCGCCAACGCTGAGGCCGAAGTAGGCGGTGAGCACGCCGCCGAGGGTGTAGGT‑GGGTACGCAGATGGACATGCCCACCCACAGGGCGGTGATGTGCCATTTGTTCCAGGTGCGCTGGTGCACCTTGGTCGGTGCTATGTCGGGGTTGTAGCGCGGGCTGTCGAGG  >  NC_002947/4546996‑4547277

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