Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,249 T→C 6.3% T392A (ACC→GCC)  kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,2490TC 6.3% 224.5 / 4.0 79T392A (ACC→GCC) kgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base C (3/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAA  >  NC_002947/1594101‑1594393
                                                                                                                                                         |                                                                                                                                                
cGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGt                                                                                                                                                  <  4:207805/149‑1 (MQ=255)
     cTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTgccgc                                                                                                                                             >  8:220035/1‑149 (MQ=255)
      tCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTgccgcc                                                                                                                                            >  8:280807/1‑149 (MQ=255)
           ggTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGccgccg                                                                                                                                           >  3:379152/1‑145 (MQ=255)
             tAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGccgccg                                                                                                                                           >  6:198909/1‑143 (MQ=255)
               aTCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCCGTGACGTACCAGGAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGa                                                                                                                                   <  1:72555/149‑1 (MQ=255)
                tCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaatacgt                                                                                                                                    <  4:346213/147‑6 (MQ=255)
                tCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaatacgt                                                                                                                                    >  3:346213/1‑142 (MQ=255)
                          ctgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAATTAGCGCTCTCCATCCGAGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGt                                                                                                                             >  1:122353/1‑148 (MQ=255)
                          ctgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                         <  1:383901/148‑1 (MQ=255)
                                ggTGTCTTTCATGGTCGCGTAAACCAGCAGCGTGGAGGGGATACAGTCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATa                                                                                                                    <  4:124391/147‑1 (MQ=255)
                                ggTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACg                                                                                                                  >  5:30762/1‑149 (MQ=255)
                                  tgtCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCTCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATa                                                                                                                    <  8:229977/145‑1 (MQ=255)
                                  tgtCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATa                                                                                                                    >  7:229977/1‑145 (MQ=255)
                                      tttCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAgg                                                                                                             >  1:91796/1‑149 (MQ=255)
                                        tCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGc                                                                                                            <  7:12134/147‑1 (MQ=255)
                                                gcgTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAg                                                                                                  >  7:142851/1‑149 (MQ=255)
                                                 cgTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAgg                                                                                                 <  1:305428/149‑1 (MQ=255)
                                                   tAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGgc                                                                                               <  2:91796/149‑1 (MQ=255)
                                                   tAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaatac                                                                                                                                      <  5:53136/110‑4 (MQ=255)
                                                   tAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaatac                                                                                                                                      >  6:53136/1‑107 (MQ=255)
                                                       ccAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGcgcg                                                                                            >  3:239996/1‑148 (MQ=255)
                                                        cagcagCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGaa                                                                                          <  8:108076/149‑1 (MQ=255)
                                                        cagcagCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGaa                                                                                          >  4:118658/1‑149 (MQ=255)
                                                          gcagcGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCCTACGGCAGGCCAACGCCCAGGGCGCGAAtt                                                                                        >  2:123131/1‑149 (MQ=255)
                                                                 gCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACCGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAg                                                                                 >  5:169751/1‑149 (MQ=255)
                                                                    ggCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAggg                                                                               >  7:261400/1‑148 (MQ=255)
                                                                           cAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                     >  4:179812/1‑135 (MQ=255)
                                                                           cAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                     <  3:179812/135‑1 (MQ=255)
                                                                            aGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCg                                                                      <  1:51744/149‑1 (MQ=255)
                                                                              gCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGc                                                                    >  3:307424/1‑149 (MQ=255)
                                                                                 gTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                 <  4:239996/149‑1 (MQ=255)
                                                                                 gTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCCGCGCTCTAGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGTCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                 <  7:300170/149‑1 (MQ=255)
                                                                                       gTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                         >  6:76821/1‑88 (MQ=255)
                                                                                       gTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                         <  5:76821/88‑1 (MQ=255)
                                                                                       gTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCCCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATAGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCTGTAGGCACCAGCTCGGCTTATACCa                                                            >  4:263092/1‑149 (MQ=255)
                                                                                        tACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCAc                                                           <  7:280807/149‑1 (MQ=255)
                                                                                         aCCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACAGCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAg                                                                                 <  5:202116/126‑1 (MQ=255)
                                                                                         aCCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAg                                                                                 >  6:202116/1‑126 (MQ=255)
                                                                                           cAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGTCCACGTATTCGGCGGTGCCGCAGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCGAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGc                                                        >  5:100659/1‑149 (MQ=255)
                                                                                           cAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCgg                                                                     >  8:29988/1‑136 (MQ=255)
                                                                                           cAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCgg                                                                     <  7:29988/136‑1 (MQ=255)
                                                                                           cAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGc                                                        >  5:254800/1‑149 (MQ=255)
                                                                                             gTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                      >  2:387231/1‑149 (MQ=255)
                                                                                                gAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGa                                                   <  4:307424/149‑1 (MQ=255)
                                                                                                   aaggCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                      <  5:277273/143‑1 (MQ=255)
                                                                                                   aaggCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                      >  6:277273/4‑146 (MQ=255)
                                                                                                      gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATc                                                 >  7:23820/1‑148 (MQ=255)
                                                                                                      gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATc                                                 >  6:95194/1‑148 (MQ=255)
                                                                                                      gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCg                                                <  2:122353/149‑1 (MQ=255)
                                                                                                       ccGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGa                                               <  4:379152/149‑1 (MQ=255)
                                                                                                        cGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAg                                              <  8:261400/149‑1 (MQ=255)
                                                                                                         gCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAgg                                             >  7:87926/1‑149 (MQ=255)
                                                                                                          ctctCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGt                                            <  5:95194/149‑1 (MQ=255)
                                                                                                                aTCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGcc                                      <  8:87926/149‑1 (MQ=255)
                                                                                                                aTCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGcc                                      >  7:299000/1‑149 (MQ=255)
                                                                                                                aTCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGCCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGAGTAGATCTAGGTATAAcc                                      >  6:347275/1‑149 (MQ=255)
                                                                                                                 tCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCg                                     <  6:30762/149‑1 (MQ=255)
                                                                                                                   ccGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGc                                    <  7:220035/148‑1 (MQ=255)
                                                                                                                     cGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCtt                                  <  5:198909/149‑1 (MQ=255)
                                                                                                                     cGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCt                                   >  6:119421/1‑148 (MQ=255)
                                                                                                                           tctTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGATTACgat                             >  1:266537/1‑148 (MQ=255)
                                                                                                                             ttGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACgatgat                          >  3:279801/1‑149 (MQ=255)
                                                                                                                             ttGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACgatgat                          >  2:44816/1‑149 (MQ=255)
                                                                                                                             ttGAACCACAGGGCCACCTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGAAGGCCAACCCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTATCACCACGGCCTAGATCGAGCTATAGCCGCTAACgatga                           >  6:3887/1‑148 (MQ=255)
                                                                                                                                 aCCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGt                      <  8:23820/149‑1 (MQ=255)
                                                                                                                                  ccACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTg                     <  6:169751/149‑1 (MQ=255)
                                                                                                                                   cacaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGtat                                          >  7:284396/1‑127 (MQ=255)
                                                                                                                                   cacaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGtat                                          <  8:284396/127‑1 (MQ=255)
                                                                                                                                    acaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCg                   <  3:118658/149‑1 (MQ=255)
                                                                                                                                    acaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCg                   >  2:98172/1‑149 (MQ=255)
                                                                                                                                       gggCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGcc                <  8:142851/149‑1 (MQ=255)
                                                                                                                                          ccACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATg             <  6:100659/149‑1 (MQ=255)
                                                                                                                                           cACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGa            >  3:9205/1‑149 (MQ=255)
                                                                                                                                           cACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGa            >  1:72195/1‑149 (MQ=255)
                                                                                                                                           cACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGa            >  1:147496/1‑149 (MQ=255)
                                                                                                                                             gatcttCGGCGGCACCGCCGAATATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCCTGATc          >  4:303149/5‑149 (MQ=255)
                                                                                                                                               tATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAg        <  2:147496/149‑1 (MQ=255)
                                                                                                                                                aTTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAg        >  8:272856/1‑148 (MQ=255)
                                                                                                                                                     gcggTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGgaagaa  <  6:254800/149‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                
CGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAA‑‑‑GCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAA  >  NC_002947/1594101‑1594393

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: