Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,874,306 G→A 6.6% R372C (CGC→TGC)  cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,874,3060GA 6.6% 185.9 / 2.9 60R372C (CGC→TGC) cysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base A (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAG  >  NC_000913/2874165‑2874453
                                                                                                                                             |                                                                                                                                                   
tCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCa                                                                                                                                              <  2:116659/149‑1 (MQ=255)
 cAGCCCACCCGTCACCGGATTTTGTTGCTAACGATCTTACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGGGGCAGGTTTTCAACTTCACGATGGGTAAGGTTATTAATATCAACCTGATAGCGATTGCCAt                                                                                                                                             <  2:788838/149‑1 (MQ=255)
  aGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATc                                                                                                                                            <  1:1093893/149‑1 (MQ=255)
  aGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATc                                                                                                                                            <  1:1340631/149‑1 (MQ=255)
  aGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGACATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATc                                                                                                                                            <  2:1223074/149‑1 (MQ=255)
       aCCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAacac                                                                                                                                       >  2:1598868/1‑149 (MQ=255)
       aCCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAacac                                                                                                                                       >  2:519635/1‑149 (MQ=255)
           gTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACAcgcgc                                                                                                                                   >  1:54912/1‑149 (MQ=255)
                 ggATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCTCGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGt                                                                                                                             <  2:195798/149‑1 (MQ=255)
                 ggATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGCTCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCACCCTGATAGCGAATGCCATCACCCCGCGCGCGCGt                                                                                                                             >  2:287114/1‑149 (MQ=255)
                  gATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGtc                                                                                                                            >  2:125040/1‑149 (MQ=255)
                  gATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGtc                                                                                                                            <  2:1433139/149‑1 (MQ=255)
                        gttgATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTa                                                                                                                      <  2:149307/149‑1 (MQ=255)
                         ttgATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTAc                                                                                                                     <  2:590735/149‑1 (MQ=255)
                          tgATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTAcc                                                                                                                    <  1:1486866/149‑1 (MQ=255)
                             tAACGATCTACCACCAGCGGCTCGTCAACAGTGAGATCCTCGAGGACGATCCCATTCAGTGGCAGGTTTTCAACTTCCCGCTGGGTAAGGTTATTAATATCACCCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGc                                                                                                                 >  1:242612/1‑149 (MQ=255)
                                   tCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAAtttt                                                                                                           <  2:1315384/149‑1 (MQ=255)
                                     tAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGa                                                                                                         <  1:1320353/149‑1 (MQ=255)
                                     tAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGa                                                                                                         <  1:1371795/149‑1 (MQ=255)
                                      aaCACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGACGATCCCATTCAGTGGCAGGTTTTCAACTTAAAGCTGCGTAAGGTTATAAATATAAAACTGATAGCGAATGCCATCAACACGCGCGCGCATCTTCTTACCGGCaaaattcat                                                                                                        >  1:1003456/1‑142 (MQ=255)
                                           cAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGt                                                                                                   >  2:545601/1‑149 (MQ=255)
                                                 cTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGACAATTTTGATGTCGTAActct                                                                                             <  1:1037499/149‑1 (MQ=255)
                                                    gTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGcc                                                                                          <  2:1286737/149‑1 (MQ=255)
                                                    gTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGcc                                                                                          >  1:136196/1‑149 (MQ=255)
                                                    gTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCAAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGcc                                                                                          >  1:1391815/1‑149 (MQ=255)
                                                    gTCAAAAGTGAGATCCAAGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGcc                                                                                          >  1:218654/1‑149 (MQ=255)
                                                       aaaaGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTg                                                                                       >  2:1542971/1‑149 (MQ=255)
                                                         aaGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGga                                                                                     <  2:1473547/149‑1 (MQ=255)
                                                         aaGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGga                                                                                     <  1:545601/149‑1 (MQ=255)
                                                         aaGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGga                                                                                     >  2:751133/1‑149 (MQ=255)
                                                            tGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTAGCGGCATTTTTGATGTCGTACCTCTGCCCTGGAGaa                                                                                  >  1:275087/1‑149 (MQ=255)
                                                                          ccGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCAAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGc                                                                    <  2:1391815/149‑1 (MQ=255)
                                                                             aTCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAATGTTATTAATATAAAACTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTATTTCCGGCAATTTTGATGTCGGAACTCTGCCCTGGAGAACGCGGCGGTTCCCccaa                                                                 >  2:791257/1‑148 (MQ=255)
                                                                                  aatCAGTGGCAGGTTATCAAATTCACGCTGGGAAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACAGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATa                                                            <  1:553907/147‑1 (MQ=255)
                                                                                  aTTCAGTGGCAGGTTTTCAACTTCACTCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGGCTACTTACCGGCTATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCTGCTGTTCCGCCATCCATa                                                            >  2:983244/1‑149 (MQ=255)
                                                                                  aTTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATa                                                            >  1:1097655/1‑149 (MQ=255)
                                                                                  aTTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCAAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATa                                                            >  2:1357848/1‑149 (MQ=255)
                                                                                   ttCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATAc                                                           >  2:1085447/1‑149 (MQ=255)
                                                                                   ttCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATAc                                                           <  2:653076/149‑1 (MQ=255)
                                                                                   ttCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATAc                                                           <  2:511137/149‑1 (MQ=255)
                                                                                     cAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACca                                                         <  1:1357343/149‑1 (MQ=255)
                                                                                           cAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATATTGATGTCGTAAATCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCa                                                   >  2:1424397/1‑149 (MQ=255)
                                                                                                tttCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGAc                                              <  1:519635/149‑1 (MQ=255)
                                                                                                tttCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGAc                                              <  1:1598868/149‑1 (MQ=255)
                                                                                                         cACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGc                                     <  1:1539710/149‑1 (MQ=255)
                                                                                                          aCGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCa                                    >  1:1172773/1‑149 (MQ=255)
                                                                                                                    aGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTaa                          >  2:1151385/1‑149 (MQ=255)
                                                                                                                    aGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTaa                          >  2:234562/1‑149 (MQ=255)
                                                                                                                     ggTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAAc                         <  1:1085447/149‑1 (MQ=255)
                                                                                                                      gTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACg                        >  2:942175/1‑149 (MQ=255)
                                                                                                                       ttattaATATCAACCTGATAGCGACTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGc                       <  1:287114/149‑1 (MQ=255)
                                                                                                                        tattaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCt                      <  1:125040/149‑1 (MQ=255)
                                                                                                                          ttaATATCAACCTGATAGCAAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTc                    <  1:1357848/149‑1 (MQ=255)
                                                                                                                              tatCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCt                <  1:608541/149‑1 (MQ=255)
                                                                                                                               atCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTg               >  1:1246258/1‑149 (MQ=255)
                                                                                                                                 cAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTgcc             >  1:820166/1‑149 (MQ=255)
                                                                                                                                   aCCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTgccgc           >  1:1435644/1‑149 (MQ=255)
                                                                                                                                         tAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCcagcag     <  1:751133/149‑1 (MQ=255)
                                                                                                                                            cGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCCCCGCCGGTAACGCTTCGTCTGCCGCcagcagcag  >  1:236057/1‑149 (MQ=255)
                                                                                                                                            cGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCcagcagcag  >  2:1524912/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAG  >  NC_000913/2874165‑2874453

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