Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,491,034 Δ1 bp 9.7% coding (697/1902 nt) cysNC → sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,491,0340G. 9.7% 192.7 / 14.9 62coding (697/1902 nt)cysNCsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base . (2/4);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

CCTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGG  >  NC_002947/1490889‑1491182
                                                                                                                                                   |                                                                                                                                                     
ccTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGt                                                                                                                                                      <  3:240986/147‑1 (MQ=255)
  tGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAg                                                                                                                                                 >  8:2727/1‑149 (MQ=255)
    taGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGGCTTACGTTGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAg                                                                                                                                                 <  1:116479/148‑1 (MQ=255)
     ccttcatcGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGtt                                                                                                                                              <  8:8877/142‑1 (MQ=255)
      gTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGttt                                                                                                                                             <  3:131442/149‑1 (MQ=255)
         cGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCTCGACGGTTTCg                                                                                                                                           <  4:331777/148‑1 (MQ=255)
          gCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCTCGACGGTTTCga                                                                                                                                          <  2:14778/148‑2 (MQ=255)
           ccGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCt                                                                                                                                        <  2:336781/149‑1 (MQ=255)
             gAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGa                                                                                                                                      <  5:231118/149‑1 (MQ=255)
              aaGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGAc                                                                                                                                     >  3:142223/1‑149 (MQ=255)
                  cATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGc                                                                                                                                  >  6:242922/1‑148 (MQ=255)
                   aTCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCaa                                                                                                                                >  4:57643/1‑149 (MQ=255)
                    tCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAAc                                                                                                                               <  2:294147/149‑1 (MQ=255)
                        ccTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCa                                                                                                                           <  8:228447/149‑1 (MQ=255)
                              cccGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGAcc                                                                                                                     >  2:284309/1‑149 (MQ=255)
                                 gAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGc                                                                                                                  <  2:45795/149‑1 (MQ=255)
                                            ccTTCG‑‑TTCCGATGTCGGCGCTAAAGGGTGACACCGTGTTCAACCACTGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTg                                                                                                       <  8:145827/148‑1 (MQ=255)
                                            aCTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTg                                                                                                       <  7:2727/149‑1 (MQ=255)
                                              ttCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTg                                                                                                       >  5:126768/1‑147 (MQ=255)
                                                 g‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGt                                                                                                   <  4:142223/148‑1 (MQ=255)
                                                     gCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTAc                                                                                                 >  6:50296/1‑148 (MQ=255)
                                                     gCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTAc                                                                                                 <  5:50296/148‑1 (MQ=255)
                                                      ccGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGt                                                                                               <  5:242922/149‑1 (MQ=255)
                                                           gTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAAcc                                                                                          >  4:250937/1‑149 (MQ=255)
                                                             cGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGt                                                                                        >  6:63703/1‑149 (MQ=255)
                                                                    cAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAAcc                                                                                 >  4:129389/1‑149 (MQ=255)
                                                                     aaGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAAc                                                                                  <  5:63703/147‑1 (MQ=255)
                                                                       gggTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGAc                                                                                                                      >  7:273790/1‑109 (MQ=255)
                                                                       gggTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGAc                                                                                                                      <  8:273790/109‑1 (MQ=255)
                                                                       gggTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGa                                                                              <  5:262846/149‑1 (MQ=255)
                                                                          tGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCACGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACt                                                                           >  1:108368/1‑149 (MQ=255)
                                                                             cAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCTCGACGGTTTCga                                                                                                                                          >  4:4616/1‑82 (MQ=255)
                                                                             cAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCTCGACGGTTTCga                                                                                                                                          <  3:4616/83‑2 (MQ=255)
                                                                                cGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCcgcg                                                                     >  6:171043/1‑149 (MQ=255)
                                                                                          cacaGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCggg                                                           <  3:250937/149‑1 (MQ=255)
                                                                                          cacaGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCggg                                                           >  2:119666/1‑149 (MQ=255)
                                                                                           acaGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGc                                                          <  3:129389/149‑1 (MQ=255)
                                                                                              gCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTAc                                                                                                 >  2:202609/1‑107 (MQ=255)
                                                                                              gCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTAc                                                                                                 <  1:202609/107‑1 (MQ=255)
                                                                                                        gCCGTGGTACGCCGGCCCGACGCTGATGGAACTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCTACCGTCCCAACCTGACCTTCCGCGGCTTCGCGGGGACCATCGCCAGTg                                             <  5:308499/149‑1 (MQ=255)
                                                                                                        gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTg                                             >  8:341520/1‑149 (MQ=255)
                                                                                                        gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTg                                             <  7:341520/149‑1 (MQ=255)
                                                                                                              gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtg                                        <  5:171043/148‑1 (MQ=255)
                                                                                                              gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtg                                        <  6:126768/148‑1 (MQ=255)
                                                                                                               tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggt                                      <  2:108368/149‑1 (MQ=255)
                                                                                                                          gACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCggg                                                           <  4:90485/117‑1 (MQ=255)
                                                                                                                          gACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCggg                                                           >  3:90485/1‑117 (MQ=255)
                                                                                                                                 aTGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggtg                                     <  4:207665/132‑1 (MQ=255)
                                                                                                                                 aTGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggtg                                     >  3:207665/1‑132 (MQ=255)
                                                                                                                                   ggAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCa                                                >  3:4828/1‑119 (MQ=255)
                                                                                                                                   ggAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCa                                                <  4:4828/119‑1 (MQ=255)
                                                                                                                                   ggAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGAtt                   <  3:57643/148‑1 (MQ=255)
                                                                                                                                      aaTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCCGTGGCGTGGTGCACAAAGGTGACGAGATTGTCg               >  7:303042/1‑149 (MQ=255)
                                                                                                                                      aaTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCg               >  1:357529/1‑149 (MQ=255)
                                                                                                                                       aTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGTACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGt              <  6:145536/149‑1 (MQ=255)
                                                                                                                                        tCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtg             >  2:311053/1‑149 (MQ=255)
                                                                                                                                           tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctg          >  8:179136/1‑149 (MQ=255)
                                                                                                                                           ccGAAACCTCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctg          >  3:183733/2‑149 (MQ=255)
                                                                                                                                           ccGAAACCTCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgct           <  4:183733/147‑1 (MQ=255)
                                                                                                                                             gAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGcc        >  1:242710/1‑149 (MQ=255)
                                                                                                                                                  cGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg   <  3:309613/149‑1 (MQ=255)
                                                                                                                                                  cGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg   >  4:309613/1‑149 (MQ=255)
                                                                                                                                                   gTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCggg  >  6:129359/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGG  >  NC_002947/1490889‑1491182

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: