Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,951,835 A→C 5.3% I569S (ATC→AGC)  mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,951,8350AC 5.3% 275.0 / 4.6 89I569S (ATC→AGC) mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (42/41);  new base C (3/2);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTGG  >  NC_002947/1951687‑1951978
                                                                                                                                                    |                                                                                                                                               
gCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGa                                                                                                                                                 <  7:135116/149‑1 (MQ=255)
  cACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATa                                                                                                                                               <  7:315550/149‑1 (MQ=255)
   aCTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAg                                                                                                                                              <  5:340928/149‑1 (MQ=255)
   aCTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCACTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAg                                                                                                                                              >  5:352077/1‑149 (MQ=255)
    cTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAgg                                                                                                                                             >  4:113570/1‑149 (MQ=255)
               cGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCCCTATCGCcgacgcgca                                                                                                                                  <  3:7445/149‑5 (MQ=255)
                 tGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCgg                                                                                                                                 >  6:407652/1‑148 (MQ=255)
                 tGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCgg                                                                                                                                 <  5:407652/148‑1 (MQ=255)
                   aTCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCa                                                                                                                               <  4:56005/148‑1 (MQ=255)
                   aTCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCAAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCGGGCCCAGTTCTGCATAGGCTTCGGCAAACGCCTGCACGTCGGCGATAGGGCCGACGATCGGCa                                                                                                                               <  1:104735/148‑1 (MQ=255)
                    tCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCcgac                                                                                                                                       <  4:297254/139‑1 (MQ=255)
                    tCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCcgac                                                                                                                                       >  3:297254/1‑139 (MQ=255)
                    tCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAgg                                                                                                                             <  1:137475/149‑1 (MQ=255)
                      aaCCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCTAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGTTCGGCAGGt                                                                                                                            >  7:16506/1‑148 (MQ=255)
                      aaCCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGt                                                                                                                            >  4:46557/1‑148 (MQ=255)
                        cccGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAAt                                                                                                                         >  6:312560/1‑149 (MQ=255)
                         ccGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATc                                                                                                                        >  7:237252/1‑149 (MQ=255)
                              tGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGc                                                                                                                   <  6:214502/149‑1 (MQ=255)
                                 gCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCa                                                                                                                               >  1:394605/1‑134 (MQ=255)
                                 gCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCa                                                                                                                               <  2:394605/134‑1 (MQ=255)
                                  ccATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGt                                                                                                               <  5:199632/149‑1 (MQ=255)
                                  ccATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGt                                                                                                               >  4:269185/1‑149 (MQ=255)
                                        ccGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCGACg                                                                                                         >  3:14301/1‑149 (MQ=255)
                                                cacaTCCACCCCGCCCACCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACgca                                                                                                  <  8:286240/148‑1 (MQ=255)
                                                 acaTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCTAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGTTCGGCAGGTAATCCGGGCTGGTGCAACGGAACgcagc                                                                                                <  8:16506/149‑1 (MQ=255)
                                                  caTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGc                                                                                                                                >  4:131258/1‑116 (MQ=255)
                                                  caTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGc                                                                                                                                <  3:131258/116‑1 (MQ=255)
                                                  caTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACgcagc                                                                                                <  3:113570/148‑1 (MQ=255)
                                                   aTACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCAc                                                                                              >  5:330520/1‑149 (MQ=255)
                                                    tACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGAGAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACg                                                                                             >  4:325796/1‑149 (MQ=255)
                                                    atcAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACg                                                                                             <  5:414318/147‑1 (MQ=255)
                                                       aaCCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGt                                                                                                               >  8:154041/1‑128 (MQ=255)
                                                       aaCCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGt                                                                                                               <  7:154041/128‑1 (MQ=255)
                                                       aaCCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTcc                                                                                          >  6:166640/1‑149 (MQ=255)
                                                       aaCCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTcc                                                                                          <  2:124861/149‑1 (MQ=255)
                                                         cccGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCt                                                                                        <  8:270476/149‑1 (MQ=255)
                                                          ccGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTg                                                                                       >  3:100430/1‑149 (MQ=255)
                                                            gCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGc                                                                                      >  6:145230/1‑148 (MQ=255)
                                                                  ccAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                              <  2:297741/102‑1 (MQ=255)
                                                                  ccAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAg                                                                                                                              >  1:297741/1‑102 (MQ=255)
                                                                    aatggcagggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCTTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGTCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCg                                                                             <  6:175596/146‑1 (MQ=255)
                                                                    aaGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCa                                                                                                                               <  4:349151/99‑1 (MQ=255)
                                                                    aaGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCa                                                                                                                               >  3:349151/1‑99 (MQ=255)
                                                                    aaGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCg                                                                             >  2:302604/1‑149 (MQ=255)
                                                                    aaGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCg                                                                             <  6:352077/149‑1 (MQ=255)
                                                                    aaGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCg                                                                             <  3:46557/149‑1 (MQ=255)
                                                                       ggcaggtcaAATCGGTTAATTGCCGCGCATCCCGTCCCATTACTGCATAGGATTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCggg                                                                           <  6:271950/148‑1 (MQ=255)
                                                                        gcagggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCAGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCCCGTCCCTGCAATTGCTCGGGGt                                                                         >  4:364398/1‑149 (MQ=255)
                                                                           gggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCa                                                                      <  8:237252/149‑1 (MQ=255)
                                                                            ggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCGACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAg                                                                     <  4:14301/149‑1 (MQ=255)
                                                                            ggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAAtt                                                                                  <  7:373741/136‑1 (MQ=255)
                                                                            ggcaCATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAg                                                                     >  4:326759/1‑149 (MQ=255)
                                                                               acaTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGtt                                                                   >  4:343539/1‑148 (MQ=255)
                                                                                caTCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTca                                                                 <  4:124351/149‑1 (MQ=255)
                                                                                  tCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTcagc                                                               >  7:267592/1‑149 (MQ=255)
                                                                                     ggTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTcagccac                                                            <  7:291639/149‑1 (MQ=255)
                                                                                             ccGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACagccagcc                                                    >  1:170921/1‑149 (MQ=255)
                                                                                             ccGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACagccagcc                                                    <  1:302604/149‑1 (MQ=255)
                                                                                                   aTCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTgcg                                              >  5:174596/1‑149 (MQ=255)
                                                                                                     cccGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTc                                                                              <  7:325697/115‑1 (MQ=255)
                                                                                                     cccGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTc                                                                              >  8:325697/1‑115 (MQ=255)
                                                                                                       cGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTgcgc                                             >  7:375061/1‑146 (MQ=255)
                                                                                                       cGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTgcgc                                             >  2:391246/1‑146 (MQ=255)
                                                                                                           ccAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATc                                                                                                                                   <  8:270456/56‑1 (MQ=255)
                                                                                                           ccAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATc                                                                                                                                   >  7:270456/1‑56 (MQ=255)
                                                                                                           ccAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTgcgc                                             <  8:330113/142‑1 (MQ=255)
                                                                                                           ccAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTgcgc                                             >  7:330113/1‑142 (MQ=255)
                                                                                                           ccAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATc                                       >  2:125464/1‑148 (MQ=255)
                                                                                                            cAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGa                                     >  1:235007/1‑149 (MQ=255)
                                                                                                             aGTACGGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGAc                                    <  3:325796/149‑1 (MQ=255)
                                                                                                               tACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTGCAGCCACAGCCAGCCGTTGCGCCAGATCGACgg                                  <  6:186642/149‑1 (MQ=255)
                                                                                                                    cATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGaaagt                             >  3:299776/1‑147 (MQ=255)
                                                                                                                    cATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAAtt                             >  1:354472/1‑149 (MQ=255)
                                                                                                                         gCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTc                        <  5:312560/149‑1 (MQ=255)
                                                                                                                         gCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTc                        <  4:100430/149‑1 (MQ=255)
                                                                                                                          ctgctgCAAACGCATGCGCTTCGGCGTTCGGGCCGACGATCGGCAGTTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTcc                       <  4:299776/144‑1 (MQ=255)
                                                                                                                               gCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCcagcag                  >  5:345957/1‑149 (MQ=255)
                                                                                                                                cAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCagcagc                 >  6:337815/1‑149 (MQ=255)
                                                                                                                                cAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCCGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCagcagc                 >  2:233178/1‑149 (MQ=255)
                                                                                                                                 aaaCGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCg                                                   <  3:69846/114‑1 (MQ=255)
                                                                                                                                 aaaCGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCg                                                   >  4:69846/1‑114 (MQ=255)
                                                                                                                                  aaCGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAgcgc               <  8:267592/149‑1 (MQ=255)
                                                                                                                                     tgccgatcgTCGGCCCTATCGCCGACGATCGGCAGGTAATCCGGGCTTGTGCAACGGAACGCACCACGTCCCTGCAATTGCTCGGGGTCCAGGTCAGCCCCACCCAGCCGTTGCGCCATATCGACGGATATTTAGTCCAGCAGCGCCAg            >  4:280164/8‑149 (MQ=255)
                                                                                                                                        cgatcgTCGGCCCTATCGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTca                                                                 >  3:326328/5‑93 (MQ=38)
                                                                                                                                        cgatcgTCGGCCCTATCGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTca                                                                 <  4:326328/89‑1 (MQ=38)
                                                                                                                                          atcgTCGGCCCTATCGCCGACGATCGGAAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTAAGCCCCAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGCTTGc       >  7:267949/3‑149 (MQ=255)
                                                                                                                                           gcgTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGcc      <  1:391246/149‑1 (MQ=255)
                                                                                                                                            cgTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGccc     >  3:169917/1‑149 (MQ=255)
                                                                                                                                              tCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTg   <  5:337815/149‑1 (MQ=255)
                                                                                                                                               cGGCGATAGGGCCGACGATCGGCAGGTAATCCTGGCTGGTGCAACGGAACGCAGCAAGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTTCGCAAGATCGACGGAAATTTCTTCCAGCAGCGCCAGGTTGCCCTgc  >  8:377366/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCATCCCGGCCCAGTACTGCATAGGCTTCGGCAAACGCCTGCGCGTCGGCGATAGGGCCGACGATCGGCAGGTAATCCGGGCTGGTGCAACGGAACGCAGCACGTCCCTGCAATTGCTCGGGGTCCAGTTCAGCCACAGCCAGCCGTTGCGCCAGATCGACGGAAATTTCGTCCAGCAGCGCCAGGTTGCCCTGG  >  NC_002947/1951687‑1951978

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