Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,874,651 T→G N257H (AAT→CAT)  cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,874,6510TG100.0% 318.1 / NA 84N257H (AAT→CAT) cysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (31/53);  total (31/53)

TGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGT  >  NC_000913/2874416‑2874885
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
gaccccgccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCgg                                                                                                                                                                                                                                         <  1:438827/235‑1 (MQ=255)
   aCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGc                                                                                                                                                                                                                                       <  1:885551/238‑1 (MQ=255)
   aCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGc                                                                                                                                                                                                                                       <  1:999193/238‑1 (MQ=255)
        cgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTa                                                                                                                                                                                                                                  <  1:1103901/238‑1 (MQ=255)
            aaCGCTTCGTCTGCCGGCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTa                                                                                                                                                                                                                             <  1:500599/239‑1 (MQ=255)
            aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGt                                                                                                                                                                                                                              >  2:876356/1‑238 (MQ=255)
            aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTa                                                                                                                                                                                                                             <  1:499115/239‑1 (MQ=255)
                       tGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGa                                                                                                                                                                                                                  <  1:1092793/239‑1 (MQ=255)
                       tGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGa                                                                                                                                                                                                                  <  1:498121/239‑1 (MQ=255)
                       tGCAGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGa                                                                                                                                                                                                                  <  2:1044027/239‑1 (MQ=255)
                       ggccgccAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGa                                                                                                                                                                                                                  <  2:604815/238‑1 (MQ=255)
                                       tCGCCACGTCTGATGTCGATCTCATCCGTCAGCCCCCGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCCCCATCAAAAGTCGCGATCCGCGCGACGTTTGATTCCCCCCCAGAGGGCAGCCCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGa                                                                                                                                                                                                    <  1:596357/237‑1 (MQ=255)
                                       tCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCt                                                                                                                                                                                                      <  1:651849/235‑1 (MQ=255)
                                         gCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCCCACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCAt                                                                                                                                                                                                <  1:893339/239‑1 (MQ=255)
                                          ccACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCAt                                                                                                                                                                                                >  1:26971/1‑238 (MQ=255)
                                             cGGCTGATGCCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATccac                                                                                                                                                                                            <  1:355018/239‑1 (MQ=255)
                                                cTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATccacc                                                                                                                                                                                           >  2:721640/1‑237 (MQ=255)
                                                    tGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGc                                                                                                                                                                                     <  1:702503/239‑1 (MQ=255)
                                                         aTCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTc                                                                                                                                                                                       >  1:11239/1‑232 (MQ=255)
                                                              aTCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCa                                                                                                                                                                           <  2:11239/239‑1 (MQ=255)
                                                               tCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTa                                                                                                                                                                                                                                  <  1:602480/183‑1 (MQ=255)
                                                               tCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTa                                                                                                                                                                                                                                  >  2:602480/1‑183 (MQ=255)
                                                                ccGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGACCTCCACg                                                                                                                                                                         >  1:1116352/1‑239 (MQ=255)
                                                                        aCCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTc                                                                                                                                                                    <  1:1062218/236‑1 (MQ=255)
                                                                         ccAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGc                                                                                                                                                                <  2:26971/239‑1 (MQ=255)
                                                                             ggggATCGCCTCTCCGGCAAAGGCTTCTTCGCGATCACCAACAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAATTCGAGATGCGGGCGATTAACGTACTGCACCGGGCAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACt                                                                                                                                                             >  2:731473/4‑238 (MQ=255)
                                                                             ggTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACt                                                                                                                                                             >  1:878904/1‑238 (MQ=255)
                                                                                 aTCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTg                                                                                                                                                                                   <  2:99397/212‑1 (MQ=255)
                                                                                 aTCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTg                                                                                                                                                                                   >  1:99397/1‑212 (MQ=255)
                                                                                 aTCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTc                                                                                                                                                                    >  2:342381/1‑227 (MQ=255)
                                                                                 aTCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTc                                                                                                                                                                    <  1:342381/227‑1 (MQ=255)
                                                                                  tCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAg                                                                                                                                                        <  1:393720/238‑1 (MQ=255)
                                                                                  tCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAg                                                                                                                                                        <  1:772870/238‑1 (MQ=255)
                                                                                   cGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGATGCCCAACTTCCACGCGACCGGATCCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCCGTTGATCTCCACGGTTTCCAGCACCTCGAGCa                                                                                                                                                      >  2:1074992/1‑239 (MQ=255)
                                                                                        ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAg                                                                                                                                                        <  2:667744/232‑1 (MQ=255)
                                                                                        ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAg                                                                                                                                                        >  1:667744/1‑232 (MQ=255)
                                                                                            ggCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGAc                                                                                                                                             <  2:878904/239‑1 (MQ=255)
                                                                                               aaaGGCTTCTTCGCGATCACCATCACAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCg                                                                                                                                           <  1:263510/238‑1 (MQ=255)
                                                                                               aaaGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCg                                                                                                                                           <  1:272925/238‑1 (MQ=255)
                                                                                                   gCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAg                                                                                                                                                        <  1:566136/221‑1 (MQ=255)
                                                                                                   gCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAg                                                                                                                                                        >  2:566136/1‑221 (MQ=255)
                                                                                                         tcGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACgg                                                                                                                                <  2:753551/239‑1 (MQ=255)
                                                                                                         tcGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACgg                                                                                                                                <  1:769681/239‑1 (MQ=255)
                                                                                                         tcGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACgg                                                                                                                                <  2:1069784/239‑1 (MQ=255)
                                                                                                           gcgATCCCCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACGCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCa                                                                                                                              <  2:447789/239‑1 (MQ=255)
                                                                                                                           gTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCGGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTTCGa                                                                                                              >  2:929991/1‑239 (MQ=255)
                                                                                                                             cACGATCCGCGCGACGTTTGATTCCACACCAGGGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGGACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATg                                                                                                            <  1:682173/239‑1 (MQ=255)
                                                                                                                              aCGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGc                                                                                                                                                                                                       <  1:18304/147‑1 (MQ=255)
                                                                                                                              aCGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGc                                                                                                                                                                                                       >  2:18304/1‑147 (MQ=255)
                                                                                                                                gATCCGCGCAGCGTTTGATTCCACACCAGATGTCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCa                                                                                                         <  2:639764/239‑1 (MQ=255)
                                                                                                                                     gcgcgACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCg                                                                                                                                           <  2:36272/200‑1 (MQ=255)
                                                                                                                                     gcgcgACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCg                                                                                                                                           >  1:36272/1‑200 (MQ=255)
                                                                                                                                      cgcgACGTTTGATTCCACACCAGAGGTCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGCCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCAc                                                                                                        >  1:639764/1‑234 (MQ=255)
                                                                                                                                                            gagGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCa                                                                                                                                                                           <  1:644530/145‑1 (MQ=255)
                                                                                                                                                            gagGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCa                                                                                                                                                                           >  2:644530/1‑145 (MQ=255)
                                                                                                                                                            gagGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCg                                                                              >  2:885316/1‑238 (MQ=255)
                                                                                                                                                               ggCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGa                                                                                                                                                                                                    <  1:1048242/117‑1 (MQ=255)
                                                                                                                                                               ggCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGa                                                                                                                                                                                                    >  2:1048242/1‑117 (MQ=255)
                                                                                                                                                                gcagcaCTTTGACACGTTGCCCGACTTGCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCgg                                                                             <  1:745482/235‑1 (MQ=255)
                                                                                                                                                                 cagcaCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTa                                                                                                                                                                                                                             >  2:267442/1‑90 (MQ=255)
                                                                                                                                                                 cagcaCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTa                                                                                                                                                                                                                             <  1:267442/90‑1 (MQ=255)
                                                                                                                                                                 cagcaCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAg                                                                                                                                                     >  1:271743/1‑162 (MQ=255)
                                                                                                                                                                 cagcaCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAg                                                                                                                                                     <  2:271743/162‑1 (MQ=255)
                                                                                                                                                                   gcaCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAg                                                                      <  2:1023733/239‑1 (MQ=255)
                                                                                                                                                                    caCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGc                                                                                                                                                                                                                                       >  2:318667/1‑77 (MQ=255)
                                                                                                                                                                    caCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGc                                                                                                                                                                                                                                       <  1:318667/77‑1 (MQ=255)
                                                                                                                                                                    caCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTa                                                                                                                                      <  1:649775/174‑1 (MQ=255)
                                                                                                                                                                    caCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTa                                                                                                                                      >  2:649775/1‑174 (MQ=255)
                                                                                                                                                                    caCTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGc                                                                     >  1:309359/1‑239 (MQ=255)
                                                                                                                                                                         tGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCg                                                                                                 >  2:176353/1‑206 (MQ=255)
                                                                                                                                                                         tGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCg                                                                                                 <  1:176353/206‑1 (MQ=255)
                                                                                                                                                                         tGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGata                                                                >  2:289045/1‑239 (MQ=255)
                                                                                                                                                                          gACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGatat                                                               <  1:819446/239‑1 (MQ=255)
                                                                                                                                                                                        cTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCgg                                                  <  1:329514/238‑1 (MQ=255)
                                                                                                                                                                                           ccACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCCCAAAGCGGATATCCAGATTACCCGGCAGc                                              >  2:1151406/1‑239 (MQ=255)
                                                                                                                                                                                              cgcgACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                                                >  1:108907/1‑234 (MQ=255)
                                                                                                                                                                                                     ggATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCaa                                    <  1:876356/239‑1 (MQ=255)
                                                                                                                                                                                                      gATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCg                                                                                                                                                <  1:569343/130‑1 (MQ=255)
                                                                                                                                                                                                      gATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCg                                                                                                                                                >  2:569343/1‑130 (MQ=255)
                                                                                                                                                                                                               gTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCaaa                            <  1:391704/237‑1 (MQ=255)
                                                                                                                                                                                                               gTTCCGGCGTAACCACGAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCaaa                            <  2:497010/237‑1 (MQ=255)
                                                                                                                                                                                                                               gAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGatat                                                               >  1:974/1‑186 (MQ=255)
                                                                                                                                                                                                                               gAAAATCGAGATGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGatat                                                               <  2:974/186‑1 (MQ=255)
                                                                                                                                                                                                                                        gaTGCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGt  <  1:929991/238‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
TGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCAAATAATCTTCACGAATACGGGTGAACGT  >  NC_000913/2874416‑2874885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: