Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,874,503 T→C 17.4% 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,5030TC17.4% 131.8 / 4.4 69G305G (GGA→GGG) cysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (40/12);  new base C (5/8);  total (46/23)
Fisher's exact test for biased strand distribution p-value = 1.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCAGCACTT  >  NC_000913/2874285‑2874729
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
ttaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGCAGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCCCCA‑CAGCCAGATCGCCCCGGCCGATGTCGATCCCATCCGCCAGCACCAGGGGGAGCG‑C‑TTCTCCGGCaaaa                                                                                                                                                                                                                               >  1:127976/1‑227 (MQ=255)
ttaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGGGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGCCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TCCTCCCGCaacg                                                                                                                                                                                                                               >  2:325025/1‑226 (MQ=255)
      tCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCTCCGGTAACGCTTCGTCTGCCGCCCGCAC‑CCGATCGCCACGGCAGATGTCGACCCCATCCGTCAGCCCCACGGCG‑TCGCC‑CTCTCCGGCAAAGGcttct                                                                                                                                                                                                                         >  1:647102/1‑228 (MQ=255)
       cAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAC‑CCGATCCCCCCGGCCGCGGTCGATCCCATCCCTCCGCCCCCGGGTGACCG‑C‑TTCACCGGCCAAGGcttcat                                                                                                                                                                                                                        >  1:558570/1‑226 (MQ=255)
       cAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACCCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCCCCCCGGCTGATGTCGATCTCATCCGTCAGCACCAGGGCGCACG‑C‑TTTCCCCGCAAAGGgttcct                                                                                                                                                                                                                        >  1:436081/1‑226 (MQ=255)
          ccTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCATAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCCGCACCGCCGGTAACCCTTCGCCTGCCGCCAGCAG‑CCGACCGCCACGACGGACGTCGATCCCATCCGTCACCACCAGGGACACCG‑C‑TTCTCACGCAAAGGccacctcaa                                                                                                                                                                                                                     >  1:62187/1‑220 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCCCGGCTGATGTCGCTCTCATCCGTCAGCCCCAGCGGGATCG‑C‑TTCTCCGGCAAAGGCTCCTTCCCGATc                                                                                                                                                                                                                 >  1:619178/1‑228 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCTTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCttcttcccccta                                                                                                                                                                                                                 >  1:644413/1‑222 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGAt                                                                                                                                                                                                                  >  2:359783/1‑227 (MQ=255)
                    aaTGCCATCCACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCACCCATCCCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTCCCGCCACCCG‑CAGAACGCCCCCGCTGAGATCGAGCACATCCATCACCACCCGGTGTATCC‑C‑TTCCCCGGCaaatgacttttaccggcccccccc                                                                                                                                                                                                           >  1:538626/1‑207 (MQ=255)
                    aaTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCCTAACTCTGCCCTGCAGAAAGCGG‑CTCTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTCACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCCCGGCTGACTCCGATCTCACCTGCCCGCCCCCTGCTCCTCC‑C‑TTATCCGCCAAAGCCCTCTTCCCTACCACCAt                                                                                                                                                                                                            >  1:387136/1‑227 (MQ=255)
                    aaTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGGCGGAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGGAGGCCTCCAGCAG‑CACATCGCAACGGCCGATGCAGATGTCATCCGTCAGCCCCACAG‑GATAG‑CGTTCTCCGGCACAGGCCGCTCCGCGTACAccca                                                                                                                                                                                                            >  1:448002/1‑225 (MQ=255)
                     aTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CCGATCGCCACGGCTGATGTCAACCTCTTCCGCCAGCACCACGGTGACCC‑C‑TCCTCCGGCAAAGGCTCTTTCCCGCTCACCaac                                                                                                                                                                                                           >  1:655966/1‑225 (MQ=255)
                            cAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTcc                                                                                                                                                                                                   >  2:370084/1‑227 (MQ=255)
                                     cgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGCCGATCTCATCCGTCAGCACCAGGGCGATCA‑C‑TTCTCCGGCAAAGGCTTCTTCCCGATCACCCTCCAAAGTCACGATCcgc                                                                                                                                                                                           >  2:89041/1‑227 (MQ=255)
                                     cgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCACAAG‑CAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCCCCAGGGAGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGACACGATCcgc                                                                                                                                                                                           >  2:655191/1‑227 (MQ=255)
                                       cgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACCGCTGATGTCGATCTCATCCATCAGCACCAGGGAGATCG‑C‑CTCCCCGGCAAAGGCTCCTTCGCGATCACCATCCAAAGTCACGATCcccgg                                                                                                                                                                                         >  1:627940/1‑226 (MQ=255)
                                       cgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAACTCGCCACGGCCGATGTCGGTCTCCTCCGTCAGCCCCAGGGTGATCG‑C‑TTCCCCGGCCAAGGCCTCTTCCCCATCACCATCAAAAGCCACCATCcgcag                                                                                                                                                                                         >  1:504295/1‑225 (MQ=255)
                                                     taataggagccagccgaaaCTCTGCCC‑GGGGAAAACGGTCTCTCCCCCCACCCATCCCACACCCACGGCCGCGCTCGGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑Cca                                                                                                                                                                         <  2:21329/210‑1 (MQ=255)
                                                      gCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑Ccac                                                                                                                                                                        >  2:13350/1‑227 (MQ=255)
                                                                           tGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCT‑C‑TTCTCCGGCAAAGGCTTCTTCGCGTTCACCATCAAAAGTCACGATCCCCGCG‑ACGTTTGATT‑CCCCACCAAAGGGCAGAACttgaa                                                                                                                                                    >  1:316848/1‑223 (MQ=255)
                                                                            gCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCCCGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCCTCTTCGCCATCACCATCCAAAGTCACCATCCCCGCG‑ACCTTTGATT‑CACAACCAGAGGGCACACCTTTGac                                                                                                                                                   >  1:380095/1‑226 (MQ=255)
                                                                            gCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGATCACCATCAAAAGTCACGATCCGCCCG‑ACGTTTGATT‑CCACACCCGCGGGCAGCACTTTGaca                                                                                                                                                  >  2:354343/1‑227 (MQ=255)
                                                                             cccTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGCTCCCCATCAAAAGCCACAGTCCCCCCG‑ACGTTTGATT‑CCCAACCCAAGGgccgtgctcttaa                                                                                                                                                   >  1:21329/1‑214 (MQ=255)
                                                                                  gagaAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑Ttct                                                                                                                                                                                                                                        >  2:83474/1‑137 (MQ=255)
                                                                                  gagaAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGCCAGCACTTTGACACGTTGc                                                                                                                                            >  2:246434/1‑227 (MQ=255)
                                                                                    gaAAGCGG‑CTCTCCCGCCACCCCTCCCACCTCCCCCGCCGCGCTCGCCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGcc                                                                                                                                           <  2:267590/226‑1 (MQ=255)
                                                                                      agcggg‑CTGTTCCGCCTCCAATCCCCCACCCACCGCCGCGCTCTGCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGAc                                                                                                                                       <  2:619178/224‑1 (MQ=255)
                                                                                       aGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCCCCATCAAAAGTCACGATCCCCGCG‑CCGTTTGCTT‑CCCCACCACAGGGCCGCACTTTGACCCGTTGCCCCACt                                                                                                                                      >  1:541391/1‑228 (MQ=255)
                                                                                       aGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGCTCACCATCAAAAGTCACGATCCGCCCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACCTTGACACAGTGCCCGACt                                                                                                                                      >  2:134029/1‑228 (MQ=255)
                                                                                       aGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCCCACCAGAGGGCAGCACTTTGAAAAGTGGCCCgcc                                                                                                                                       >  1:267590/1‑225 (MQ=255)
                                                                                             cTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGCTGTCGATCTCCTCCGTCAGCCCCAGGGTGACCG‑C‑TTCTCCGGAAAAGGCTTCTTCGCGATCACCACCAAAAGTCACCACCCGCCCG‑ACGTTTGATT‑CCCCACCGAAGGGCAGCACTTTGACCCGTTGCCCCACTTTCaa                                                                                                                                 >  1:649396/1‑227 (MQ=255)
                                                                                             cTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCCCG‑CCGTTTGTTT‑CCACACCCGAGGGCAGCCCTTTGACACGTTGCCCGACTTCCa                                                                                                                                  >  2:378264/1‑227 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGTCAAAGGCTTCTTCGCGATCACCATCACAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCATAGGTCAGCACTTTGACACGTTGCCCGCCTTCCACGCGCCCGGTGGCCAGCGTTcc                                                                                                            >  2:184213/1‑228 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGATCCCCGACTTCCACGCGACCGGATTCCACCGTTcc                                                                                                            >  2:347477/1‑228 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCCGAGGGAAGCACTTTCACACGTTGCCCGACTTCCACGCGACCGGTAGCCGACGTTcc                                                                                                            >  1:309460/1‑228 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCACCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCCCGATCACCATCAAAAGTCACGATCCGCCCG‑ACGTTTGATTCCCACCCCA‑AGGGCAGCACTTTTACACGTTTCCCGCGTTCCACCCACCCGGAGTCCAGCGTTcc                                                                                                            >  1:596749/1‑228 (MQ=255)
                                                                                                                         gACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCACG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCCACCGGATGCCCGCGTTCCGGccgaa                                                                                                     >  1:526843/1‑224 (MQ=255)
                                                                                                                            gcgcTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCT‑CCGTTTGATT‑CCACACCCGAGGGCCGCACTTTGCCCCGTTGCCCGCCTTCCCCGCGACCGATGGCCACCGTTCCGCCGTAcccc                                                                                                  >  2:198785/1‑227 (MQ=255)
                                                                                                                              gttcggccccgccgGTACCGCTTCGGCTGCCGCCAGCAG‑CAGTTCGCCAAGGCTGATTACGATCCCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACg                                                                                                <  2:316848/220‑1 (MQ=255)
                                                                                                                                         ccgGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGCAAATCGAGAt                                                                                     >  2:88515/1‑228 (MQ=255)
                                                                                                                                          cgGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCAACCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCCCGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGAAGCACTTTGACACGTTGCCCGACTTCCACCAGACCGGGTGCCAGCGTTCCGGGGCAACACCGAAAAtgcaaat                                                                                     >  1:433251/1‑221 (MQ=255)
                                                                                                                                             tAACGCTTCGTCTGCCGGCAGGAG‑CAGATCGCCCCGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGAACG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCgg                                                                                 <  1:13350/228‑1 (MQ=255)
                                                                                                                                                  cTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGGCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGAt                                                                            <  2:436081/228‑1 (MQ=255)
                                                                                                                                                    tccgcTGCCTCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTa                                                                          <  2:433251/224‑1 (MQ=255)
                                                                                                                                                    tCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGGGATTCGGGCGAtt                                                                           >  2:199347/1‑227 (MQ=255)
                                                                                                                                                         tGCCGCCAGCCG‑CAGATCGCCACGGCTGATGGCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:127976/228‑1 (MQ=255)
                                                                                                                                                         tGCCGCCAGCAG‑CAGATGGCCCCGGCTGATGCCGATCTCATCCGTCAGCTCCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:309460/228‑1 (MQ=255)
                                                                                                                                                          gccggcAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCCCCGGCAAAGGCTTCTGCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:380095/227‑1 (MQ=255)
                                                                                                                                                          gccgccAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:373452/227‑1 (MQ=255)
                                                                                                                                                              ccAGCAG‑CAGACCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCa                                                                <  2:196289/228‑1 (MQ=255)
                                                                                                                                                                        gATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGGGATTAACCGAATGCCCCGGGAAgc                                                       >  1:638406/1‑228 (MQ=255)
                                                                                                                                                                           cGCCACGGCTGATGTCGATTGCATCCGTCAGCACCAGGGTGATCG‑C‑TTCGCCGGCAAAGGCTTCCTCGCGATCACCATCAAAAGTCACGATCCGCTCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑Ca                                                   <  2:596749/228‑1 (MQ=255)
                                                                                                                                                                                    ttattcgcattcccgccGGCAGCCCCAGGGTGATCG‑C‑TTCCCCGGCAAAGGCTTCCTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑GCGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGa                                           <  1:88515/212‑1 (MQ=255)
                                                                                                                                                                                         acggatgagatcGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTCACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACTGAAATCGACATTCGGGCGATTAAGGTCCTTGACC‑GGAAGAGCCATTGGCTGAGCCT‑cc                                    >  1:374057/12‑228 (MQ=255)
                                                                                                                                                                                             tctcCTCCCTCCGCACCAGGGTGATCG‑G‑TTCCTCGGCAAAGGCTTTCTCGCGATCACCAACAAAAGTCACCATCCG‑GCGCACTTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑ccacca                                <  1:359783/228‑1 (MQ=255)
                                                                                                                                                                                             tctcACCCGTCAGCACCAGGGTGAGCG‑C‑TTCTCCGGCAAAGGCTGCTTCGCGATCACCAGCAAAAGTCACGAACCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑ccacca                                <  1:89041/228‑1 (MQ=255)
                                                                                                                                                                                              cccacccGTCAACACCATGGTTATCG‑C‑TCTACCGGAAAAGGCTTTTCCGCTATCACGATCAAAGGTCACGATCCGCGCG‑CTGTTTGATC‑CCACACCAGAGGGCCGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑ccacca                                <  1:347477/222‑1 (MQ=255)
                                                                                                                                                                                                            ccAGGGTTATCG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCC‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCa                 <  2:644413/228‑1 (MQ=255)
                                                                                                                                                                                                             cAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAACTCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCa                 >  2:122286/1‑227 (MQ=255)
                                                                                                                                                                                                              aGGGTGATCG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACg               <  2:37861/228‑1 (MQ=255)
                                                                                                                                                                                                              aGGGTGATCG‑C‑TTCCCCCGCAAAGGCTTCTTCCCGATCCCCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACCCCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACg               <  2:638406/228‑1 (MQ=255)
                                                                                                                                                                                                               gggTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCTGGCGATTAACGTACTGCACCGGGAAGGG‑CATTGTCCGAGCAT‑CCACCACTCGCGGGAACCCCACgg              >  1:89142/1‑228 (MQ=255)
                                                                                                                                                                                                                    aTCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCTGGCGATTAACGTACTGCACCGGGAATCG‑CATTGGCTGTGCAT‑CCCCCACTCGCTGGATATCCACGGTTTcc         >  2:21444/1‑228 (MQ=255)
                                                                                                                                                                                                                    aTCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTc          >  2:481791/1‑227 (MQ=255)
                                                                                                                                                                                                                    aTCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGAAATTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACTACCTGCACCCGGAACCG‑CATTGGCACA‑CATACCACCACTCGCTGGTtccacccgtgttc          >  1:92298/1‑216 (MQ=255)
                                                                                                                                                                                                                     tCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCCCCGCGACGCG‑CATTGGCTGGGCAT‑CCCCCACTCCCTGACTACCCCCGGTTTccc        >  1:46976/1‑227 (MQ=255)
                                                                                                                                                                                                                      tccc‑TTCCCCGGCAAAGGTTTCTTCGCGATCCCCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGCCACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCa        <  2:243018/225‑1 (MQ=255)
                                                                                                                                                                                                                      cG‑C‑TTCCCCGGCAAAGGCTTCTTCGCGATCCCCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCa        <  2:541391/227‑1 (MQ=255)
                                                                                                                                                                                                                      cG‑C‑TTCACCGGCAAAGGCTTCTGCGCGATCACCATCAAAAGCCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCa                 <  2:46976/218‑1 (MQ=255)
                                                                                                                                                                                                                             ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑cc                                                                                                                                                                          >  1:247456/1‑62 (MQ=255)
                                                                                                                                                                                                                             ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGCGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAACTCGCGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CAACCACTCGCTGGATCTCCCCGGGTTcaaaacctc  >  1:551590/1‑220 (MQ=255)
                                                                                                                                                                                                                             ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGCCTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGCTTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CAACCACTCGCTGGATCTCCACGGTTTcacgaactt  >  2:164457/1‑220 (MQ=255)
                                                                                                                                                                                                                              tcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑Cca                                                                                                                                                                         <  2:247456/62‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
TTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGG‑CTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAG‑CAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCG‑C‑TTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCG‑ACGTTTGATT‑CCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCG‑CATTGGCTGAGCAT‑CCACCACTCGCTGGATCTCCACGGTTTCCAGCACTT  >  NC_000913/2874285‑2874729

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