Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 895,210 G→C 5.5% intergenic (+50/+14) PP_0777 → / ← PP_0778 glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947895,2100GC 5.5% 235.0 / 3.1 72intergenic (+50/+14)PP_0777/PP_0778glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base C (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

GTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAG  >  NC_002947/895065‑895358
                                                                                                                                                  |                                                                                                                                                     
gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCgcgggcg                                                                                                                                                    >  1:377082/1‑149 (MQ=255)
gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCgcgg                                                                                                                                                       <  3:214060/146‑1 (MQ=255)
    gCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGgcgc                                                                                                                                               >  5:397949/1‑149 (MQ=255)
     cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                              >  5:407537/1‑149 (MQ=255)
     cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                              <  8:267509/149‑1 (MQ=255)
      aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTg                                                                                                                                             >  2:415865/1‑149 (MQ=255)
      aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTg                                                                                                                                             <  2:324468/149‑1 (MQ=255)
       aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTgg                                                                                                                                            <  6:397949/149‑1 (MQ=255)
        ggATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGa                                                                                                                                           <  1:379238/149‑1 (MQ=255)
          aTGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAtc                                                                                                                                         <  1:333730/149‑1 (MQ=255)
          aTGGCCGGGTGCTGGCGCGTTTCGCGCCGCGCAAGGCGCCGGATGACCCAGCGGTTGTGCAGGCCATCGACAAACCCCTGGCCTGAATAGTGCCTTAGCGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAt                                                                                                                                          <  1:108475/149‑1 (MQ=255)
               cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                      >  4:292969/1‑147 (MQ=255)
               cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTcaca                                                                                                                                    <  1:317442/149‑1 (MQ=38)
                gggTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGTGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAg                                                                                                                                   <  2:274502/149‑1 (MQ=255)
                  gTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGgc                                                                                                                                 <  3:18833/149‑1 (MQ=255)
                  gTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGgc                                                                                                                                 <  7:24476/149‑1 (MQ=255)
                    gCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGgcgc                                                                                                                               <  3:23604/149‑1 (MQ=255)
                      tGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTg                                                                                                                             <  8:362784/149‑1 (MQ=255)
                        gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTg                                                                                                                             <  6:16281/147‑1 (MQ=255)
                        gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                            >  4:361979/1‑148 (MQ=255)
                        gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTca                                                                                                                                      <  2:123817/138‑1 (MQ=38)
                        gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTca                                                                                                                                      >  1:123817/1‑138 (MQ=38)
                          gcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCaaa                                                                                                                         <  2:377082/149‑1 (MQ=255)
                           cgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCaaa                                                                                                                         >  8:312229/1‑148 (MQ=255)
                           cgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAAc                                                                                                                        <  6:407537/149‑1 (MQ=255)
                            gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTgg                                                                                                                                            >  5:100675/1‑128 (MQ=255)
                            gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTgg                                                                                                                                            <  6:100675/128‑1 (MQ=255)
                                   gccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCCGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                                    <  1:40157/146‑1 (MQ=255)
                                    ccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                                  >  7:201622/1‑146 (MQ=255)
                                      gcgcACGGCGCCGTTTGACCCAGCGGTGGTGCAGGCAAGCCTGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGTATCTCAAAGGCGCTGCAAACATTGCGCAGaa                                                                                                             <  8:201622/149‑1 (MQ=255)
                                         cACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGc                                                                                                          <  1:112634/149‑1 (MQ=255)
                                              cgcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGc                                                                                                     >  2:17191/1‑149 (MQ=255)
                                                 ctgATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGt                                                                                                  <  7:312229/147‑1 (MQ=255)
                                                 cGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCg                                                                                                   >  5:46810/1‑148 (MQ=255)
                                                   gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTCGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAg                                                                                                <  5:379922/149‑1 (MQ=255)
                                                         ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                          <  6:2196/149‑1 (MQ=255)
                                                           aGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGt                                                                                        >  1:303171/1‑149 (MQ=255)
                                                                tggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcca                                                                                   <  2:303171/149‑2 (MQ=255)
                                                                tggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCg                                                                                   <  1:17191/149‑1 (MQ=255)
                                                                     cAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                                    >  3:335588/1‑111 (MQ=255)
                                                                     cAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                                    <  4:335588/111‑1 (MQ=255)
                                                                     cAGGCCATAGAGAAAGCCCTGGCCTGAAGACTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATAGCGCAGAACGCCCAGCCGTAGCACCAGGATGCCggactt                                                                              >  5:47093/1‑145 (MQ=255)
                                                                              gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                       >  5:276148/1‑147 (MQ=255)
                                                                              gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCa                                                                     <  6:46810/149‑1 (MQ=255)
                                                                               agaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAg                                                                    <  8:93829/149‑1 (MQ=255)
                                                                                 aaaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGa                                                                   >  7:7250/1‑148 (MQ=255)
                                                                                  aaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGAt                                                                  <  3:361979/148‑1 (MQ=255)
                                                                                       cTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCt                                                            >  2:422050/1‑149 (MQ=255)
                                                                                            cTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAAc                                                       >  2:289790/1‑149 (MQ=255)
                                                                                            cTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAAc                                                       <  3:292969/149‑1 (MQ=255)
                                                                                               aaGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGt                                                    >  4:109833/1‑149 (MQ=255)
                                                                                                   gTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGa                                                >  1:21901/1‑149 (MQ=255)
                                                                                                              tGATGTTTTGGCAACTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCgg                                      >  1:161870/1‑149 (MQ=255)
                                                                                                                  gTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAg                                  >  2:299361/1‑149 (MQ=255)
                                                                                                                  gTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAg                                  >  4:169905/1‑149 (MQ=255)
                                                                                                                   ttttGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGc                                 >  8:6290/1‑149 (MQ=255)
                                                                                                                        gCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg                            >  3:123582/1‑149 (MQ=255)
                                                                                                                        gCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg                            >  3:411702/1‑149 (MQ=255)
                                                                                                                        gCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg                            >  4:324308/1‑149 (MQ=255)
                                                                                                                        gCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg                            >  7:312430/1‑149 (MQ=255)
                                                                                                                         cACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggc                           >  3:406616/1‑149 (MQ=255)
                                                                                                                         cACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggc                           <  4:406616/149‑1 (MQ=255)
                                                                                                                          aCCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcg                          >  3:27787/1‑149 (MQ=255)
                                                                                                                          aCCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcg                          >  3:74369/1‑149 (MQ=255)
                                                                                                                            cTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcggc                        >  7:197937/1‑149 (MQ=255)
                                                                                                                                 agaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCa                   <  5:198835/149‑1 (MQ=255)
                                                                                                                                  gaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAt                  <  2:259700/149‑1 (MQ=255)
                                                                                                                                  gaTCCAGCGCCGCCCGCGCGGCACTCGATCTCACAGGTGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAt                  >  7:151554/1‑149 (MQ=38)
                                                                                                                                         tGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGACATGCTGTcg           <  2:161870/149‑1 (MQ=255)
                                                                                                                                           ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGc                    <  2:23510/138‑1 (MQ=255)
                                                                                                                                           ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGc                    >  1:23510/1‑138 (MQ=255)
                                                                                                                                            cgcgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgcg         >  8:336801/1‑148 (MQ=255)
                                                                                                                                             gcgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCa       >  2:70740/1‑149 (MQ=255)
                                                                                                                                             ccgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGCCATCCCCTCTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCa       >  2:181943/2‑149 (MQ=255)
                                                                                                                                                 gggcggCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCa   <  8:197937/149‑1 (MQ=255)
                                                                                                                                                  ggcggcGCTGGATCTGACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAg  <  2:318323/149‑1 (MQ=255)
                                                                                                                                                  ggcggcGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAg  <  8:312430/149‑1 (MQ=255)
                                                                                                                                                  ggcggcGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAg  <  3:169905/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAG  >  NC_002947/895065‑895358

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: