Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,491,038:1 +C 9.8% coding (701/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,0381.C 9.8% 219.5 / 19.6 61coding (701/1902 nt)cysNCsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base . (25/30);  new base C (2/4);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG  >  NC_002947/1490891‑1491186
                                                                                                                                                      |                                                                                                                                                    
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                                                               >  2:119666/1‑149 (MQ=255)
                                                                                        cacaGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCggg                                                               <  3:250937/149‑1 (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                                                 <  7:341520/149‑1 (MQ=255)
                                                                                                      gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTg                                                 >  8:341520/1‑149 (MQ=255)
                                                                                                            gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtg                                            <  5:171043/148‑1 (MQ=255)
                                                                                                            gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtg                                            <  6:126768/148‑1 (MQ=255)
                                                                                                             tACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggt                                          <  2:108368/149‑1 (MQ=255)
                                                                                                                        gACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCggg                                                               >  3:90485/1‑117 (MQ=255)
                                                                                                                        gACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCggg                                                               <  4:90485/117‑1 (MQ=255)
                                                                                                                               aTGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggtg                                         >  3:207665/1‑132 (MQ=255)
                                                                                                                               aTGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggtg                                         <  4:207665/132‑1 (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)
                                                                                                                                                    gAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAACGTGACGAGATTGTCGTGCTGCCGTCGGGCaa   <  2:357529/149‑1 (MQ=255)
                                                                                                                                                     aGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCATCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAag  <  8:303042/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCG‑‑TGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAG  >  NC_002947/1490891‑1491186

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: