Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,706,254 T→C 5.9% T347A (ACC→GCC)  phaK ← phenylacetic acid‑specific porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,706,2540TC 5.9% 199.2 / 2.5 69T347A (ACC→GCC) phaKphenylacetic acid‑specific porin
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base C (2/2);  total (34/35)
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.27e-01

GGAGCGGAAGGTGACATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACGGGTCGCCGCCGTTGACGAACGGGAAGGCGCTG  >  NC_002947/3706108‑3706402
                                                                                                                                                  |                                                                                                                                                    
ggAGCGGAAGGTGACATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCa                                                                                                                                                    >  6:328520/1‑149 (MQ=255)
   gCGGAAGGTGACATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGa                                                                                                                                                 <  2:92439/149‑1 (MQ=255)
      gAAGGTGACATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAAc                                                                                                                                              >  8:298612/1‑149 (MQ=255)
      gAAGGTGACATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAAc                                                                                                                                              >  8:142766/1‑149 (MQ=255)
         ggTGACATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCa                                                                                                                                           <  2:110173/149‑1 (MQ=255)
             aCATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGcc                                                                                                                                       >  2:328917/1‑149 (MQ=255)
              cATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCg                                                                                                                                      <  7:142766/149‑1 (MQ=255)
               aTTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCgg                                                                                                                                     >  2:46908/1‑149 (MQ=255)
                ttGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTTTAGCGGGGCATGAAACTCAGGCCggg                                                                                                                                    >  7:255578/1‑149 (MQ=255)
                    gcCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATg                                                                                                                                >  1:144767/1‑149 (MQ=255)
                    gcCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACAAACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATAGGTCACGCTGGCAGTTTCCGCGTGGCGACCCTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATg                                                                                                                                >  3:315424/1‑149 (MQ=255)
                      cAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGcc                                                                                                                              >  3:96889/1‑149 (MQ=255)
                      cAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGcc                                                                                                                              >  3:371823/1‑149 (MQ=255)
                         cgcAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCAt                                                                                                                           <  7:78869/149‑1 (MQ=255)
                           cAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATgg                                                                                                                         <  4:305678/149‑1 (MQ=255)
                            aGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATgg                                                                                                                         >  1:35226/1‑148 (MQ=255)
                              gCTGACGTCCTTGAACGTGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCaa                                                                                                                      <  5:218981/149‑1 (MQ=255)
                                 gACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCTCGGTGGCAGTTTCCGTGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACg                                                                                                                    >  1:240879/1‑148 (MQ=255)
                                    gTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagt                                                                                                                <  2:144767/149‑1 (MQ=255)
                                     tCCTTGATCGGGCCACTCTGGATGACATACGTGTTGTGGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGAGCGTCGGTATAGCGGGTCATGAATCTCAGGCCGGGGATGCCCATGGCAACGAagtc                                                                                                               <  4:315424/149‑1 (MQ=255)
                                          gAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtcgta                                                                                                            <  7:298612/147‑1 (MQ=255)
                                          gAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtc                                                                                                               <  1:151410/144‑1 (MQ=255)
                                           aaCGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtcgtagtc                                                                                                         >  4:93719/1‑149 (MQ=255)
                                            aCGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtcgtagtcg                                                                                                        >  6:128288/1‑149 (MQ=255)
                                               ggCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACg                                                                                                                    <  8:211627/134‑1 (MQ=255)
                                               ggCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACg                                                                                                                    >  7:211627/1‑134 (MQ=255)
                                               ggCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGtcgtagtcgtag                                                                                                     >  6:30047/1‑149 (MQ=255)
                                               ggCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGtcgtagtcgtag                                                                                                     >  7:297419/1‑149 (MQ=255)
                                                gCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAgc                                                                                                    <  4:96889/149‑1 (MQ=255)
                                                       tGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTg                                                                                             >  3:316433/1‑149 (MQ=255)
                                                         gATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAACCTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTg                                                                                             <  7:122222/147‑1 (MQ=255)
                                                         gATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGcc                                                                                           <  1:46908/149‑1 (MQ=255)
                                                          aTGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGcc                                                                                           >  5:385533/1‑148 (MQ=255)
                                                           tGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtc                                                                                                               >  2:209319/1‑127 (MQ=255)
                                                           tGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtc                                                                                                               <  1:209319/127‑1 (MQ=255)
                                                             aCATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAaga                                                                                       <  2:314290/149‑1 (MQ=255)
                                                              cATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCTCGGTGGCAGTTTCCGTGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAaga                                                                                       <  2:240879/148‑1 (MQ=255)
                                                               aTACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtcgtagt                                                                                                          >  2:301795/1‑128 (MQ=255)
                                                               aTACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtcgtagt                                                                                                          <  1:301795/128‑1 (MQ=255)
                                                                   gTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAg                                                                                 >  3:138300/1‑149 (MQ=255)
                                                                          cGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCgg                                                                                                                                     <  2:35722/90‑1 (MQ=255)
                                                                            gtgtCGCGCTCGCGCTCACGGCCATTGGTCTCGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGcgcg                                                                        >  1:35368/1‑149 (MQ=255)
                                                                            gtgtCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGcgcg                                                                        >  1:66394/1‑149 (MQ=255)
                                                                             tgtCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCgg                                                                       <  5:128288/149‑1 (MQ=255)
                                                                             tgtCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCgg                                                                       <  5:30047/149‑1 (MQ=255)
                                                                             tgtCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCgg                                                                       <  8:297419/149‑1 (MQ=255)
                                                                               tcgcgctcgcgcTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtc                                                                                                               >  4:123755/1‑107 (MQ=255)
                                                                               tcgcgctcgcgcTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAagtc                                                                                                               <  3:123755/107‑1 (MQ=255)
                                                                                   gctcgcgctcACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAg                                                                 <  4:138300/149‑1 (MQ=255)
                                                                                           tcACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGtgttg                                                           >  8:350833/1‑147 (MQ=255)
                                                                                           tcACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGCCGAGGGTCAAGGtgttg                                                           >  7:225920/1‑147 (MQ=255)
                                                                                                        tCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGtt                                             >  3:69862/1‑148 (MQ=255)
                                                                                                        tCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGCCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAgtcgtagtcgta                                                                                                      >  5:319445/1‑91 (MQ=255)
                                                                                                        tCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGCCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAgtcgtagtcgta                                                                                                      <  6:314177/91‑1 (MQ=255)
                                                                                                        tCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGCCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAgtcgtagtcgta                                                                                                      <  6:319445/91‑1 (MQ=255)
                                                                                                        tCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGCCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAgtcgtagtcgta                                                                                                      >  5:314177/1‑91 (MQ=255)
                                                                                                          aCGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGGCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGAc                                          <  2:35368/149‑1 (MQ=255)
                                                                                                           cGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACg                                         <  4:69862/149‑1 (MQ=255)
                                                                                                             gTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGtgttg                                                           >  6:24603/1‑129 (MQ=255)
                                                                                                             gTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGtgttg                                                           <  5:24603/129‑1 (MQ=255)
                                                                                                             gTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACTGt                                       <  4:316433/149‑1 (MQ=255)
                                                                                                                gCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGaa                                    <  5:328520/149‑1 (MQ=255)
                                                                                                                   ggttCCGCGTGGCGACCGTCGGTATAGCGGGGCATGCACCTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACgg                                 <  8:255578/147‑1 (MQ=255)
                                                                                                                             ggCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGa                                                    <  2:81051/120‑1 (MQ=255)
                                                                                                                             ggCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGa                                                    >  1:81051/1‑120 (MQ=255)
                                                                                                                                gACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACGGGTCGCCGCCGtt                     >  1:155655/1‑148 (MQ=255)
                                                                                                                                 aCCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACGGGTCGCCGCCGTTGa                   <  2:66394/149‑1 (MQ=255)
                                                                                                                                            aGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACGGGTCGCCGCCGTTGACGAACGGGAAg        >  7:278123/1‑149 (MQ=255)
                                                                                                                                                 gTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACGGGTCGCCGCCGTTGACGAACGGGAAGGCGCt   <  4:371823/149‑1 (MQ=255)
                                                                                                                                                  tCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACGGGTCGCCGCCGTTGACGAACGGGAAGGCGCTg  <  1:328917/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGAGCGGAAGGTGACATTGCGCCAGCGCAGGCTGACGTCCTTGAACGGGCCACTCTGGATGACATACGTGATGTCGGTGTCGCGCTCGCGCTCACGGCCATTGGTCACGGTGGCAGTTTCCGCGTGGCGACCGTCGGTATAGCGGGTCATGAAACTCAGGCCGGGGATGCCCATGGCAACGAAGTCGTAGTCGTAGCGCACTTGCCAAGAGTCCAGCCCGGCGCGGGTGAAGGTGTTGTAGGTGACCAGGTTGACGGTGAACGGGTCGCCGCCGTTGACGAACGGGAAGGCGCTG  >  NC_002947/3706108‑3706402

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: