BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,573,096	A→G	pseudogene (683/1503 nt)	yjiT →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,096	0	A	G	99.2%	479.1 / NA	124	pseudogene (683/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (71/52); total (71/53)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

AGCTGTCATAGAAAAATCAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTG  >  NC_000913/4572950‑4573243

agCTGTCATAGAAAAATCAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGtt                                                                                                                                                 >  2:517330/1‑151 (MQ=255)

  cTGTCATAGAAAAATCAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTc                                                                                                                                                >  1:272115/1‑150 (MQ=255)

        tAGAAAAATCAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTAc                                                                                                                                         >  2:344820/1‑151 (MQ=255)

          gAAAAATCAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCt                                                                                                                                       >  4:222778/1‑151 (MQ=255)

             aaaTCAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGc                                                                                                                                    >  3:129692/1‑151 (MQ=255)

                 cAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGaa                                                                                                                                 <  4:32385/150‑1 (MQ=255)

                      cTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGt                                                                                                                           >  1:261428/1‑151 (MQ=255)

                      cTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGt                                                                                                                           >  4:423545/1‑151 (MQ=255)

                       tGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGtt                                                                                                                          >  2:1665138/1‑151 (MQ=255)

                        gCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGttt                                                                                                                         >  4:107681/1‑151 (MQ=255)

                                tttttAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATAcc                                                                                                                 >  4:1362002/1‑151 (MQ=255)

                                  tttAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCg                                                                                                                <  4:1205292/150‑1 (MQ=255)

                                   ttAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGtt                                                                                                              >  2:1569878/1‑151 (MQ=255)

                                   ttAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGt                                                                                                               <  2:1018415/150‑1 (MQ=255)

                                     aGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTa                                                                                                             <  2:284423/150‑1 (MQ=255)

                                     aGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTa                                                                                                             <  2:364856/150‑1 (MQ=255)

                                       cGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGtt                                                                                                              >  2:265726/1‑147 (MQ=255)

                                       cGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGtt                                                                                                              >  3:262206/1‑147 (MQ=255)

                                       cGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAgat                                                                                                          <  3:394656/151‑1 (MQ=255)

                                        gAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAgatg                                                                                                         >  4:1171065/1‑151 (MQ=255)

                                        gAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAgatg                                                                                                         >  2:1392598/1‑151 (MQ=255)

                                         aaGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAgatga                                                                                                        <  3:368625/151‑1 (MQ=255)

                                           gATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGatgatg                                                                                                      <  4:387651/151‑1 (MQ=255)

                                            aTACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGatgatg                                                                                                      >  3:93907/1‑150 (MQ=255)

                                            aTACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGatgatg                                                                                                      >  2:668446/1‑150 (MQ=255)

                                              aCCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGaa                                                                                                    <  2:1121914/150‑1 (MQ=255)

                                              aCCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGaa                                                                                                    >  2:212645/1‑150 (MQ=255)

                                               ccTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAAc                                                                                                  <  2:1392232/151‑1 (MQ=255)

                                                 tctGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAAc                                                                                                  >  2:635305/1‑149 (MQ=255)

                                                 tctGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAAc                                                                                                  >  4:1613473/1‑149 (MQ=255)

                                                 tctGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAAc                                                                                                  >  4:859909/1‑149 (MQ=255)

                                                   tGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAgg                                                                                               >  2:873033/1‑150 (MQ=255)

                                                    gTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAGCCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGaa                                                                                             <  2:955372/151‑1 (MQ=255)

                                                    gTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAg                                                                                                <  1:232235/148‑1 (MQ=255)

                                                      aGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTCCTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCCCCCTACCTGGCGAAGCGAGTTTCCAATCCCG‑TAGATGATGAAACAGGAAc                                                                                            <  4:1115940/149‑1 (MQ=255)

                                                      aGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAAc                                                                                            <  2:1219619/150‑1 (MQ=255)

                                                        aaCTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGc                                                                                          >  3:432698/1‑150 (MQ=255)

                                                        aaCTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGc                                                                                          >  4:629076/1‑150 (MQ=255)

                                                           ttATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCAttt                                                                                      >  2:1543389/1‑151 (MQ=255)

                                                             aTCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCAttttt                                                                                    >  4:835415/1‑151 (MQ=255)

                                                               cAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCAtttttt                                                                                   <  4:880166/150‑1 (MQ=255)

                                                               cAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTg                                                                                  <  2:907332/151‑1 (MQ=255)

                                                                aGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTg                                                                                  >  2:1353043/1‑150 (MQ=255)

                                                                aGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGa                                                                                 <  4:1381256/151‑1 (MQ=255)

                                                                      atGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGAc                                                                           >  2:1417707/1‑151 (MQ=255)

                                                                      atGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGAc                                                                           >  4:409554/1‑151 (MQ=255)

                                                                          gtGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCt                                                                       >  2:371743/2‑151 (MQ=255)

                                                                             aTAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATg                                                                    <  1:351940/151‑1 (MQ=255)

                                                                               aaCCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATgcg                                                                  >  4:1338118/1‑151 (MQ=255)

                                                                                 ccTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCa                                                                <  4:843997/151‑1 (MQ=255)

                                                                                   tCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGc                                                              >  2:993194/1‑151 (MQ=255)

                                                                                    cAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGcc                                                             <  2:6550/151‑1 (MQ=255)

                                                                                         ccTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                        >  2:436704/1‑151 (MQ=255)

                                                                                           ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                        >  1:88662/1‑149 (MQ=255)

                                                                                           ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                        >  2:467445/1‑149 (MQ=255)

                                                                                            tAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                        >  2:1579930/1‑148 (MQ=255)

                                                                                             aGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                    <  2:1013759/151‑1 (MQ=255)

                                                                                              gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAgg                                                   >  2:425195/1‑151 (MQ=255)

                                                                                                   tactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                               <  3:20298/150‑1 (MQ=255)

                                                                                                   tactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                               <  2:534697/150‑1 (MQ=255)

                                                                                                   tactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTa                                                >  1:310050/1‑149 (MQ=255)

                                                                                                       cacaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTc                                          <  3:261443/151‑1 (MQ=255)

                                                                                                        acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                         <  3:209218/151‑1 (MQ=255)

                                                                                                        acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                         >  4:1208261/1‑151 (MQ=255)

                                                                                                          aaCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGAc                                       >  2:1336442/1‑151 (MQ=255)

                                                                                                           aCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACt                                      >  4:1558432/1‑151 (MQ=255)

                                                                                                           aCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACt                                      <  4:711815/151‑1 (MQ=255)

                                                                                                            ccTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                     <  2:724917/151‑1 (MQ=255)

                                                                                                                ataaATCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaa                                 <  1:207440/151‑1 (MQ=255)

                                                                                                                  aaaTCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaa                               <  4:1569214/151‑1 (MQ=255)

                                                                                                                  aaaTCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaa                                <  4:232645/150‑1 (MQ=255)

                                                                                                                  aaaTCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaa                                >  4:669431/1‑150 (MQ=255)

                                                                                                                   aaTCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaa                               >  2:1623994/1‑150 (MQ=255)

                                                                                                                    aTCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaaaa                             <  4:660087/151‑1 (MQ=255)

                                                                                                                    aTCACAAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaaa                              >  3:359113/1‑150 (MQ=255)

                                                                                                                      cacaAAGAACCAGTACAACAGCTTGAGAGGGTTCGCCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTTAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAAt                            >  2:1249052/1‑150 (MQ=255)

                                                                                                                      cacaAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAAt                            <  3:41887/150‑1 (MQ=255)

                                                                                                                       acaAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATa                           >  2:1307221/1‑150 (MQ=255)

                                                                                                                        caAAGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATa                           >  1:414526/1‑149 (MQ=255)

                                                                                                                           aGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGt                      <  4:384930/151‑1 (MQ=255)

                                                                                                                           aGAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGt                      >  2:897858/1‑151 (MQ=255)

                                                                                                                            gAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTa                     <  2:1281630/151‑1 (MQ=255)

                                                                                                                            gAACCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTa                     <  2:734734/151‑1 (MQ=255)

                                                                                                                             aaCCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAc                    >  4:1234674/1‑151 (MQ=255)

                                                                                                                             aaCCAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAc                    >  2:778666/1‑151 (MQ=255)

                                                                                                                                cAGTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgcg                 >  2:1617149/1‑151 (MQ=255)

                                                                                                                                  gTACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGtt               <  2:1336443/151‑1 (MQ=255)

                                                                                                                                   tACAACAGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg              >  4:1441829/1‑151 (MQ=255)

                                                                                                                                       acaGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa          <  2:714065/151‑1 (MQ=255)

                                                                                                                                        caGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAt         <  3:157034/151‑1 (MQ=255)

                                                                                                                                        caGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAt         >  4:159659/1‑151 (MQ=255)

                                                                                                                                        caGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAt         >  4:1547037/1‑151 (MQ=255)

                                                                                                                                         aGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt        >  1:306450/1‑151 (MQ=255)

                                                                                                                                         aGCTTGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt        >  4:1050014/1‑151 (MQ=255)

                                                                                                                                            ttGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCt      >  1:6918/1‑150 (MQ=255)

                                                                                                                                            ttGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCt      >  3:44684/1‑150 (MQ=255)

                                                                                                                                             tGAGAGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTa     <  3:110325/150‑1 (MQ=255)

                                                                                                                                              gagaGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt   >  3:342151/1‑151 (MQ=255)

                                                                                                                                              gagaGGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt   <  2:1199074/151‑1 (MQ=255)

                                                                                                                                               agagGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg  <  4:102898/151‑1 (MQ=255)

AGCTGTCATAGAAAAATCAGCGCTGCCAACTGTTTTTAGCGAAGATACCTCTGTAGAACTTATCAGTCATATGGCTGATAACCTCAATTCCTTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTG  >  NC_000913/4572950‑4573243

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