Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,441,038 G→A intergenic (‑254/‑20) ydbK ← / → ydbJ putative pyruvate‑flavodoxin oxidoreductase/DUF333 domain‑containing protein YdbJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,441,0380GA83.0% 124.0 / 2.9 53intergenic (‑254/‑20)ydbK/ydbJputative pyruvate‑flavodoxin oxidoreductase/DUF333 domain‑containing protein YdbJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/26);  minor base T (9/0);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACGGCGATTGGGATGTGTG‑CATTACCCAACGGCAAACGCTG  >  NC_000913/1440796‑1441281
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
acagtGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGGGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑a                                                                                                                                                                                                                                                     <  1:585413/246‑1 (MQ=255)
aCATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑a                                                                                                                                                                                                                                                     <  2:63552/250‑1 (MQ=255)
                     aaCGAAATTGCTGCGCGCCCAGTAGTAATCTTTCACTTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCCGTTGGCTTAAATAT‑C‑AGATGGAGAAAAGAT‑ggcagcgg                                                                                                                                                                                                                              >  2:277093/2‑244 (MQ=255)
                        gAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑t                                                                                                                                                                                                                         <  1:94465/251‑1 (MQ=255)
                         aaaTTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTTAAGAT‑C‑AGAGGGAGAAAAGAT‑GCGAG‑CAGCggt                                                                                                                                                                                                                         >  1:333124/1‑249 (MQ=255)
                          aaTTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGGG‑CAGCG‑ttt                                                                                                                                                                                                                       >  1:77139/1‑251 (MQ=255)
                           aTTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCCGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGACGTTTGGCTTTACGAT‑C‑AGCGGAAAAAAAGAT‑GCGag‑ccggggttt                                                                                                                                                                                                                       >  1:325574/1‑243 (MQ=255)
                           aTTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACAACAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑tttt                                                                                                                                                                                                                      <  2:454489/251‑1 (MQ=255)
                              gctgcGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCACACAGAGGGCAAAAAAATAGCCATGAAAAAGTAAAAATGCTGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCTTTTTGCTTACATTT‑C‑AGAGGAATAAAAAAT‑TCGAG‑CAGCGGTTTTgg                                                                                                                                                                                                                    >  2:446031/1‑251 (MQ=255)
                                 gcgcgcCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTTCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGt                                                                                                                                                                                                                  <  1:535634/249‑1 (MQ=255)
                                  cgcgcCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTGAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAgg                                                                                                                                                                                                               <  1:326620/251‑1 (MQ=255)
                                     gcCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGgtg                                                                                                                                                                                                            <  2:365095/251‑1 (MQ=255)
                                                        aTTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCg                                                                                                                                                                                         <  1:242113/251‑1 (MQ=255)
                                                          tttAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGATTC‑A‑CGGGAGAAAAGAT‑GCGGGGCA‑CG‑GTTTGGGGAGGGGGTGCCGCTTTTTtgttgttggg                                                                                                                                                                                       >  2:483823/1‑247 (MQ=255)
                                                                   tGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑GTTTGGGTAGGGGGTGCCGCTTTATTGTTGTCGGCggggagt                                                                                                                                                                                >  1:266861/1‑243 (MQ=255)
                                                                    ggCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTagtag‑t                                                                                                                                                                            <  1:179522/251‑1 (MQ=255)
                                                                          cttctGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGGTTGGCTTAAAGAT‑C‑AGAGGGAGAAAAGAT‑GCGAG‑AAGCG‑GTTTGGGGTGGGGGTGCCGCTTTATTTTTGTCGGGGGgtagttaggaaact                                                                                                                                                                       >  2:20751/1‑242 (MQ=255)
                                                                                     tCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAgg                                                                                                                                                            <  1:516501/250‑1 (MQ=255)
                                                                                             gtgtCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTAAAGAT‑C‑CGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑GTTTGGGGAGGGGGTGCCGCCTTATTGTTGTTGGCGGGTGGGAG‑TGGACCTTTTTAGCAGGCGACTGccc                                                                                                                                                   >  1:560895/1‑250 (MQ=255)
                                                                                               gtCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAACAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTTAAGCT‑CCAG‑GGAAGAAAAAATTGC‑AG‑GCGGG‑TTTTGGGTAGGGGGTGCCGCCTTTTTTTTGTTGGCGGGGTGTGG‑TGAACCTGGTCAGCAGGCggacgccgcg                                                                                                                                                 >  2:495594/1‑242 (MQ=255)
                                                                                               gtCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGGCGGTTGGGTTTAAGAT‑T‑AGAGGGAGAAAAGAT‑GCGTG‑GAGCG‑TTTTGGGTAGGGGGTGCCGCTTTTTTTTTGTCGGGGGGTAGGAG‑TGGACCTGTTCCGGAGGGGACTGCcgcg                                                                                                                                                 >  1:375002/1‑251 (MQ=255)
                                                                                                tCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAA‑GCGGG‑CAGCG‑TTTTGGGGAGGGTGTGCCGCTTTTTTGTTGTCGGCGGGGAGGAG‑TGTAACTGTTCAGCAGGCGACTGCCgcgc                                                                                                                                                >  1:42955/1‑251 (MQ=255)
                                                                                                   cccACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACg                                                                                                                                             <  2:568500/251‑1 (MQ=255)
                                                                                                                aatgaatgTTTTGATCAAACAGAGGGCCAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGGTTGGCTTTAAGAT‑T‑CGAGGAAGAAAAAAA‑GCGAG‑CAGCG‑GTTTGGGGGGGGGGTGCCGCTTTTTTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGCGGACTGCCGCGCCCGCAGCGCCACG‑ttt                                                                                                                               >  2:157547/1‑251 (MQ=255)
                                                                                                                   gaatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTaa                                                                                                                             <  2:178881/250‑1 (MQ=255)
                                                                                                                    aatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGGGGCGGTTGGCTTTAAAAT‑C‑AGAGGGAGAAAAAAT‑GCGGG‑CAGGG‑GTTTTGGTAGGGGGTGCCGCTTTATTGTTTTCGGCGGGTAGTCT‑GGAACCTGTTTAGGATGCGCCTGCCGCGCCCGTAGCGCCAGG‑GTTaaa                                                                                                                            >  2:188556/1‑250 (MQ=255)
                                                                                                                     atgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAg                                                                                                                          <  1:583628/251‑1 (MQ=255)
                                                                                                                      tgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGGTTTGCTTTAAGAT‑C‑AGAGGGAGAAAAGAT‑GCGAG‑GAGCGGTTTTGGG‑GGGGGGTGCCGCTTTTTTGTTGTTGGCGGGTTGGAG‑TGGACCTGGTCAGCAGGCGACTGCCGCGCCCGtgaggcccgggttcaa                                                                                                                             >  2:182621/1‑233 (MQ=255)
                                                                                                                         tttGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTTTAGTTTGGCTTTAAGAT‑T‑AGAGGAAGAAAAGAT‑GCGAG‑GAGGG‑TTTTGGGGAGGGTGTGCCGCTTTATTTGTGGCGGAGGGTAGGAG‑TGAACCTGTTCAGGAGGGGGCTGCCGCGCCCGCAGCGCCAGG‑TTTAAAAGCGt                                                                                                                       >  1:252201/1‑250 (MQ=255)
                                                                                                                         tttGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGt                                                                                                                       <  1:364300/250‑1 (MQ=255)
                                                                                                                                 aaCAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAA‑GCGAG‑GAGCG‑TTTTTGGTAGGGGGTGCCGCTTTTTTGTTGTCGGCGGGGTGTAG‑GGAACCTGTTCAGCAGGGGACTGCCGCGCCCGTAGGGGCAGGGTTTAAAA‑CGTCGATGTccc                                                                                                              >  1:63243/1‑250 (MQ=255)
                                                                                                                                  aCAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCag                                                                                                             <  2:316422/251‑1 (MQ=255)
                                                                                                                                          cAAAAAAAGAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCACTGTGCTGCCTAGTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGGGGGCCGCTTTATTGTTGGCGGCGTGGAGTAG‑TGAACCTGTGCAGTAGGCGACTGCCGCGCAAGTAGCGCCAGG‑GTTAAAAGCGGCGATGTCGAGTAGTGgag                                                                                                     <  2:549207/251‑1 (MQ=255)
                                                                                                                                                         aaaaaaGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGGCGCTTTATTGTTGTCGGCGTGGTGTAGGTGA‑CCTGTTTAGCAGGCGGCTGCCGCGCACGTAGCGCCCGG‑TTTTAAAAAGGCGGTGTGCAGTAGTGGAGAAGCAAATTGTGcc                                                                                       >  1:52976/1‑249 (MQ=255)
                                                                                                                                                         aaaaaaGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCaa                                                                                      <  2:77139/251‑1 (MQ=255)
                                                                                                                                                              aGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGGGTCGGTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGGG‑GAGCGATTTTGGG‑AGGGTGGGCCGCCTTATTGTTTTCGGGGTGTAGTAG‑TGAACCTGGTTAGCAGGCGACTGCCGCGCACGTAGGGCCAGG‑TTTTAAAACATCGGTGTTCAGTAGTGGGGAAGCAAATTGTTCaaatgat                                                                                 >  1:541647/1‑246 (MQ=255)
                                                                                                                                                                    aaTGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑t                                                                          <  2:28801/251‑1 (MQ=255)
                                                                                                                                                                            gCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTAGTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTc                                                                  <  1:547693/251‑1 (MQ=255)
                                                                                                                                                                             cAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGGTTAAAGAT‑C‑CGAGGGAGAAAAGAT‑GGGGG‑CAGCG‑TTTTTGGTAGGGGGGGCCGCTTTTTTGTTGTCGGGGGGTAGTAG‑TGAAACTGTTTAGCAGGGGCCTGCCGCCCCCGGAGCGCCCGG‑GTTAAAAGCGTCGATGTCCCGTAGGGGAGAAGCAAATTGTGCAATGATcggcgggtggcctttt                                                                  >  2:31890/1‑241 (MQ=255)
                                                                                                                                                                              aaCGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGGTTGGCTTAAAGAT‑C‑CGAGGGAGGAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGGGGGGCGCTTTATTGTTGTCGGCGTGGAGTAG‑TGAACCTGTTCAGCAGGCGCCCGCCGCGCACGTAGCGCCAGG‑TTTAAAAACGTCGAGGTCCAGTAGTGGAGAAGGAAATTTGGCAATGATcggggggttgcgttttt                                                                 >  2:271576/1‑237 (MQ=255)
                                                                                                                                                                              aaCGCGCAGGACAGCTCTAGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGACGGTTGGCTTAAAGAT‑C‑AGAGGGAGGAAAGAT‑GCGAA‑CAGCG‑GTTTGGGTAGGGGGTGGCGCTGTATTGTTGTCGGAGGGGTGGAG‑GGGACCTGTTCAGCAGGCGCATGCCGCGCACGTAGCGCCAGG‑GTTAAAAGCGGCGACGTCCAGGAGGGGGGAAACAAATTGGGAAAAGATCGGCGG‑GTGGCtttttt                                                                >  1:464174/1‑248 (MQ=255)
                                                                                                                                                                                   gcAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑GTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑GGAACCTGTTCAGCAGGCGACTGCCGCGCACGGAGCGCCAGGGTTTAAAA‑CGTCGATGTCCAGTAGGGGGGAAGCAAATTGTGGAATGATCGGCGG‑GTTGCttttttttttc                                                           >  1:314728/1‑243 (MQ=255)
                                                                                                                                                                                             tctcGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTTAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGGGTGGCGCTTTATTGTTGTCGGCGGGTAGGAG‑GGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCCGG‑GTTAAAAGCGGCGATGTCCAGTAGGGGAGAAGCAAAATGTGCAATGATCGGCGG‑GTTGCtttttttttgcccgccacagg                                                 >  1:337591/1‑233 (MQ=255)
                                                                                                                                                                                                cGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑TCGAG‑CAGCG‑TTTTGGGGTGGGTGTGCCGCTTTTTTGTTTTCGGCGTGTAGGAG‑TGGAACTGTTCAGCAGGCGACTGCCGCGCACGGAGCGCCCGG‑GTTTAAAGAGTCGATGTCCCGTAGTGGGGAAGAAAAATGTGCAATGATCGGCGGGTTCCCTTTTTTTTGccccccacctggtg                                               >  2:430712/1‑240 (MQ=255)
                                                                                                                                                                                                cGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGAt                                               <  1:405785/250‑1 (MQ=255)
                                                                                                                                                                                                     ccgccCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACg                                         <  2:570069/251‑1 (MQ=255)
                                                                                                                                                                                                       gccCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACGGc                                       <  1:7006/251‑1 (MQ=255)
                                                                                                                                                                                                                       attattTATGGCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGAGCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACGGCGATTGTGATGTGTG‑c                       <  2:342447/250‑1 (MQ=255)
                                                                                                                                                                                                                                ggCGGTGTCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACGGCGATTGGGATGTGTG‑CATTACCCaa              <  1:535433/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gtgtCGTTTGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACGGCGATTGGGATGTGTG‑CATTACCCAACGGCa         <  2:333124/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tttGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTTGCTTTCTGTTGCCCGTCAACTGGATGGTACGGGGGTTGGGGTGTGGG‑GATTACCCCACGGGaaaagcc   >  1:300523/1‑247 (MQ=255)
                                                                                                                                                                                                                                           ttGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGGCCGGCAAATGGGTGGGAAGGCGCTTGGGATG‑GTGGCATTACCCCACGGCCAACGccg  >  2:86190/1‑249 (MQ=255)
                                                                                                                                                                                                                                           ttGGCTTAAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACGGCGATTGGGATGTGTG‑CATTACCCAACGGCAAACGCt   <  2:560895/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
ACATTGCGCAAATGAGGGGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGAT‑C‑AGAGGAAGAAAAGAT‑GCGAG‑CAGCG‑TTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCGTGTAGTAG‑TGAACCTGTTCAGCAGGCGACTGCCGCGCACGTAGCGCCAGG‑TTTAAAAGCGTCGATGTCCAGTAGTGGAGAAGCAAATTGTGCAATGATCGGCGG‑TTCGCTTTCTGTTGCCCGTCAACTGGATGGTACGGCGATTGGGATGTGTG‑CATTACCCAACGGCAAACGCTG  >  NC_000913/1440796‑1441281

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