Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,980,282 C→G 8.4% D117H (GAC→CAC)  cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,980,2820CG 8.4% 198.2 / 5.1 65D117H (GAC→CAC) cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base G (2/3);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTC  >  NC_002947/3980136‑3980426
                                                                                                                                                  |                                                                                                                                                
aTAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCAt                                                                                                                                                <  5:21269/149‑1 (MQ=255)
 tAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCATAGCCGCTCGGCCTGCTCACCGGGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATc                                                                                                                                               <  7:240663/149‑1 (MQ=255)
 tAGTCGCGGCCCAGCCAATGGCTGGCCAGGCACTTGAACAGGTTGATTGAGTTGGCCTTGCCGCCCTGGCGCAGGTAATACCAAAGCCGCTAGGCCATCACACCCCGGACCTTGCCCAGGCTGCTCAGATCCGAGTCCGGGCGGTCATc                                                                                                                                               >  8:240663/1‑149 (MQ=255)
    tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGATCAGCTCCGGGTCCGGGCGATGACCGcc                                                                                                                                            <  2:182183/149‑1 (MQ=255)
      gcgGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                                                                                                                                           <  5:75677/148‑1 (MQ=255)
           ccAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCAcc                                                                                                                                     <  3:30209/149‑1 (MQ=255)
            cAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCAc                                                                                                                                      >  6:193018/1‑147 (MQ=255)
            cAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCAc                                                                                                                                      >  6:110024/1‑147 (MQ=255)
            cAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGGCCGGGCGATGACCGCCCggacccg                                                                                                                                    <  7:82870/149‑6 (MQ=255)
            cAGCCATTGGCTGGCCAGGCAGTTGAACAGGTGGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCa                                                                                                                                    <  3:235640/149‑1 (MQ=255)
              gCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAg                                                                                                                                   >  1:161415/1‑148 (MQ=255)
               ccATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGaa                                                                                                                                 >  4:367589/1‑149 (MQ=255)
                cATNTGCTGGCCANGCAGTTGTANAGGTTGATGGCGGTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAAt                                                                                                                                <  6:26328/149‑1 (MQ=255)
                         gcCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCAc                                                                                                                       >  6:162955/1‑149 (MQ=255)
                         gcCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCAc                                                                                                                       <  8:264087/149‑1 (MQ=255)
                            aGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGcc                                                                                                                    <  4:205689/149‑1 (MQ=255)
                              gCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAGAGCCGCTACGCCTTCTAAACGCGGACGTTGCCCAGGCTGGTAAGCTACGGGTACCGGCGGTGATAGCCCGGCAACAGAATAAGCTTCACGCCaa                                                                                                                  >  1:98809/1‑148 (MQ=255)
                                aGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACgg                                                                                                                <  3:367589/149‑1 (MQ=255)
                                aGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACgg                                                                                                                <  5:110024/149‑1 (MQ=255)
                                  ttGAACAGGTTGATGGCGTTGGCCTTGCCGCCCAGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGggc                                                                                                              <  8:39100/149‑1 (MQ=255)
                                    gAACCGGTATATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGggcg                                                                                                             <  1:319480/148‑1 (MQ=255)
                                    gAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTTATCGCCCGGCACCAGAATCAGCTGCACGCCACGggcg                                                                                                             <  7:310378/148‑1 (MQ=255)
                                    gAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGggcg                                                                                                             <  2:161415/148‑1 (MQ=255)
                                         ggTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGc                                                                                                       <  5:193018/149‑1 (MQ=255)
                                             gATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGa                                                                                                   >  3:246796/1‑149 (MQ=255)
                                                   gTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCAc                                                                                                                       >  8:37186/1‑123 (MQ=255)
                                                   gTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCAc                                                                                                                       <  7:37186/123‑1 (MQ=255)
                                                   gTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgtt                                                                                             <  6:355394/149‑1 (MQ=255)
                                                   gTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgtt                                                                                             >  7:49516/1‑149 (MQ=255)
                                                   gTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgtt                                                                                             <  4:307325/149‑1 (MQ=255)
                                                               gccCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGATGACCGCCCg                                                                                                                                          <  4:79625/92‑1 (MQ=255)
                                                                         aGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGCCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAg                                                                         >  7:3890/1‑147 (MQ=255)
                                                                           gTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGc                                                                     >  7:197723/1‑149 (MQ=255)
                                                                            tAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGcc                                                                    <  4:246796/149‑1 (MQ=255)
                                                                             aaTGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTG‑CCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCaa                                                                  >  7:30674/1‑149 (MQ=255)
                                                                                   aaaGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCa                                                             >  4:138813/1‑149 (MQ=255)
                                                                                         gCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCaa                                                                  >  4:338584/1‑138 (MQ=255)
                                                                                         gCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCaa                                                                  <  3:338584/138‑1 (MQ=255)
                                                                                         gCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTg                                                        >  7:293678/1‑148 (MQ=255)
                                                                                          cTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCa                                                      >  2:333689/1‑149 (MQ=255)
                                                                                               ccTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGac                                                 <  2:98809/149‑1 (MQ=255)
                                                                                                   cTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACAcc                                              <  8:197723/148‑1 (MQ=255)
                                                                                                    tCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAg                                            >  8:75874/1‑149 (MQ=255)
                                                                                                     cACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAg                                            >  2:371893/1‑148 (MQ=255)
                                                                                                     cACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAg                                            >  1:18992/1‑148 (MQ=255)
                                                                                                      cccGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGGTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGa                                          >  2:61260/2‑149 (MQ=255)
                                                                                                      aCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCCTCGACCGGCACAAGCATCAGCTGCACGCCAAGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGAGACACCAGGa                                          >  4:119671/1‑149 (MQ=255)
                                                                                                      aCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCTCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGa                                          <  1:106184/149‑1 (MQ=255)
                                                                                                        cgcgGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                         >  4:43455/1‑148 (MQ=255)
                                                                                                                ttGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggc                                <  4:61507/149‑1 (MQ=255)
                                                                                                                   cccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtg                             <  8:49516/149‑1 (MQ=255)
                                                                                                                   cccAGGCTGGTCAACTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtg                             >  7:165704/1‑149 (MQ=255)
                                                                                                                    ccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                            >  6:197923/1‑149 (MQ=255)
                                                                                                                    ccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                            <  8:30674/149‑1 (MQ=255)
                                                                                                                      aGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtggcg                          >  5:32894/1‑149 (MQ=255)
                                                                                                                        gCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCgca                        <  5:371083/149‑1 (MQ=255)
                                                                                                                             tCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCAcc                   >  7:305127/1‑149 (MQ=255)
                                                                                                                               aGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTg                 <  8:3890/149‑1 (MQ=255)
                                                                                                                                acTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTgg                <  8:165704/148‑1 (MQ=255)
                                                                                                                                 cTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGt               >  3:123439/1‑149 (MQ=255)
                                                                                                                                   ccgggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTc              >  5:55322/1‑148 (MQ=255)
                                                                                                                                    cgggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCg             <  1:333689/148‑1 (MQ=255)
                                                                                                                                     gggtccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGa            >  2:278075/1‑148 (MQ=255)
                                                                                                                                        gccgggCGATGACCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATa        >  8:82372/2‑149 (MQ=255)
                                                                                                                                              cGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTc  <  7:75874/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTC  >  NC_002947/3980136‑3980426

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