BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,950,467	2 bp→TT	intergenic (+47/‑39)	ilvL → / → ilvX	ilvG operon leader peptide/uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,950,467	0	G	T	94.1%	423.7 / 13.1	120	intergenic (+47/‑40)	ilvL/ilvX	ilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (3/4); new base T (52/59); total (57/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01
*	NC_000913	3,950,468	0	G	T	99.2%	466.6 / NA	120	intergenic (+48/‑39)	ilvL/ilvX	ilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base T (57/62); total (57/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

ATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAACAC  >  NC_000913/3950322‑3950612

||

aTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑tttttt                                                                                                                                              >  2:670302/1‑151 (MQ=255)

      gCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGa                                                                                                                                        >  2:53648/1‑151 (MQ=255)

      gCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGGGGTGATTATTATCCCACCGTCCGGGGCTGCTCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCTACACGTCCGGG‑TTTTTTTTTTGa                                                                                                                                        >  1:1133084/1‑151 (MQ=255)

        ccTTCTACGAGTGATGAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGAcc                                                                                                                                      <  1:253681/151‑1 (MQ=255)

          ttCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCtt                                                                                                                                    >  2:1251775/1‑151 (MQ=255)

          ttCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGGCCGGG‑TTTTTTTTTTGACCtt                                                                                                                                    >  1:87011/1‑151 (MQ=255)

               cGAGGGATTAGCCCGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCCCTTGGACTAGGAAAGTTTTAGAGATCAAGCTTTACCGAACTAAAGCCCCCGCACCCAAAGGTCCGGT‑TTTTTTTTTTGACCTTaaaaa                                                                                                                               <  1:1243013/151‑1 (MQ=255)

                   tGATTAGCCTGGTCGTGTTTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGCCCGGG‑ATTTTTTTTTGACCTTAAAAACCTa                                                                                                                           >  1:476243/1‑151 (MQ=255)

                   tGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATa                                                                                                                           >  2:1388616/1‑151 (MQ=255)

                     aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAAc                                                                                                                         <  2:1275846/151‑1 (MQ=255)

                     aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCAAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAAc                                                                                                                         <  2:635126/151‑1 (MQ=255)

                      ttAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAAcc                                                                                                                        <  2:1032922/151‑1 (MQ=255)

                              gTCGTGATTAGCGTGGTGGTGCTTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTTACCTTAAAAACCTAACCGAGGAGCa                                                                                                                >  1:922404/1‑151 (MQ=255)

                              gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                                                                >  2:236181/1‑151 (MQ=255)

                                  tGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGAc                                                                                                             <  2:642281/150‑1 (MQ=255)

                                   gATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                                                                           >  1:1413245/1‑151 (MQ=255)

                                           tggtggtgATTATTATCCCACGGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                                                                   <  1:1115000/151‑1 (MQ=255)

                                           tggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                                                                   <  2:257863/151‑1 (MQ=255)

                                           tggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                                                                   <  1:1277789/151‑1 (MQ=255)

                                           tggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                                                                   <  2:736844/151‑1 (MQ=255)

                                             gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGACCTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                                                                 >  1:593348/1‑151 (MQ=255)

                                             gtggtcATTCTTATCCCACCGTGCGTGGCAGCACTTGGACGAGGAAAGGCTTTGAGATCAAGCCTTAACGAACCAAGACCCCCGCACCGCAAGGTCCGGG‑TTTTTTTTTGTCCCTTCAAAAC‑TAATCCATGAAAAGACAATGATTAAcc                                                                                                  >  1:503454/1‑148 (MQ=255)

                                              tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                                                                 <  2:959202/150‑1 (MQ=255)

                                              tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                                                                 <  2:1480171/150‑1 (MQ=255)

                                               ggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAAACAATGAATAACAGca                                                                                               >  1:929697/1‑151 (MQ=255)

                                                    ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGCGGAGCAGACAATGAATAACATCACaaaa                                                                                          >  1:750525/1‑151 (MQ=255)

                                                    ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                                                                          >  1:245275/1‑151 (MQ=255)

                                                    ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                                                                          >  2:330371/1‑151 (MQ=255)

                                                    ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                                                                         >  1:1012842/1‑151 (MQ=255)

                                                      attatCCCACCGTGCGGGGCTGGACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTCAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTATGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                                                                                        >  1:591732/1‑151 (MQ=255)

                                                      attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                                                                                        >  1:1451512/1‑151 (MQ=255)

                                                       ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGGTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                                                                                        >  2:1148862/1‑150 (MQ=255)

                                                       ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                                                                                        >  1:797322/1‑150 (MQ=255)

                                                       ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                                                                                        >  1:1467657/1‑150 (MQ=255)

                                                       ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                                                                                        >  2:1421320/1‑150 (MQ=255)

                                                         aacccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGTCTTAACGAACTAAGACCCCCGAACCGAACGCTCCGGG‑TG‑TTTTTTTGACCTTAAAAACATAACCGGGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                                                                    <  1:1089274/149‑1 (MQ=255)

                                                              aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                                                                  >  2:1146087/1‑149 (MQ=255)

                                                               ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                                                                  >  2:129376/1‑149 (MQ=255)

                                                                 gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACATACCCGAGGAGCAGACAATGAATAACAGCACAACATTCTGTTtctaa                                                                              >  1:956026/1‑149 (MQ=255)

                                                                 gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                                                                             >  1:556006/1‑151 (MQ=255)

                                                                 gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                                                                             >  1:1282142/1‑151 (MQ=255)

                                                                 gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                                                                             >  1:248113/1‑151 (MQ=255)

                                                                 gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTAAAAACATAACCGGGGGGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                                                                            >  2:1403973/1‑151 (MQ=255)

                                                                      gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt                                                                         >  2:1418337/1‑151 (MQ=255)

                                                                      gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt                                                                         >  1:1328557/1‑150 (MQ=255)

                                                                      gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt                                                                         >  2:935651/1‑150 (MQ=255)

                                                                       ggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                                                                        <  1:1475807/150‑1 (MQ=255)

                                                                            cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                                                                   <  1:1385690/151‑1 (MQ=255)

                                                                            cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                                                                   <  1:758916/150‑1 (MQ=255)

                                                                            cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                                                                   <  1:218179/150‑1 (MQ=255)

                                                                            cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                                                                   <  2:766879/150‑1 (MQ=255)

                                                                            cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg                                                                  >  2:558756/1‑151 (MQ=255)

                                                                                tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa                                                              <  2:621859/151‑1 (MQ=255)

                                                                                         aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta                                                      <  1:525362/150‑1 (MQ=255)

                                                                                            ggtttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                                  >  2:988012/4‑151 (MQ=255)

                                                                                            ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                                  >  2:422209/1‑151 (MQ=255)

                                                                                            ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                                  >  1:681156/1‑151 (MQ=255)

                                                                                            ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                                  >  1:489549/1‑151 (MQ=255)

                                                                                            ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                                  >  1:873220/1‑151 (MQ=255)

                                                                                            ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                                  >  1:1078926/1‑151 (MQ=255)

                                                                                            ggCTTAGAGATCAAGCCTCAACGGACTGAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                                  >  2:107719/1‑151 (MQ=255)

                                                                                             gCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAAt                                                 >  2:1077429/1‑151 (MQ=255)

                                                                                             gCTTAGAGAGCAAGCCTTAACGAACTAAGACGCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGa                                                   <  1:1424732/149‑1 (MQ=255)

                                                                                              cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAAt                                                 <  2:351856/150‑1 (MQ=255)

                                                                                              cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAAt                                                 <  2:1246126/150‑1 (MQ=255)

                                                                                              cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAAt                                                 <  1:362636/150‑1 (MQ=255)

                                                                                               ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATgg                                               <  1:236113/151‑1 (MQ=255)

                                                                                               ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATgg                                               <  2:830295/151‑1 (MQ=255)

                                                                                               ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATgg                                               <  1:1247235/151‑1 (MQ=255)

                                                                                               ttAGAGATCAAGCCGTAACGAACTAAGGCCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAGCCGAGGAGCAGACAATGAATAACAGCAGAAAATTCTGTTTCTCAAGATTCAGGCCGGGGAACTAACTAGGAATgg                                               <  1:964357/151‑1 (MQ=255)

                                                                                                    gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGca                                           >  1:483300/1‑150 (MQ=255)

                                                                                                    gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGCCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACACCACAACATTCTGTTTTTCAGGATTCAGGACGGGGAACTAACTATGAATggtgta                                           >  1:663702/1‑146 (MQ=255)

                                                                                                     aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGcac                                          <  1:1352788/151‑1 (MQ=255)

                                                                                                      tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAg                                        <  1:1060972/151‑1 (MQ=255)

                                                                                                      tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAgt                                       <  2:166483/151‑1 (MQ=255)

                                                                                                          gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAgtggg                                    >  2:1005817/1‑151 (MQ=255)

                                                                                                             ttAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGggtggt                                <  2:841731/151‑1 (MQ=255)

                                                                                                              tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGggtggt                                <  2:1468426/151‑1 (MQ=255)

                                                                                                              tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGggtggt                                <  1:1104074/151‑1 (MQ=255)

                                                                                                                   aaCTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATg                           >  2:469634/1‑151 (MQ=255)

                                                                                                                   aaCTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATg                           <  1:550343/151‑1 (MQ=255)

                                                                                                                    aCTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATg                           >  2:942116/1‑150 (MQ=255)

                                                                                                                    aCTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATg                           >  2:1437530/1‑150 (MQ=255)

                                                                                                                     cTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCg                         <  2:794056/151‑1 (MQ=255)

                                                                                                                     cTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCg                         <  1:1029423/151‑1 (MQ=255)

                                                                                                                     cTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAAACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCg                         <  2:930178/151‑1 (MQ=255)

                                                                                                                          aCCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTg                      <  1:1129581/149‑1 (MQ=255)

                                                                                                                           cccccGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCgg                   >  1:447941/1‑151 (MQ=255)

                                                                                                                               cGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGcac               <  2:338952/151‑1 (MQ=255)

                                                                                                                                 cACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAg             <  1:929767/151‑1 (MQ=255)

                                                                                                                                  aCCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAgg            <  2:175105/151‑1 (MQ=255)

                                                                                                                                   ccGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAggg           <  2:1007127/151‑1 (MQ=255)

                                                                                                                                       aaGGTCCGGG‑TTTTTTTTTTGACCTTAAAATCATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTAAGGACGGGGAACTAACTATGAATGGCGCACAGGGGGTGGTACATGCGTTGCGGGCGCAGGGtgtg       <  1:1365337/151‑1 (MQ=255)

                                                                                                                                          gTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGaa     >  2:434340/1‑150 (MQ=255)

                                                                                                                                           tCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAac    <  1:1398353/150‑1 (MQ=255)

                                                                                                                                           tCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAac    <  2:1326027/150‑1 (MQ=255)

                                                                                                                                           tCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAac    <  1:848180/150‑1 (MQ=255)

                                                                                                                                           tCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAac    <  1:1225939/150‑1 (MQ=255)

                                                                                                                                           tCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAac    <  1:558512/150‑1 (MQ=255)

                                                                                                                                           tCCGGG‑TTTTTGTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAac    <  1:428944/150‑1 (MQ=255)

||

ATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCACAGGGTGTGAACAC  >  NC_000913/3950322‑3950612

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