Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,645,887 G→C 6.5% S490R (AGC→AGG)  gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,645,8870GC 6.5% 195.2 / 2.9 61S490R (AGC→AGG) gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTG  >  NC_002947/1645739‑1646033
                                                                                                                                                    |                                                                                                                                                  
cTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCggctgg                                                                                                                                                    <  4:196765/149‑1 (MQ=255)
       ttGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGa                                                                                                                                             >  3:193501/1‑149 (MQ=255)
            gggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTcc                                                                                                                                        <  1:79176/149‑1 (MQ=255)
            gggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTc                                                                                                                                         >  3:148169/1‑148 (MQ=255)
             ggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCa                                                                                                                                       >  1:245918/1‑149 (MQ=255)
              gtggggATCTCGCGGAGCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCAt                                                                                                                                      <  6:105194/149‑1 (MQ=255)
                ggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGt                                                                                                                                    <  5:79905/149‑1 (MQ=255)
                 gggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGGTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTc                                                                                                                                   >  8:267667/1‑149 (MQ=255)
                  ggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTcc                                                                                                                                  >  2:100157/1‑149 (MQ=255)
                       tcGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAg                                                                                                                             >  8:5272/1‑149 (MQ=255)
                       tcGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCTGCCGAccc                                                                                                                                          <  1:271309/136‑3 (MQ=255)
                       tcGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCTGCCGAccc                                                                                                                                          >  2:271309/1‑134 (MQ=255)
                       tcGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccct                                                                                                                                         >  1:248161/1‑134 (MQ=255)
                       tcGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccct                                                                                                                                         <  2:248161/137‑4 (MQ=255)
                       tcGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACACGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAg                                                                                                                             >  7:208750/1‑149 (MQ=255)
                             aTCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgtg                                                                                                                       <  2:314512/149‑1 (MQ=255)
                                 gCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGtggtgg                                                                                                                   >  3:170371/1‑149 (MQ=255)
                                 gCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGtggtgg                                                                                                                   <  1:165183/149‑1 (MQ=255)
                                 gCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGtggtgg                                                                                                                   >  2:144377/1‑149 (MQ=255)
                                 gCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGtggtgg                                                                                                                   >  2:165183/1‑149 (MQ=255)
                                  cACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCCCAGCCTCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGAAGCCGACGTCCATGTCGCACAGGTCgtggtgcg                                                                                                                  >  5:324098/1‑147 (MQ=255)
                                   aCGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTggtgggt                                                                                                                 <  2:172301/149‑1 (MQ=255)
                                       gtgtcGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCCCCGGTGTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                             <  8:322614/145‑1 (MQ=255)
                                          gCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                             >  5:302435/1‑146 (MQ=255)
                                          gCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                             <  6:302435/146‑1 (MQ=255)
                                             tACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAg                                                                                                       <  7:267667/149‑1 (MQ=255)
                                                cGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTc                                                                                                    >  3:277939/1‑149 (MQ=255)
                                                cGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTc                                                                                                    <  5:284656/149‑1 (MQ=255)
                                                  tCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTTGTGGTGGGTGAACTGATAGTTc                                                                                                    >  2:219848/1‑147 (MQ=255)
                                                  tCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTc                                                                                                    >  6:134850/1‑147 (MQ=255)
                                                        cggTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGc                                                                                            <  2:245918/149‑1 (MQ=255)
                                                         ggTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGcc                                                                                           >  5:146854/1‑149 (MQ=255)
                                                         ggTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGGGTACCAGGGTCGGCTTGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTTAACTGATAGTTCCAGCGTGcc                                                                                           <  2:32901/149‑1 (MQ=255)
                                                           tCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCtt                                                                                         <  8:35468/149‑1 (MQ=255)
                                                             gAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGc                                                                                       >  4:40185/1‑149 (MQ=255)
                                                               tcACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCg                                                                                     <  8:208750/148‑1 (MQ=255)
                                                                  cGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTg                                                                                  <  1:144377/149‑1 (MQ=255)
                                                                   gTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTgg                                                                                 >  5:352059/1‑149 (MQ=255)
                                                                         gCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAgg                                                                            >  7:164024/1‑148 (MQ=255)
                                                                               ttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAggg                                                                     >  4:39171/1‑149 (MQ=255)
                                                                               ttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAggg                                                                     >  7:305033/1‑149 (MQ=255)
                                                                               ttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAgg                                                                      >  7:250014/1‑148 (MQ=255)
                                                                                 gcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAgg                                                                      <  8:250014/146‑1 (MQ=255)
                                                                                   ttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTGGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGGGAc                                                                 >  8:60514/1‑149 (MQ=255)
                                                                                     ggTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGAc                                                                 <  5:134850/147‑1 (MQ=255)
                                                                                          gTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGt                                                                           >  1:78579/1‑132 (MQ=255)
                                                                                          gTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGt                                                                           <  2:78579/132‑1 (MQ=255)
                                                                                             cGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGgcg                                                       <  7:5272/149‑1 (MQ=255)
                                                                                             cGATGATCGCCGGTTTCACACCATCGGCGGTGTTCAGGTGGACCAGGGGCGGCGGGCGGCCGACGGCCATGGCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGgc                                                        <  5:305104/148‑1 (MQ=255)
                                                                                              gatgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGa                                                                  <  4:35496/137‑1 (MQ=255)
                                                                                              gatgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGa                                                                  >  3:35496/1‑137 (MQ=255)
                                                                                               atgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCt                                                     >  2:355790/1‑149 (MQ=255)
                                                                                                      ccGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTtc                                              <  8:164024/149‑1 (MQ=255)
                                                                                                              acacCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGgcgc                                       <  4:170371/148‑1 (MQ=255)
                                                                                                              acacCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGcc                                      >  5:216515/1‑149 (MQ=255)
                                                                                                               cacCATCGTCGGTTTTCAGGTGTACCAGGGTCAGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGccc                                     >  4:87712/1‑149 (MQ=255)
                                                                                                                      gtcgGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGc                              <  4:277939/149‑1 (MQ=255)
                                                                                                                         ggTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGt                           <  3:39171/149‑1 (MQ=255)
                                                                                                                           ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGccc                                     >  3:257215/1‑137 (MQ=255)
                                                                                                                           ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGccc                                     <  4:257215/137‑1 (MQ=255)
                                                                                                                           ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCg                         <  4:193501/149‑1 (MQ=255)
                                                                                                                                       ccAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTg              <  4:148169/148‑1 (MQ=255)
                                                                                                                                                gctggctgCCGACTTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTc    <  8:305033/149‑1 (MQ=255)
                                                                                                                                                    gctgCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTg  <  1:355790/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTG  >  NC_002947/1645739‑1646033

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