BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,776,262	G→T	6.9%	S40Y (TCC→TAC) ‡	yibT ←	protein YibT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,776,262	0	G	T	6.9%	341.1 / 9.7	102	S40Y (TCC→TAC) ‡	yibT	protein YibT
Reads supporting (aligned to +/- strand): ref base G (49/46); new base T (2/5); total (51/51)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 109 total aligned reads displayed.

CTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCGGG  >  NC_000913/3776113‑3776410

ctgaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcgg                                                                                                                                                      >  1:604799/1‑150 (MQ=255)

 taaGCTGATTTTGGAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcgga                                                                                                                                                     >  1:314695/3‑150 (MQ=255)

  gaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcggac                                                                                                                                                    <  2:1295485/150‑1 (MQ=255)

    gCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcggacgg                                                                                                                                                  >  1:67103/1‑150 (MQ=255)

     cTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcggacgga                                                                                                                                                 <  1:1123334/150‑1 (MQ=255)

         tttCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGTCCGGAAtac                                                                                                                                             <  2:323465/150‑3 (MQ=255)

           tCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGTCCGGAAtaccc                                                                                                                                           <  1:985477/150‑5 (MQ=255)

           tCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGTCCGGAAtaccc                                                                                                                                           <  2:387273/150‑5 (MQ=255)

            ccAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGt                                                                                                                                          >  2:1563734/1‑150 (MQ=255)

             cAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGTCCGGAAtacccga                                                                                                                                         <  1:104555/150‑7 (MQ=255)

             cAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGtt                                                                                                                                         <  2:604799/150‑1 (MQ=255)

             cAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGtt                                                                                                                                         <  2:1017160/150‑1 (MQ=255)

               aaaTCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTAc                                                                                                                                       <  1:560341/150‑1 (MQ=255)

                   cAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACATGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACgggg                                                                                                                                   >  1:387848/1‑150 (MQ=255)

                     gCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGggag                                                                                                                                 >  1:627842/1‑150 (MQ=255)

                     gCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGggag                                                                                                                                 >  1:1331120/1‑150 (MQ=255)

                      cTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGggagg                                                                                                                                >  2:786026/1‑150 (MQ=255)

                      cTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGggagg                                                                                                                                <  2:164622/150‑1 (MQ=255)

                      cTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGggagg                                                                                                                                >  2:1254563/1‑150 (MQ=255)

                       ttAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGggagga                                                                                                                               >  1:164839/1‑150 (MQ=255)

                           aaaaaGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATATTACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACgg                                                                                                                                     <  1:294213/140‑1 (MQ=255)

                                gCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTc                                                                                                                      <  2:1313208/150‑1 (MQ=255)

                                 ccTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCa                                                                                                                     <  1:764588/150‑1 (MQ=255)

                                  cTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAg                                                                                                                    <  1:1477686/150‑1 (MQ=255)

                                   tttCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGa                                                                                                                   >  1:532264/1‑150 (MQ=255)

                                   tttCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGa                                                                                                                   <  1:1319300/150‑1 (MQ=255)

                                     tCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGCTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATa                                                                                                                 >  1:117421/1‑150 (MQ=255)

                                      gggcaGTAAGCAGACGT‑TGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATAc                                                                                                                >  1:1619053/5‑150 (MQ=255)

                                       ggcaGTAAGCAGACGT‑TGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACt                                                                                                               >  1:779271/4‑150 (MQ=255)

                                       aaaaGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTACGGACGGAATAGCGTTACGGGGAGGAAGTTTTTTCAGATACt                                                                                                               >  1:754573/1‑150 (MQ=255)

                                        gcaGTAAGCAGACGT‑TGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTc                                                                                                              >  1:1006347/3‑150 (MQ=255)

                                        aaaGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTc                                                                                                              >  2:957187/1‑150 (MQ=255)

                                          aGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTccc                                                                                                            >  2:1211773/1‑150 (MQ=255)

                                          aGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTccc                                                                                                            >  2:265451/1‑150 (MQ=255)

                                            tAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCgg                                                                                                          <  2:1042862/150‑1 (MQ=255)

                                                cA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGGCGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACg                                                                                                      >  1:1448486/1‑150 (MQ=255)

                                                      tCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTgg                                                                                                 >  1:573512/1‑150 (MQ=255)

                                                        tGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCt                                                                                               <  1:16366/150‑1 (MQ=255)

                                                             aCTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATg                                                                                          <  1:24292/150‑1 (MQ=255)

                                                               tGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGcc                                                                                        <  1:1422856/150‑1 (MQ=255)

                                                                gCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCa                                                                                       >  1:619784/1‑150 (MQ=255)

                                                                  cccTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATg                                                                                     >  2:171332/1‑150 (MQ=255)

                                                                     tctACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGaaa                                                                                  <  1:1254563/150‑1 (MQ=255)

                                                                       tACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATc                                                                                <  2:67103/150‑1 (MQ=255)

                                                                        aCCTGATTCGGCCTATAAAGACGACGATAATACTGCAGACTTTGTATATATAACGGCACGCTTTGATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCt                                                                               <  2:1043562/150‑1 (MQ=255)

                                                                         ccTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTa                                                                              >  2:593403/1‑150 (MQ=255)

                                                                         ccTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTa                                                                              >  2:822625/1‑150 (MQ=255)

                                                                           tGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACg                                                                            <  2:164839/150‑1 (MQ=255)

                                                                                    cGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCg                                                                   >  2:887584/1‑150 (MQ=255)

                                                                                    cGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCg                                                                   >  1:769668/1‑150 (MQ=255)

                                                                                    cGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCg                                                                   <  1:593403/150‑1 (MQ=255)

                                                                                      attaaGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGAt                                                                 <  2:1092121/147‑1 (MQ=255)

                                                                                       tAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATa                                                                <  1:379842/150‑1 (MQ=255)

                                                                                       tAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATa                                                                <  2:890834/150‑1 (MQ=255)

                                                                                           gCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATaagaa                                                            >  1:953360/1‑150 (MQ=255)

                                                                                              gacgaTAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCg                                                         >  2:894459/1‑150 (MQ=255)

                                                                                              gacgaTAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCg                                                         <  2:387848/150‑1 (MQ=255)

                                                                                               acgaTAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCgg                                                        <  1:1203961/150‑1 (MQ=255)

                                                                                                   taataCTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACTGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACg                                                    <  2:117421/150‑1 (MQ=255)

                                                                                                   taataCTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACg                                                    >  1:1118898/1‑150 (MQ=255)

                                                                                                      taCTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGttt                                                 <  2:1331120/150‑1 (MQ=255)

                                                                                                       aCTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGtttt                                                >  2:775665/1‑150 (MQ=255)

                                                                                                       aCTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGtttt                                                >  2:769952/1‑150 (MQ=255)

                                                                                                         tCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCa                                              <  2:1437672/150‑1 (MQ=255)

                                                                                                         tCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCa                                              <  1:948703/150‑1 (MQ=255)

                                                                                                         tCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCa                                              <  1:171332/150‑1 (MQ=255)

                                                                                                         tCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCa                                              >  1:640844/1‑150 (MQ=255)

                                                                                                          ccAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCAc                                             >  1:1171814/1‑150 (MQ=255)

                                                                                                             gACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGcc                                                                                        >  2:989004/1‑104 (MQ=255)

                                                                                                             gACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGcc                                                                                        <  1:989004/104‑1 (MQ=255)

                                                                                                                  tGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTg                                     <  2:627842/150‑1 (MQ=255)

                                                                                                                    tAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACTGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGcc                                   >  1:1282198/1‑150 (MQ=255)

                                                                                                                    tAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGcc                                   <  1:751252/150‑1 (MQ=255)

                                                                                                                     aGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGccc                                  >  1:562135/1‑150 (MQ=255)

                                                                                                                      gATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCa                                 <  2:497227/150‑1 (MQ=255)

                                                                                                                        tataACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCatg                               >  1:1485541/1‑150 (MQ=255)

                                                                                                                         ataACGGCACGCTTTCATCGGGTATTACGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCAGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAACCGGAACGTTATCACCTAATTTGCCCaaga                              >  2:1077752/1‑150 (MQ=255)

                                                                                                                          taACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCatgat                             <  1:1012058/150‑1 (MQ=255)

                                                                                                                          taACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCatgat                             >  1:1086959/1‑150 (MQ=255)

                                                                                                                          taACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCatgat                             >  1:1623010/1‑150 (MQ=255)

                                                                                                                          taACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTCACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCatgat                             <  2:1059179/150‑1 (MQ=255)

                                                                                                                           aaCGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATAGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCatgatg                            <  2:203083/150‑1 (MQ=255)

                                                                                                                            aCGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCAtgatga                           >  1:1229927/1‑150 (MQ=255)

                                                                                                                             cGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATgatgat                          >  2:157628/1‑150 (MQ=255)

                                                                                                                             cGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATgatgat                          >  1:1126715/1‑150 (MQ=255)

                                                                                                                              ggCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGata                         <  1:775665/150‑1 (MQ=255)

                                                                                                                              ggCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGata                         <  1:769952/150‑1 (MQ=255)

                                                                                                                                  cGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATaag                                                              >  1:556567/1‑109 (MQ=255)

                                                                                                                                  cGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATaag                                                              <  2:556567/109‑1 (MQ=255)

                                                                                                                                  cGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATAtctc                     >  1:203305/1‑150 (MQ=255)

                                                                                                                                  cGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATAtctc                     >  1:578877/1‑150 (MQ=255)

                                                                                                                                      ttCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTa                 >  1:1113624/1‑150 (MQ=255)

                                                                                                                                            gggTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCACATGATGATATATCCGTATTAccc           >  2:304844/1‑150 (MQ=255)

                                                                                                                                            cgcgaTTCCGTCCGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTAccc           >  1:1429706/5‑150 (MQ=255)

                                                                                                                                            cgcgaTTCCGTCCGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTAccc           >  2:428816/5‑150 (MQ=255)

                                                                                                                                              gTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCcg         <  2:421876/150‑1 (MQ=255)

                                                                                                                                                  tCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGcc     <  1:128165/150‑1 (MQ=255)

                                                                                                                                                   ccGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCg    >  2:723635/1‑150 (MQ=255)

                                                                                                                                                    cggacggaATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCgg   <  1:723635/150‑1 (MQ=255)

                                                                                                                                                     ggacggaATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCggg  >  1:116943/1‑150 (MQ=255)

CTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCGGG  >  NC_000913/3776113‑3776410

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