Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,959 A→C 18.8% V133G (GTA→GGA)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9590AC18.8% 122.5 / 4.6 71V133G (GTA→GGA) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/21);  new base C (2/12);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 2.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAA  >  NC_000913/989723‑990196
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
gggAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTAc                                                                                                                                                                                                                                                 <  1:992390/238‑1 (MQ=255)
     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACCAATTTGCTTTCCAGCGTGTCGGAGGGCTGGGCC‑AAC‑TC‑GGCCTGCGGTGAGCGCTAGCGCTGTACCActt                                                                                                                                                                                                                                            >  2:897698/1‑238 (MQ=255)
     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTAGACTAATTTGCTTTCCAGCGTGTC‑GAGGGGTTGCGC‑CAAGTC‑GGCCTGCGGGGAGCGAGAGCGATGTAcccattt                                                                                                                                                                                                                                           >  1:148373/1‑234 (MQ=255)
                    tCGTATTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGA‑CGTTTAACTATGGTTTCTTTTACTACGAGAAAAAAGCGGCGACGGGTGCCAAAACCCATCAG‑ATGGAAGCGCGACATATCAAGCCCGGGGGTTAACGCCGAGAGGGCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGTTGGCGC‑AACGCC‑GGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑cc                                                                                                                                                                                                                            <  2:270343/237‑1 (MQ=255)
                            gcggtggCCCACAAACCAGGATTCCGGTTGCAGA‑CGTTTAACGATGTTTTCTTTTACTTCGCGAAAAAAAGGGCGACGGGTCCCAAAATCCACCAG‑ATGGAAGCCCGACATACAAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGaa                                                                                                                                                                                                                  <  2:148373/233‑1 (MQ=255)
                             cTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCCGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGGGATGGACCATTTTACCG‑ATAACCG‑CCAGCAAAAGa                                                                                                                                                                                                                   >  1:1209956/1‑237 (MQ=255)
                              tGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GGGGGCTTGCGC‑CACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTTACTG‑ATAACCC‑CCCGCACAAGAACt                                                                                                                                                                                                                >  2:428772/1‑239 (MQ=255)
                               ggTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTGCTTTCCAGCGTGTC‑GGGGGGGTGCGC‑AAACTC‑GGCCTGCGGGGGGCGATAGCGCTGTCGCAATTTACTG‑ATAAACGCCCACCAAGGGAACt                                                                                                                                                                                                                >  2:44251/1‑239 (MQ=255)
                                  gCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGGACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCGACAAATCGACTAATTTTCTTTCCCGCCTGGT‑GAGGGCTTAAGC‑AAAGTC‑CGCCTGCGGTGAGCGAAAAGGGTTGACCCTTTTCCTG‑AGACACG‑CCAGCAAAGGGAATTcgc                                                                                                                                                                                                            >  2:649585/1‑237 (MQ=255)
                                  gCCCACAAACCAGGATTCCTGTGGCAGA‑CGCTTAACGATGTTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAG‑ATGGAAGCGCGACATATCAAGCCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTcc                                                                                                                                                                                                             <  1:1212552/238‑1 (MQ=255)
                                       caAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTTCAACGGGGC‑GAGGGCCTGCGC‑AACGCC‑GGCCCGCGGTGAGCGATAACGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCCGAgg                                                                                                                                                                                                       >  1:136552/1‑239 (MQ=255)
                                         aaCCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCAAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTACCGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGAt                                                                                                                                                                                                     <  1:428772/239‑1 (MQ=255)
                                          aCCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAAGTGACTAATTTGCTTTTCAACCGGTC‑GAGGGGCTGCGCAAACGTG‑GGCC‑GCGGGGGGTGGAAGAGGTGTACCATTTAACTG‑ATAACCC‑CCAGGAAAGAAACATCCCAGAGGaag                                                                                                                                                                                                    >  2:934169/1‑237 (MQ=255)
                                               gATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTTCAACGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGGCCGCGGGGGGCGATAGCGGTTTACCATTTCACCG‑ATAAACG‑CCCGCAAAGGAAATTTCCCGAGGGTGACtt                                                                                                                                                                                                >  2:519363/1‑238 (MQ=255)
                                                aTTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                                               <  1:484931/238‑1 (MQ=255)
                                                 tttccgTTGCAGGCCGCTTAAGGATGGTTTCTTCTACGTCGCGAAAAAACCGGCGACGGGTGCAAAAAGCCATAAG‑ATGGAAGGCCGAAATATCAAGACCGGGGTTTAACGCCGAAAAGTCGATTAATTTGTTTCCCAGCGTGTC‑GAGGGCTTGCGC‑ACCGTC‑GGCCTGCGGTGGGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                                               <  2:784158/233‑1 (MQ=255)
                                                    cTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTTCTTTTCAGCGGGTC‑GAGGGCTTGCGC‑AACGGC‑GGCCTGCGGTGAGCGGTAACGATGTAACATTTTACTT‑ATAAACG‑CCCGCAAAGGAACTTCCCCGAGGGTGACTTTACGc                                                                                                                                                                                           >  1:532900/1‑238 (MQ=255)
                                                    cTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCCAGAAGTCGACTCATTTTCTTTCCCGCGGGTC‑GGGGGGTTGCGC‑AACGCC‑CGCCCGCGGGGGGGGGTAGCGGTGTACCATTTTCATG‑ATAAACC‑CCCGCAAAAGGACATCCCAGGGGGAGACtttccgc                                                                                                                                                                                           >  2:182688/1‑233 (MQ=255)
                                                                           ttttCTTGTACTGCGGAAGAAAAAGGGCGAGGGGTGCCAAAATCCA‑CAGAAGGGAAGCGCGACATATCAAGACGGGCGTTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGA‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAatat                                                                                                                                                                     <  1:967773/237‑1 (MQ=255)
                                                                                 tGTCCTTGGGGAAAAAAGCGGGGACGGGGCCAAAAACCCATCAG‑ATGGAAGCCGGACATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTGGCTTTCCAGCGTGTC‑GAGGGCTTGCGA‑ACCGTC‑GGCCTGCGGGGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCtt                                                                                                                                                              <  2:961553/239‑1 (MQ=255)
                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTC‑GGGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCAATG‑ATAACCG‑CCCGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGc                                                                                                                                                                                    >  1:1235500/1‑210 (MQ=255)
                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCCGCGGTGGGCGATAGCGATGTAACATTTCACTG‑CTAAACG‑CCAGCAAAAGAACTTCCCCCAGGATGACTTCACGCCACCGGCCGCTTAAACGCATATCCAACTTTCCCttt                                                                                                                                                       >  1:1001106/1‑238 (MQ=255)
                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGGACAACAAATTGCTTTACAGAGGGTCGGAGGGC‑TGCGC‑GAAATC‑GGGCTTCAGGGGGCGATAGACATTTCACATTTCCTCG‑CAAAACG‑CCCGCCAAGGGACTTACCAGAAGAAGGGTTCCAGCCACCGGGCGCCTAAACCGATATTCCACTTGCCACTa                                                                                                                                                       >  2:210005/1‑239 (MQ=255)
                                                                                                aaccGGCGAGGGGGGCAAAAACCCATCAG‑ATGGAAGCGCGACATACCAAGACCGGCGGTAAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGa                                                                                                                                               <  1:296234/236‑1 (MQ=255)
                                                                                                              gcaaaaacacATCAG‑ATGGAAGCGCGACATATCAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCGGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCGGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  2:532900/230‑1 (MQ=255)
                                                                                                              gacAAAACCCATAAA‑GTGGAAGGCCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTCTTTTCCAGCGTGTC‑GGGGGTTGGCGC‑ACCGCC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  2:849192/237‑1 (MQ=255)
                                                                                                               ccAAAATCCATCAG‑AGAAAAGCGGAAAATATAAGGACCGGCTGTTAACGCGAATAAGTCGACTATTTGGCTTCCCGGCGTGTG‑TGGGGCTTGCGC‑ACCGTG‑CGCCGGCGGTGGTCGAGAGCGAGTTCCCATTTCACTG‑ATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  2:943213/238‑1 (MQ=255)
                                                                                                               ccAAAACCCATCAG‑ATGGAAGCGCGACATATCAAGGCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  1:519363/238‑1 (MQ=255)
                                                                                                                         tCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCCTTTCACTG‑ATAACCG‑CCCGCCAAGGAACTTCCCAAGAGATGACTTCACGCCAACGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCATTGGACACAGGGGCTTG‑TTCCGGGGTaaa                                                                                                                      >  2:572585/1‑238 (MQ=255)
                                                                                                                                ggAAGCGCGAAAGATCAAGACCGGGGGTTAACGCCGAGAAGCCGACTAATTTGCTTCCCAGCGTGTC‑GAGGTCTTGCCC‑AACGTC‑GGCCGGGGGTGAGCAATAGCGATGTACCTTTCCCCTG‑AAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGaa                                                                                                                <  2:1209956/238‑1 (MQ=255)
                                                                                                                                        aaCATATCAAGCCCGGGGGTTAACGCCGGGAGGTCGACTATTTTGTTTCCCAGCGTGTC‑GGGGGTTTGCGC‑AACGGC‑GGCCTGCGGGGGGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                                       <  2:289057/238‑1 (MQ=255)
                                                                                                                                         aCATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                                       <  1:44251/238‑1 (MQ=255)
                                                                                                                                            tatCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTTCTTTCCACCGTGTC‑GAGGGCTTGCGC‑CACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCCTTTCACTG‑ATAACCG‑CCCGCAAACGCACTTCCCCGAGGATGACTTCCCCCCACCGGCCGCTTAAACGAATATCCCGCCTGCCCTTATTCTTGTACACGGTGGCTTGTTTCCGGGTTAAA‑CGGACAACGCGGAACCac                                                                                                   >  2:689178/1‑238 (MQ=255)
                                                                                                                                               cAAGCCCGGGGTTTAACGTCGAAACGTCGGCTAGTTGGCTTCCCAGGGTGTG‑GAGGGTTTGCCC‑AAGGTC‑GGCGCGGGGGGACGGATAGCGTGGTCCCATTTCTCGA‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTa                                                                                                <  2:1226803/239‑1 (MQ=255)
                                                                                                                                                      gggggTAAACGCCGGGAAGGCGAATAATTTGCTTTCCAGCTGGTC‑GAGGGCATGCCA‑ACCGTC‑GGCTGGCGGTGACGAATAGCGATGTACAATTGCCCTG‑AGAACCG‑CCAGCAAAGGAACTGCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCGTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                                         <  2:157560/236‑1 (MQ=255)
                                                                                                                                                       gcggTTTAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTC‑GGGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCCTTTTACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCCCGGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCGTTGGACACGGGGGATTG‑TTCCGGGGTAAAGCGGAACTCCGGGAAACAGTTaaaaaaaac                                                                                        >  2:686493/1‑232 (MQ=255)
                                                                                                                                                       gcggTTAACGCCGAGAAGTCGACTAATTTTCTTTCCAGCGTGTC‑GGGGGCTTGCGC‑AACGGC‑GGCCTGCGGGGAGCGAGAGCGGTGTACCATTTTACTG‑ATAAACG‑CCCGCCAAGGAACTTCCCCGAGGGTGACTTTCCCCCACCGGCCGCTTAAACGAAAATACACCTTGCCCTTATCGTTTGACACGGTGACCTG‑TTTCGGGTTCAAGCGGAACTCCCGTAAACCGTTACGATAGt                                                                                         >  2:631283/1‑238 (MQ=255)
                                                                                                                                                        cggTTACCGCCGAGAAGTCGACTAATTTGCTTCCCAGGGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGAGATTTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                                        <  1:740722/238‑1 (MQ=255)
                                                                                                                                                        cggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCGGCGTTGC‑GTGTGCTTGCGA‑ACCGGC‑TGCGTGCGGGGATCGATAGCGATGTCCCATTTCCCTA‑AAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                                        <  2:136552/238‑1 (MQ=255)
                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                                           >  2:379034/1‑239 (MQ=255)
                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                                           >  2:116645/1‑239 (MQ=255)
                                                                                                                                                                      aaGTCGACTAATTGGCTTTCCAGCGTGTC‑GGGGGCTTGCCA‑ACCGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                                         <  1:901034/239‑1 (MQ=255)
                                                                                                                                                                       agtcatagaaTTTGTTTCCCACTGTGCC‑GAGGGTTTGCGC‑CCCGTG‑GGCCTGCGGGGAGAGTTAGTGATGTACCATTGCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                                         <  1:1158402/230‑1 (MQ=255)
                                                                                                                                                                        gTCGGCTAATTTGCTTTCGTTGGTGGG‑GAGGTCTTGCGC‑ATGGCC‑GCCCTGCGTGGAGGGAGAGCGATGTCCTTTCTCATTG‑ATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                                        <  2:379445/238‑1 (MQ=255)
                                                                                                                                                                         tCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag                                                                      >  2:1218794/1‑239 (MQ=255)
                                                                                                                                                                            aCTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcagggc                                                                   >  1:204556/1‑239 (MQ=255)
                                                                                                                                                                                   tgcttccgggCGTGTG‑GGGGGCTTGCG‑‑AACGGCGGGCGAGCGGTGGGGGATGGCGGTGTTCCATTTCGCTA‑ATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                                            <  2:20286/230‑1 (MQ=255)
                                                                                                                                                                                       ttCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                                           >  1:707889/1‑236 (MQ=255)
                                                                                                                                                                                          cAGTGTGTC‑GGGGGTTTCGGC‑AATCGC‑GGCCTGCGTGGGCCAATAGCTATGTACCATTTCACTT‑ATATCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                                     <  1:865714/239‑1 (MQ=255)
                                                                                                                                                                                              gtgtC‑GAGGGTTTGCGC‑AACGCC‑GGCCGGCGGGGAGCGATAGCGATGTACAATTTCACAG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                                                  <  1:416752/238‑1 (MQ=255)
                                                                                                                                                                                              gtggC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                                                  <  1:1218794/238‑1 (MQ=255)
                                                                                                                                                                                               tgtC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtca                                                >  1:819315/1‑239 (MQ=255)
                                                                                                                                                                                                           tGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  1:731040/1‑239 (MQ=255)
                                                                                                                                                                                                            gcgc‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  1:1091819/1‑238 (MQ=255)
                                                                                                                                                                                                            gcgc‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  2:249316/1‑238 (MQ=255)
                                                                                                                                                                                                                         gCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgct                          >  2:334793/1‑238 (MQ=255)
                                                                                                                                                                                                                           gggcgtgaggggatagcgatGT‑CCCTTTCACCGAAAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:548900/228‑1 (MQ=255)
                                                                                                                                                                                                                           cggCGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:1152536/237‑1 (MQ=255)
                                                                                                                                                                                                                           cgcgggggAGCGGAAGCGATTTAACTTTCCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:600998/232‑1 (MQ=255)
                                                                                                                                                                                                                           cTGGGGTGAGCGATACCGTTGTGCCATTTCACTA‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:680435/239‑1 (MQ=255)
                                                                                                                                                                                                                           cTGCGGTGGGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:182688/239‑1 (MQ=255)
                                                                                                                                                                                                                           cTGCGGTGAGGGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgctgc                        <  1:1117046/238‑1 (MQ=255)
                                                                                                                                                                                                                              cGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc                     >  2:678394/1‑238 (MQ=255)
                                                                                                                                                                                                                              cGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc                     >  1:467246/1‑238 (MQ=255)
                                                                                                                                                                                                                                   agcgatagcgatGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaa                                                                               >  1:490444/1‑175 (MQ=255)
                                                                                                                                                                                                                                        tagcgatGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGaa              >  1:988364/1‑235 (MQ=255)
                                                                                                                                                                                                                                           cgatGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTc       >  2:1020312/1‑239 (MQ=255)
                                                                                                                                                                                                                                            gttGTACCATTTCCCTG‑ACAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTc       <  2:386943/236‑1 (MQ=255)
                                                                                                                                                                                                                                             atGTACCTTTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                                   <  1:747319/161‑1 (MQ=255)
                                                                                                                                                                                                                                             atGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                                   >  2:747319/1‑161 (MQ=255)
                                                                                                                                                                                                                                                tACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGaa  >  1:290137/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 aCCATTTCACTG‑AAAACCG‑CCATCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGaa  <  2:669835/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
GGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAA  >  NC_000913/989723‑990196

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