Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,250 G→A 10.5% G391G (GGC→GGT)  kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,2500GA10.5% 178.1 / 8.8 67G391G (GGC→GGT) kgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base A (5/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACG  >  NC_002947/1594109‑1594395
                                                                                                                                             |                                                                                                                                                 
aGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGccgccg                                                                                                                                             >  1:343338/1‑148 (MQ=255)
   ggTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGa                                                                                                                                         <  8:205397/149‑1 (MQ=255)
    gTAATCGGCGAGTGCTTCTTGGTGGCTTTCATGGTCGCGGAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTTCCAGTCGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGa                                                                                                                                         <  8:210851/148‑1 (MQ=255)
     tAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATc                                                                                                                                       <  7:42855/149‑1 (MQ=255)
      aaTCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaatacgt                                                                                                                                      >  8:254804/1‑144 (MQ=255)
       aTCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGa                                                                                                                                     <  5:24799/149‑1 (MQ=255)
         ccGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGaca                                                                                                                                   >  3:112245/1‑149 (MQ=255)
         ccGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaata                                                                                                                                         >  6:144268/1‑141 (MQ=255)
         ccGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaata                                                                                                                                         <  5:144268/143‑3 (MQ=255)
          cgcgAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGacac                                                                                                                                  >  7:274015/1‑149 (MQ=255)
                 gctgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                           <  4:79988/149‑1 (MQ=255)
                  ctgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                           <  8:67992/148‑1 (MQ=255)
                  ctgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                           <  4:115502/148‑1 (MQ=255)
                   tgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAggg                                                                                                                          >  8:284990/1‑148 (MQ=255)
                   tgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                           <  7:221651/147‑1 (MQ=255)
                   tgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                           >  8:221651/1‑147 (MQ=255)
                     ctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCAt                                                                                                                       <  8:274015/149‑1 (MQ=255)
                              tttCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGc                                                                                                              >  6:302654/1‑149 (MQ=255)
                              tttCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGc                                                                                                              >  2:352864/1‑149 (MQ=255)
                                                    agcGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTc                                                                                        <  7:284990/149‑1 (MQ=255)
                                                    agcGAGCAGGCGATACAGGCGGTGACGTACCAGTACAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTc                                                                                        <  6:23500/149‑1 (MQ=255)
                                                     gcGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                       <  4:8901/149‑1 (MQ=255)
                                                        aGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCgg                                                                                      >  4:59646/1‑147 (MQ=255)
                                                          cAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAgg                                                                                  >  5:335886/1‑149 (MQ=255)
                                                            ggCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTAGGTAggg                                                                                 >  8:332680/1‑148 (MQ=255)
                                                              cGATACAGGCGGTGAGGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGTATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAAc                                                                              <  1:258373/149‑1 (MQ=255)
                                                               gATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCTGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGAAGGCCAACGCCCAGGGCGCGAATTTAGGTAGGGaa                                                                               >  4:204091/1‑147 (MQ=255)
                                                                  aCAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCt                                                                          >  5:220602/1‑149 (MQ=255)
                                                                   cAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAATGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTc                                                                         <  1:49965/149‑1 (MQ=255)
                                                                         gTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                   <  7:332680/149‑1 (MQ=255)
                                                                              gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCCACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCa                                                              >  6:291457/1‑149 (MQ=255)
                                                                              gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAAttt                                                                                         <  6:34221/122‑1 (MQ=255)
                                                                              gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAAttt                                                                                         >  5:34221/1‑122 (MQ=255)
                                                                              gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCa                                                              >  5:84883/1‑149 (MQ=255)
                                                                               tACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCAc                                                             >  4:193370/1‑149 (MQ=255)
                                                                                aCCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACg                                                            >  3:134135/1‑149 (MQ=255)
                                                                                    gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                        <  3:95279/149‑1 (MQ=255)
                                                                                     tAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACg                                                                                                         >  2:227992/1‑99 (MQ=255)
                                                                                     tAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACg                                                                                                         <  1:227992/99‑1 (MQ=255)
                                                                                      aGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGtt                                                       <  6:220602/148‑1 (MQ=255)
                                                                                        aaGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGAt                                                    >  3:133919/1‑149 (MQ=255)
                                                                                              ctctCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGa                                                     <  5:114738/142‑1 (MQ=255)
                                                                                              ctctCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGt                                              <  2:36691/149‑1 (MQ=255)
                                                                                              ctcaCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGa                                                     >  6:114738/1‑142 (MQ=255)
                                                                                                   cATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGc                                         >  7:287873/1‑149 (MQ=255)
                                                                                                       ccGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTAAGGGCATACGGCAGTCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGc                                      <  3:204091/148‑1 (MQ=255)
                                                                                                            gcTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACga                                >  1:260348/1‑149 (MQ=255)
                                                                                                            gcTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACga                                >  7:251268/1‑149 (MQ=255)
                                                                                                              tctTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACgat                               >  2:88728/1‑148 (MQ=255)
                                                                                                                  gAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                       >  7:239572/1‑88 (MQ=255)
                                                                                                                  gAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                       <  8:239572/88‑1 (MQ=255)
                                                                                                                   aaCCACAGGGCCACGTATTCGGCGGTGCCACCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTa                                   >  1:301491/1‑139 (MQ=255)
                                                                                                                      cacaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGc                      >  5:46773/1‑149 (MQ=255)
                                                                                                                       acaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCg                     <  4:134135/149‑1 (MQ=255)
                                                                                                                            gCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCAt                >  7:6377/1‑149 (MQ=255)
                                                                                                                                gatcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGAt             >  8:97233/5‑148 (MQ=255)
                                                                                                                                 gTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCa           <  4:133919/149‑1 (MQ=255)
                                                                                                                                 atcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGAt             >  5:20669/4‑147 (MQ=255)
                                                                                                                                    ttCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCAGGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGga        >  1:305212/1‑149 (MQ=255)
                                                                                                                                    ttCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAgg         <  4:112245/148‑1 (MQ=255)
                                                                                                                                    ttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAgg         >  8:296216/1‑148 (MQ=255)
                                                                                                                                      cggcggTTCCGGCGAAGTTCTAGAGCGTCAGGTTATACGGCTGGAGAACGCTGTGGGCTCGTATTTCGGTAGGGTATAGCTCGTCTTTGATCACGGCGTTGATCGAGGTATATCCGCTTACTATGATGAGTGCTGCCATTATCAGgaa       <  4:154855/148‑1 (MQ=255)
                                                                                                                                      cggcggCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAg          <  7:97233/145‑1 (MQ=255)
                                                                                                                                       ggcggTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGgaag      >  1:41188/1‑148 (MQ=255)
                                                                                                                                        gcggTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGgaagaa    >  5:315324/1‑149 (MQ=255)
                                                                                                                                        gcggTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGAGTAGATCTAGGTATACCCGCTTACAATAATAAGTGAGGCCATGATCAGgaagaa    >  8:24252/1‑149 (MQ=255)
                                                                                                                                          ggTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACg  <  3:59646/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
AGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACG  >  NC_002947/1594109‑1594395

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