Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,508,565 A→G 5.2% G72G (GGT→GGC cysG ← multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,508,5650AG 5.2% 223.9 / 2.1 77G72G (GGT→GGCcysGmultifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGG  >  NC_002947/4508419‑4508710
                                                                                                                                                    |                                                                                                                                                 
ggcgAACGGTCCACGACCGCCGGGAACATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACCCCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAAcc                                                                                                                                                 <  1:114227/149‑1 (MQ=255)
   gAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACGACGCCGCCAACCgc                                                                                                                                               >  1:152013/1‑148 (MQ=255)
   gAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGAc                                                                                                                                              <  4:144981/149‑1 (MQ=255)
                tCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGcc                                                                                                                                 <  6:60360/149‑1 (MQ=255)
                  gCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCt                                                                                                                                >  5:171500/1‑148 (MQ=255)
                  gCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAAGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGcc                                                                                                                                 <  5:106145/147‑1 (MQ=255)
                   ccGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                              <  2:406055/149‑1 (MQ=255)
                   ccGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                              >  1:86273/1‑149 (MQ=255)
                   ccGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGCTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCCGGTCGGCTGCCTGa                                                                                                                              >  7:74695/1‑149 (MQ=255)
                    cGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGa                                                                                                                              >  3:77151/1‑148 (MQ=255)
                     gggAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAGCCGACCTGGTCGGTTGccggctg                                                                                                                            <  2:157773/149‑6 (MQ=255)
                     gggAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAGCCGACCTGGTCGGTTGccggctg                                                                                                                            <  7:48281/149‑6 (MQ=255)
                            tCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACg                                                                                                                      <  4:77151/148‑1 (MQ=255)
                                cGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCa                                                                                                                 <  6:228522/149‑1 (MQ=255)
                                    cacaGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCAc                                                                                                             <  6:171500/149‑1 (MQ=255)
                                     acaGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCAc                                                                                                             >  4:284279/1‑148 (MQ=255)
                                         ggCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCAc                                                                                                        <  7:148065/149‑1 (MQ=255)
                                          gCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCAcc                                                                                                       <  8:393395/149‑1 (MQ=255)
                                           ccGGCGCATCTACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCCCACCg                                                                                                      <  4:366461/149‑1 (MQ=255)
                                           ccGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAg                                                                                                                <  7:374743/139‑1 (MQ=255)
                                           ccGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAg                                                                                                                >  8:374743/1‑139 (MQ=255)
                                            cGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGc                                                                                                     <  2:86273/149‑1 (MQ=255)
                                             ggCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGcc                                                                                                    <  7:399535/149‑1 (MQ=255)
                                                  aTCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCtt                                                                                               >  7:407743/1‑149 (MQ=255)
                                                      accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCg                                                                                           >  1:114848/1‑149 (MQ=255)
                                                       ccacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCgg                                                                                          >  8:160839/1‑149 (MQ=255)
                                                       ccacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCgg                                                                                          >  4:349191/1‑149 (MQ=255)
                                                        cacaTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTAAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGCCAGATAGCCACGCACCCGCACCTCACCGCCACCTTACCggg                                                                                         >  5:274205/1‑149 (MQ=255)
                                                           aTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAg                                                                                                                <  8:270498/123‑1 (MQ=255)
                                                           aTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAg                                                                                                                >  7:270498/1‑123 (MQ=255)
                                                            ttGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCCGCACCTCACCGCCACCTTCCCGGGCCaa                                                                                     >  2:420226/1‑149 (MQ=255)
                                                               aCCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAAcgc                                                                                  >  3:198977/1‑149 (MQ=255)
                                                               aCCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAAcg                                                                                   <  5:341096/148‑1 (MQ=255)
                                                                 cGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCCCCAGCACCTCACCGCCACCTTCCCGGGCCACCGCgg                                                                                >  3:42551/1‑149 (MQ=255)
                                                                  ggCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGc                                                                               <  2:318106/149‑1 (MQ=255)
                                                                     aGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCGCCAGCCCCTCACCGCCACCTTCCCGGGGCAAcagatcg                                                                              >  4:135957/1‑141 (MQ=255)
                                                                     aGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCaa                                                                                     <  3:135957/140‑1 (MQ=255)
                                                                       gCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACt                                                                          >  5:56216/1‑149 (MQ=255)
                                                                         tAAGGGCCTGCGCATCGGCCGACACCTGGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCgg                                                                                >  1:260697/1‑142 (MQ=255)
                                                                         tAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCgg                                                                                <  2:260697/141‑1 (MQ=255)
                                                                         tAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTgg                                                                        >  2:247007/1‑149 (MQ=255)
                                                                         tAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTgg                                                                        >  3:125952/1‑149 (MQ=255)
                                                                              gCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGACGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTc                                                                   >  2:312103/1‑149 (MQ=255)
                                                                                   cgcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACg                                                               <  6:361264/148‑1 (MQ=255)
                                                                                     cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTcc                                                            >  3:28012/1‑149 (MQ=255)
                                                                                      aTCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCg                                                           >  1:344461/1‑149 (MQ=255)
                                                                                      aTCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCg                                                           >  5:143813/1‑149 (MQ=255)
                                                                                            cGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGcca                                                     <  4:28012/149‑1 (MQ=255)
                                                                                              acacCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGccac                                                    >  4:386747/1‑148 (MQ=255)
                                                                                                acCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCcaccac                                                 >  2:339180/1‑149 (MQ=255)
                                                                                                   tGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCAcgcg                                              >  8:57474/1‑149 (MQ=255)
                                                                                                     gggCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCAcgcg                                              >  6:169626/1‑148 (MQ=255)
                                                                                                      ggCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACgcgc                                             <  6:143813/148‑1 (MQ=255)
                                                                                                      ggCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgca                                            <  3:284279/149‑1 (MQ=255)
                                                                                                        cATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCCCCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgcag                                           >  1:312269/1‑148 (MQ=255)
                                                                                                        cATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgcagc                                          <  2:114848/149‑1 (MQ=255)
                                                                                                              gCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTgg                                                                        >  1:419043/1‑113 (MQ=255)
                                                                                                              gCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTgg                                                                        <  2:419043/113‑1 (MQ=255)
                                                                                                               cccAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCAcgccgc                                   >  4:64246/1‑149 (MQ=255)
                                                                                                               cccAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCAcgccgc                                   >  8:246834/1‑149 (MQ=255)
                                                                                                                 cAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTGCGGCGCCACCACGCGCAGCACGccgccg                                 <  1:242692/149‑1 (MQ=255)
                                                                                                                 cAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGccgccg                                 <  2:152013/149‑1 (MQ=255)
                                                                                                                   gggTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGc                               <  5:169626/149‑1 (MQ=255)
                                                                                                                     gtcgtcAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCAt                             >  4:412392/1‑149 (MQ=255)
                                                                                                                        gtcAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCg                           <  8:407743/148‑1 (MQ=255)
                                                                                                                        gtcAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCg                           <  7:359286/148‑1 (MQ=255)
                                                                                                                                   cgATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgc               <  2:344461/149‑1 (MQ=255)
                                                                                                                                    gATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACACCGCCGGCATCGGCCAGCAGgcgcg              >  6:417336/1‑149 (MQ=255)
                                                                                                                                      tCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGcc            <  1:247007/149‑1 (MQ=255)
                                                                                                                                        aCCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCt           <  1:339180/148‑1 (MQ=255)
                                                                                                                                         ccAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTg         >  6:27961/1‑149 (MQ=255)
                                                                                                                                              ccGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCa     <  3:64246/149‑1 (MQ=255)
                                                                                                                                                 gCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAgg  >  6:307203/1‑149 (MQ=255)
                                                                                                                                                 gCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCAACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGTGCAAgg  >  6:310890/1‑149 (MQ=255)
                                                                                                                                                 gCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGccacca                                                  <  8:200081/101‑1 (MQ=255)
                                                                                                                                                 gCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGccacca                                                  >  7:200081/1‑101 (MQ=255)
                                                                                                                                                 gCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAgg  <  3:50870/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGAACGGTCCACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGG  >  NC_002947/4508419‑4508710

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: