Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,874,505 T→C 15.6% G305G (GGA→GGG)  cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,874,5050TC15.6% 130.4 / 3.4 66G305G (GGA→GGG) cysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (40/11);  new base C (4/7);  total (45/21)
Fisher's exact test for biased strand distribution p-value = 9.73e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCAGCACTT  >  NC_000913/2874287‑2874731
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                      
ttaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGCAGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCCCCA‑CAGCCAGATCGCCCCGGCCGATGTCGATCCCATCCGCCAGCACCAGGGGGAGCG‑C‑TTCTCCGGCaaaa                                                                                                                                                                                                                               >  2:127976/1‑227 (MQ=255)
ttaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGGGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGCCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TCCTCCCGCaacg                                                                                                                                                                                                                               >  1:325025/1‑226 (MQ=255)
      tCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCTCCGGTAACGCTTCGTCTGCCGCCCGCAC‑CCGATCGCCACGGCAGATGTCGACCCCATCCGTCAGCCCCACGGCG‑TCGCC‑CTCTCCGGCAAAGGcttct                                                                                                                                                                                                                         >  2:647102/1‑228 (MQ=255)
       cAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAC‑CCGATCCCCCCGGCCGCGGTCGATCCCATCCCTCCGCCCCCGGGTGACCG‑C‑TTCACCGGCCAAGGcttcat                                                                                                                                                                                                                        >  2:558570/1‑226 (MQ=255)
       cAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACCCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCCCCCCGGCTGATGTCGATCTCATCCGTCAGCACCAGGGCGCACG‑C‑TTTCCCCGCAAAGGgttcct                                                                                                                                                                                                                        >  2:436081/1‑226 (MQ=255)
          ccTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCATAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCCGCACCGCCGGTAACCCTTCGCCTGCCGCCAGCAG‑CCGACCGCCACGACGGACGTCGATCCCATCCGTCACCACCAGGGACACCG‑C‑TTCTCACGCAAAGGccacctcaa                                                                                                                                                                                                                     >  2:62187/1‑220 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCCCGGCTGATGTCGCTCTCATCCGTCAGCCCCAGCGGGATCG‑C‑TTCTCCGGCAAAGGCTCCTTCCCGATc                                                                                                                                                                                                                 >  2:619178/1‑228 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCTTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCttcttcccccta                                                                                                                                                                                                                 >  2:644413/1‑222 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGAt                                                                                                                                                                                                                  >  1:359783/1‑227 (MQ=255)
                    aaTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCCTAACTCTGCCCTGCAGAAAGCGGCTCTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTCACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCCCGGCTGACTCCGATCTCACCTGCCCGCCCCCTGCTCCTCC‑C‑TTATCCGCCAAAGCCCTCTTCCCTACCACCAt                                                                                                                                                                                                            >  2:387136/1‑227 (MQ=255)
                    aaTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGGCGGAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGGAGGCCTCCAGCAG‑CACATCGCAACGGCCGATGCAGATGTCATCCGTCAGCCCCACAG‑GATAG‑CGTTCTCCGGCACAGGCCGCTCCGCGTACAccca                                                                                                                                                                                                            >  2:448002/1‑225 (MQ=255)
                     aTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CCGATCGCCACGGCTGATGTCAACCTCTTCCGCCAGCACCACGGTGACCC‑C‑TCCTCCGGCAAAGGCTCTTTCCCGCTCACCaac                                                                                                                                                                                                           >  2:655966/1‑225 (MQ=255)
                            cAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTcc                                                                                                                                                                                                   >  1:370084/1‑227 (MQ=255)
                                     cgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGCCGATCTCATCCGTCAGCACCAGGGCGATCA‑C‑TTCTCCGGCAAAGGCTTCTTCCCGATCACCCTCCAAAGTCACGATCcgc                                                                                                                                                                                           >  1:89041/1‑227 (MQ=255)
                                     cgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCACAAG‑CAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCCCCAGGGAGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGACACGATCcgc                                                                                                                                                                                           >  1:655191/1‑227 (MQ=255)
                                       cgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACCGCTGATGTCGATCTCATCCATCAGCACCAGGGAGATCG‑C‑CTCCCCGGCAAAGGCTCCTTCGCGATCACCATCCAAAGTCACGATCcccgg                                                                                                                                                                                         >  2:627940/1‑226 (MQ=255)
                                       cgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAACTCGCCACGGCCGATGTCGGTCTCCTCCGTCAGCCCCAGGGTGATCG‑C‑TTCCCCGGCCAAGGCCTCTTCCCCATCACCATCAAAAGCCACCATCcgcag                                                                                                                                                                                         >  2:504295/1‑225 (MQ=255)
                                                      gCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑Ccac                                                                                                                                                                        >  1:13350/1‑227 (MQ=255)
                                                                           tGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCT‑C‑TTCTCCGGCAAAGGCTTCTTCGCGTTCACCATCAAAAGTCACGATCCCCGCG‑ACGTTTGATT‑CCCCACCAAAGGGCAGAACttgaa                                                                                                                                                    >  2:316848/1‑223 (MQ=255)
                                                                            gCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCCCGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCCTCTTCGCCATCACCATCCAAAGTCACCATCCCCGCG‑ACCTTTGATT‑CACAACCAGAGGGCACACCTTTGac                                                                                                                                                   >  2:380095/1‑226 (MQ=255)
                                                                            gCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGATCACCATCAAAAGTCACGATCCGCCCG‑ACGTTTGATT‑CCACACCCGCGGGCAGCACTTTGaca                                                                                                                                                  >  1:354343/1‑227 (MQ=255)
                                                                             cccTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGCTCCCCATCAAAAGCCACAGTCCCCCCG‑ACGTTTGATT‑CCCAACCCAAGGgccgtgctcttaa                                                                                                                                                   >  2:21329/1‑214 (MQ=255)
                                                                                  gagaAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑Ttct                                                                                                                                                                                                                                        >  1:83474/1‑137 (MQ=255)
                                                                                  gagaAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGCCAGCACTTTGACACGTTGc                                                                                                                                            >  1:246434/1‑227 (MQ=255)
                                                                                    gaAAGCGGCTCTCCCGCCACCCCTCCCACCTCCCCCGCCGCGCTCGCCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGcc                                                                                                                                           <  1:267590/226‑1 (MQ=255)
                                                                                      agcgggCTGTTCCGCCTCCAATCCCCCACCCACCGCCGCGCTCTGCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGAc                                                                                                                                       <  1:619178/224‑1 (MQ=255)
                                                                                       aGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCCCCATCAAAAGTCACGATCCCCGCG‑CCGTTTGCTT‑CCCCACCACAGGGCCGCACTTTGACCCGTTGCCCCACt                                                                                                                                      >  2:541391/1‑228 (MQ=255)
                                                                                       aGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGCTCACCATCAAAAGTCACGATCCGCCCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACCTTGACACAGTGCCCGACt                                                                                                                                      >  1:134029/1‑228 (MQ=255)
                                                                                       aGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCCCACCAGAGGGCAGCACTTTGAAAAGTGGCCCgcc                                                                                                                                       >  2:267590/1‑225 (MQ=255)
                                                                                            cTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGCTGTCGATCTCCTCCGTCAGCCCCAGGGTGACCG‑C‑TTCTCCGGAAAAGGCTTCTTCGCGATCACCACCAAAAGTCACCACCCGCCCG‑ACGTTTGATT‑CCCCACCGAAGGGCAGCACTTTGACCCGTTGCCCCACTTTCaa                                                                                                                                 >  2:649396/1‑227 (MQ=255)
                                                                                            cTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCCCG‑CCGTTTGTTT‑CCACACCCGAGGGCAGCCCTTTGACACGTTGCCCGACTTCCa                                                                                                                                  >  1:378264/1‑227 (MQ=255)
                                                                                                                 aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGTCAAAGGCTTCTTCGCGATCACCATCACAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCATAGGTCAGCACTTTGACACGTTGCCCGCCTTCCACGCGCCCGGTGGCCAGCGTTcc                                                                                                            >  1:184213/1‑228 (MQ=255)
                                                                                                                 aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGATCCCCGACTTCCACGCGACCGGATTCCACCGTTcc                                                                                                            >  1:347477/1‑228 (MQ=255)
                                                                                                                 aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCCGAGGGAAGCACTTTCACACGTTGCCCGACTTCCACGCGACCGGTAGCCGACGTTcc                                                                                                            >  2:309460/1‑228 (MQ=255)
                                                                                                                 aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCACCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGATCACCATCAAAAGTCACGATCCGCCCG‑ACGTTTGATTCCCACCCCA‑AGGGCAGCACTTTTACACGTTTCCCGCGTTCCACCCACCCGGAGTCCAGCGTTcc                                                                                                            >  2:596749/1‑228 (MQ=255)
                                                                                                                        gACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCACG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCCACCGGATGCCCGCGTTCCGGccgaa                                                                                                     >  2:526843/1‑224 (MQ=255)
                                                                                                                           gcgcTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCT‑CCGTTTGATT‑CCACACCCGAGGGCCGCACTTTGCCCCGTTGCCCGCCTTCCCCGCGACCGATGGCCACCGTTCCGCCGTAcccc                                                                                                  >  1:198785/1‑227 (MQ=255)
                                                                                                                             gttcggccccgccgGTACCGCTTCGGCTGCCGCCAGCAG‑CAGTTCGCCAAGGCTGATTACGATCCCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACg                                                                                                <  1:316848/220‑1 (MQ=255)
                                                                                                                                        ccgGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGCAAATCGAGAt                                                                                     >  1:88515/1‑228 (MQ=255)
                                                                                                                                         cgGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCAACCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCCCGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGAAGCACTTTGACACGTTGCCCGACTTCCACCAGACCGGGTGCCAGCGTTCCGGGGCAACACCGAAAAtgcaaat                                                                                     >  2:433251/1‑221 (MQ=255)
                                                                                                                                            tAACGCTTCGTCTGCCGGCAGGAG‑CAGATCGCCCCGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGAACG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCgg                                                                                 <  2:13350/228‑1 (MQ=255)
                                                                                                                                                 cTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGGCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGAt                                                                            <  1:436081/228‑1 (MQ=255)
                                                                                                                                                   tccgcTGCCTCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTa                                                                          <  1:433251/224‑1 (MQ=255)
                                                                                                                                                   tCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGGGATTCGGGCGAtt                                                                           >  1:199347/1‑227 (MQ=255)
                                                                                                                                                        tGCCGCCAGCCG‑CAGATCGCCACGGCTGATGGCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  1:127976/228‑1 (MQ=255)
                                                                                                                                                        tGCCGCCAGCAG‑CAGATGGCCCCGGCTGATGCCGATCTCATCCGTCAGCTCCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  1:309460/228‑1 (MQ=255)
                                                                                                                                                         gccggcAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCCCCGGCAAAGGCTTCTGCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  1:380095/227‑1 (MQ=255)
                                                                                                                                                         gccgccAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  1:373452/227‑1 (MQ=255)
                                                                                                                                                             ccAGCAG‑CAGACCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCa                                                                <  1:196289/228‑1 (MQ=255)
                                                                                                                                                                       gATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGGGATTAACCGAATGCCCCGGGAAgc                                                       >  2:638406/1‑228 (MQ=255)
                                                                                                                                                                          cGCCACGGCTGATGTCGATTGCATCCGTCAGCACCAGGGTGATCG‑C‑TTCGCCGGCAAAGGCTTCCTCGCGATCACCATCAAAAGTCACGATCCGCTCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑Ca                                                   <  1:596749/228‑1 (MQ=255)
                                                                                                                                                                                        acggatgagatcGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTCACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACTGAAATCGACATTCGGGCGATTAAGGTCCTTGACC‑GGAAGAGCCATTGGCTGAGCCT‑cc                                    >  2:374057/12‑228 (MQ=255)
                                                                                                                                                                                            tctcCTCCCTCCGCACCAGGGTGATCG‑G‑TTCCTCGGCAAAGGCTTTCTCGCGATCACCAACAAAAGTCACCATCCG‑GCGCACTTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑ccacca                                <  2:359783/228‑1 (MQ=255)
                                                                                                                                                                                            tctcACCCGTCAGCACCAGGGTGAGCG‑C‑TTCTCCGGCAAAGGCTGCTTCGCGATCACCAGCAAAAGTCACGAACCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑ccacca                                <  2:89041/228‑1 (MQ=255)
                                                                                                                                                                                             cccacccGTCAACACCATGGTTATCG‑C‑TCTACCGGAAAAGGCTTTTCCGCTATCACGATCAAAGGTCACGATCCGCGCG‑CTGTTTGATC‑CCACACCAGAGGGCCGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑ccacca                                <  2:347477/222‑1 (MQ=255)
                                                                                                                                                                                                           ccAGGGTTATCG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCC‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCa                 <  1:644413/228‑1 (MQ=255)
                                                                                                                                                                                                            cAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAACTCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCa                 >  1:122286/1‑227 (MQ=255)
                                                                                                                                                                                                             aGGGTGATCG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACg               <  1:37861/228‑1 (MQ=255)
                                                                                                                                                                                                             aGGGTGATCG‑C‑TTCCCCCGCAAAGGCTTCTTCCCGATCCCCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACCCCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACg               <  1:638406/228‑1 (MQ=255)
                                                                                                                                                                                                              gggTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCTGGCGATTAACGTACTGCACCGGGAAGGG‑CATTGTCCGAGCAT‑CCACCACTCGCGGGAACCCCACgg              >  2:89142/1‑228 (MQ=255)
                                                                                                                                                                                                                   aTCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCTGGCGATTAACGTACTGCACCGGGAATCG‑CATTGGCTGTGCAT‑CCCCCACTCGCTGGATATCCACGGTTTcc         >  1:21444/1‑228 (MQ=255)
                                                                                                                                                                                                                   aTCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTc          >  1:481791/1‑227 (MQ=255)
                                                                                                                                                                                                                   aTCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGAAATTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACTACCTGCACCCGGAACCG‑CATTGGCACA‑CATACCACCACTCGCTGGTtccacccgtgttc          >  2:92298/1‑216 (MQ=255)
                                                                                                                                                                                                                    tCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCCCCGCGACGCG‑CATTGGCTGGGCAT‑CCCCCACTCCCTGACTACCCCCGGTTTccc        >  2:46976/1‑227 (MQ=255)
                                                                                                                                                                                                                     tccc‑TTCCCCGGCAAAGGTTTCTTCGCGATCCCCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGCCACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCa        <  1:243018/225‑1 (MQ=255)
                                                                                                                                                                                                                     cG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCCCCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCa        <  1:541391/227‑1 (MQ=255)
                                                                                                                                                                                                                     cG‑C‑TTCACCGGCAAAGGCTTCTGCGCGATCACCATCAAAAGCCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCa                 <  1:46976/218‑1 (MQ=255)
                                                                                                                                                                                                                            ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑cc                                                                                                                                                                          >  2:247456/1‑62 (MQ=255)
                                                                                                                                                                                                                            ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGCGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAACTCGCGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CAACCACTCGCTGGATCTCCCCGGGTTcaaaacctc  >  2:551590/1‑220 (MQ=255)
                                                                                                                                                                                                                            ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGCCTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGCTTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CAACCACTCGCTGGATCTCCACGGTTTcacgaactt  >  1:164457/1‑220 (MQ=255)
                                                                                                                                                                                                                             tcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑Cca                                                                                                                                                                         <  1:247456/62‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                      
TTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCAGCACTT  >  NC_000913/2874287‑2874731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑