BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,776,263	A→C	6.9%	S40A (TCC→GCC) ‡	yibT ←	protein YibT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,776,263	0	A	C	6.9%	302.7 / 13.2	102	S40A (TCC→GCC) ‡	yibT	protein YibT
Reads supporting (aligned to +/- strand): ref base A (49/46); new base C (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 = 5.61e-01

Only 100 of 108 total aligned reads displayed.

TGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCGGG  >  NC_000913/3776114‑3776410

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                                                                                                                                           <  2:387273/150‑5 (MQ=255)

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

          tCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGTCCGGAAtaccc                                                                                                                                           <  1:985477/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‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACgggg                                                                                                                                   <  1:393014/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                   ttCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGAt                                                                                                                  >  2:1273393/1‑150 (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)

                                                                    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                                                                   <  1:593403/150‑1 (MQ=255)

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

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

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

                                                                                        aaGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATaag                                                              <  1:1563734/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                                              <  1:948703/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                                    aGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCC‑GGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCa                                 <  1:1211773/150‑1 (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)

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

                                                                                                                          aaCGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATAGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCatgatg                            <  2:203083/150‑1 (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                                                              <  2:556567/109‑1 (MQ=255)

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

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

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

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

                                                                                                                                           gggTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTAccc           <  1:957187/150‑1 (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)

TGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAACGCCTGGCTTGATGCCATGAAATCTACGGCTTTATCGATAAGAAGCGGAACGTTTTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCGGG  >  NC_000913/3776114‑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: `‑`