Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,853,629 T→A 6.1% H491L (CAC→CTC)  dnaX ← DNA polymerase III subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,853,6290TA 6.1% 190.1 / 2.3 65H491L (CAC→CTC) dnaXDNA polymerase III subunit gamma
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base A (2/2);  total (29/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 = 6.35e-01

CCTGAACATGTTCGGCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTCGACCACTTCGGCGACCGGCTCAGGT  >  NC_002947/4853483‑4853775
                                                                                                                                                  |                                                                                                                                                  
ccTGAACATGTTCGGCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTTGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCgtgt                                                                                                                                                  >  6:81296/1‑149 (MQ=255)
  tGAACATGTTCGGCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGc                                                                                                                                                >  1:14892/1‑149 (MQ=255)
  tGAACATGTTCGGCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGc                                                                                                                                                >  3:193724/1‑149 (MQ=255)
            cGGCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGggcggc                                                                                                                                       <  4:193724/148‑1 (MQ=255)
             ggCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCt                                                                                                                                     <  2:81452/149‑1 (MQ=255)
              gCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCtg                                                                                                                                    >  7:4233/1‑149 (MQ=255)
               ccAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCtgt                                                                                                                                   >  6:106449/1‑149 (MQ=255)
               ccAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCtgt                                                                                                                                   <  8:4233/149‑1 (MQ=255)
               ccAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCtgt                                                                                                                                   >  6:146638/1‑149 (MQ=255)
               ccAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCtgt                                                                                                                                   <  2:39841/149‑1 (MQ=255)
                   ttCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCATGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTgg                                                                                                                               <  8:270767/149‑1 (MQ=255)
                      gTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGgc                                                                                                                             <  6:112289/148‑1 (MQ=255)
                             tAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCgg                                                                                                                     <  5:74625/149‑1 (MQ=255)
                                 tGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGc                                                                                                                  >  2:45019/1‑148 (MQ=255)
                                  gAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGcc                                                                                                                                      <  6:256313/127‑1 (MQ=255)
                                  gAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCaa                                                                                                                >  2:244229/1‑149 (MQ=255)
                                        tGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCt                                                                                                          >  3:159399/1‑149 (MQ=255)
                                         gCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTc                                                                                                         <  5:146638/149‑1 (MQ=255)
                                          ccGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCg                                                                                                        >  8:233017/1‑149 (MQ=255)
                                            ggATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCggg                                                                                                      <  7:73123/149‑1 (MQ=255)
                                               tCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGgcgc                                                                                                   <  1:244229/149‑1 (MQ=255)
                                                 gCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCa                                                                                                  <  8:103842/148‑1 (MQ=255)
                                                  cTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAgg                                                                                                >  1:47837/1‑149 (MQ=255)
                                                  cTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAgg                                                                                                >  5:125905/1‑149 (MQ=255)
                                                  cTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGAGTGCCGGGGCGGCCTGTGTGGGCGCCACCAGGGCAACAGGCTCGGGCGCAgg                                                                                                >  2:100748/1‑149 (MQ=255)
                                                   tttCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGt                                                                                               >  3:254509/1‑149 (MQ=255)
                                                        ccACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCgg                                                                                          >  3:111614/1‑149 (MQ=255)
                                                        ccACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCgg                                                                                          <  8:25273/149‑1 (MQ=255)
                                                           ccGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGtt                                                                                        >  5:175508/1‑148 (MQ=255)
                                                           ccGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTc                                                                                       <  7:233017/149‑1 (MQ=255)
                                                           ccGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTc                                                                                       <  3:137102/149‑1 (MQ=255)
                                                               aaTAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCTGGTTCAg                                                                                     >  2:85394/1‑147 (MQ=255)
                                                                   gTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAg                                                                               <  4:90466/149‑1 (MQ=255)
                                                                        cGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTc                                                                           <  1:307387/148‑1 (MQ=255)
                                                                        cGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGAGTGCCGGGGCGGCCTGTGTGGGCGCCACCAGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGtt                                                                            <  1:100748/147‑1 (MQ=255)
                                                                          tCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGc                                                                        >  6:126091/1‑149 (MQ=255)
                                                                          tCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGc                                                                        <  5:15087/149‑1 (MQ=255)
                                                                            aaGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGt                                                                                               <  2:295896/124‑1 (MQ=255)
                                                                            aaGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGt                                                                                               >  1:295896/1‑124 (MQ=255)
                                                                            aaGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGAAGCGTAGGCAGACGGGTCGAAAGGCGGCTAGTCGTGTGCCGGGGAGGCCTGTGTGGGCGCCACCGGGGCAACAGGCACGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCg                                                                      >  7:197124/1‑149 (MQ=255)
                                                                                 gcggcTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAAc                                                                 >  6:250893/1‑149 (MQ=255)
                                                                                        aTCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCgg                                                          <  1:166787/149‑1 (MQ=255)
                                                                                         tcgtcgCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCggg                                                         <  6:175508/149‑1 (MQ=255)
                                                                                           gtcgCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAACCGCGGCGCGTCGTGTGCCCGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGCGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCggggg                                                       <  4:111614/149‑1 (MQ=255)
                                                                                           gtcgCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCgggg                                                        <  5:250893/148‑1 (MQ=255)
                                                                                               cgcTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAAc                                                   >  2:210816/1‑149 (MQ=255)
                                                                                                 cTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACgg                                                 <  7:194260/149‑1 (MQ=255)
                                                                                                  tCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACgg                                                 >  1:280064/1‑148 (MQ=255)
                                                                                                    cATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCa                                               <  1:45019/148‑1 (MQ=255)
                                                                                                      tCCCCAGGTCAGCGTAGGCTGACTGGTCGAAAGGCTGCTCGTCGTGTGCCGGGGCGGCCTTTGTGGGCGCGAGCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCATGGGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGc                                            <  8:197124/149‑1 (MQ=255)
                                                                                                      tCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGc                                            <  5:81296/149‑1 (MQ=255)
                                                                                                         ccACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGctct                                         >  8:208992/1‑149 (MQ=255)
                                                                                                         ccACGGAAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGAGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGACGCTAACGGCGCGGGGGCAACGGCAGGctct                                         >  1:39490/1‑149 (MQ=255)
                                                                                                          cACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCtt                                        <  6:125905/149‑1 (MQ=255)
                                                                                                           aCGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTc                                       >  1:188079/1‑149 (MQ=255)
                                                                                                             ggCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCGGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTccc                                     >  6:145441/1‑149 (MQ=255)
                                                                                                               cAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCAc                                   <  2:80542/149‑1 (MQ=255)
                                                                                                               cAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCAc                                   <  5:126091/149‑1 (MQ=255)
                                                                                                                   gTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACgg                                 >  1:247952/1‑147 (MQ=255)
                                                                                                                      ggCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGt                            <  2:14892/149‑1 (MQ=255)
                                                                                                                      ggCAGACGGGTCGAAAGGCGGCTCGTCGAGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAgg                             <  2:39490/148‑1 (MQ=255)
                                                                                                                       gCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCCCCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTc                           >  8:52438/1‑149 (MQ=255)
                                                                                                                                  gAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTCGACCACTTCg                 <  4:159399/148‑1 (MQ=255)
                                                                                                                                      ggcggcTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTCGACCACTTCGGCGAc            <  2:280064/149‑1 (MQ=255)
                                                                                                                                      ggcggcTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTCGACCACTTCGGCGAc            <  4:254509/149‑1 (MQ=255)
                                                                                                                                         ggcTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTCGACCACTTCGGCGACCgg         >  4:159744/1‑149 (MQ=255)
                                                                                                                                                 gtgtGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTCGACCACTTCGGCGACCGGCTCAGGt  <  5:106449/148‑1 (MQ=255)
                                                                                                                                                  tgtgCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTccc                                     >  2:149133/1‑112 (MQ=255)
                                                                                                                                                  tgtgCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTccc                                     <  1:149133/112‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCTGAACATGTTCGGCCAATTCGTCCAGGTAGCTGAAACCTGCCGGATCGCTTTCGCCACCGTAATAGTCTTCGTCAAGCGGCGGCTCATCGTCGCGCTCCATCCCCACGGCAGCGTAGGCAGACGGGTCGAAAGGCGGCTCGTCGTGTGCCGGGGCGGCCTGTGTGGGCGCCACCGGGGCAACAGGCTCGGGCGCAGGTGCCGGTTCAGGCGCAGGTTCAGCCGCTAACGGCGCGGGGGCAACGGCAGGCTCTTCCCACGGCAGGTCGACCACTTCGGCGACCGGCTCAGGT  >  NC_002947/4853483‑4853775

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: