BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,213,234	A→C	5.5%	M1M (ATG→CTG) †	yjaA →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,213,234	0	A	C	5.5%	1788.9 / 28.7	531	M1M (ATG→CTG) †	yjaA	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base A (258/234); new base C (3/35); total (262/269)
Fisher's exact test for biased strand distribution p-value = 2.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 537 total aligned reads displayed.

ACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTT  >  NC_000913/4213085‑4213383

aCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCa                                                                                                                                                       <  2:3911284/150‑1 (MQ=255)

  tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATg                                                                                                                                                     >  2:6491143/1‑150 (MQ=255)

  tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATg                                                                                                                                                     >  2:5807548/1‑150 (MQ=255)

   cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                    <  1:7106036/150‑1 (MQ=255)

   cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                    <  2:1275642/150‑1 (MQ=255)

   cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                    <  2:3106891/150‑1 (MQ=255)

   cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                                                                                    <  2:3819332/150‑1 (MQ=255)

    aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGta                                                                                                                                                   >  1:3395531/1‑149 (MQ=255)

     gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCa                                                                                                                                                  >  1:1273264/1‑150 (MQ=255)

     gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCa                                                                                                                                                  <  2:8068286/150‑1 (MQ=255)

       aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGt                                                                                                                                                >  1:6248102/1‑150 (MQ=255)

         tGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTc                                                                                                                                              <  2:5734793/150‑1 (MQ=255)

                 cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAAACTCCTGTCAGTTCTGTATATc                                                                                                                                      <  1:6680830/150‑1 (MQ=255)

                  gTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTATATcc                                                                                                                                     <  1:4329483/150‑1 (MQ=255)

                  gTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATcc                                                                                                                                     >  1:5691527/1‑150 (MQ=255)

                          gATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAATGGAATCTCATGTCAGTTCTGTATATCCAAATtcgt                                                                                                                             <  2:2721991/150‑1 (MQ=255)

                          gATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgt                                                                                                                             <  1:372040/150‑1 (MQ=255)

                           aTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCCTAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                                                                            <  2:4468882/150‑1 (MQ=255)

                           aTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                                                                            <  2:4555058/150‑1 (MQ=255)

                                aGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAAt                                                                                                                       >  1:2683951/1‑150 (MQ=255)

                                  tgtgGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAATTCTGTATATCCAAATTCGTCGTAATCa                                                                                                                     >  2:6139686/1‑150 (MQ=255)

                                    tgGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCaaa                                                                                                                   <  2:5280024/150‑1 (MQ=255)

                                    tgGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCaaa                                                                                                                   >  1:1483479/1‑150 (MQ=255)

                                     ggggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAt                                                                                                                  <  2:872855/150‑1 (MQ=255)

                                      gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAtt                                                                                                                 >  2:1008783/1‑150 (MQ=255)

                                       ggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTa                                                                                                                <  2:4697274/150‑1 (MQ=255)

                                        gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTAc                                                                                                               >  2:7856647/1‑150 (MQ=255)

                                           tcCACATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACACCGGACCATAAGGCCGGGGTGTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTAAATCCAAATTCCTCGTAATCAAATTACTGt                                                                                                            >  2:5661895/1‑150 (MQ=39)

                                              ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg                                                                                                         >  1:3991950/1‑150 (MQ=255)

                                                  gCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGAt                                                                                                     <  1:6201513/150‑1 (MQ=255)

                                                  gCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTAGTTGCAGGAGGAAGGGAATCTCATGTCCGTTCTGTATATCCAAATTCGTCGGAATCAAATTACGCTTCGCGAt                                                                                                     >  2:5009154/1‑150 (MQ=255)

                                                    gagaGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCt                                                                                                   <  2:3269646/150‑1 (MQ=255)

                                                        gTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCAATAAGGCCGGGGTTTTTTGAATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaa                                                                                               <  2:2550437/150‑1 (MQ=255)

                                                        gTAGGGAACTGCCAGGCATCAAATTAGAAAAAACCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaa                                                                                               <  2:3369710/150‑1 (MQ=255)

                                                        gTAGGGAACTGCCAGGCATCAAAATAGAAAAACCCCGGTCCATAACGCCGGGGATTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaa                                                                                               <  1:7887114/150‑1 (MQ=255)

                                                         tAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                                                                              <  1:3645819/150‑1 (MQ=255)

                                                         tAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                                                                              <  1:4671774/150‑1 (MQ=255)

                                                         tAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                                                                              <  2:1053879/150‑1 (MQ=255)

                                                           gggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGc                                                                                            >  1:5575929/1‑150 (MQ=255)

                                                            ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCa                                                                                           >  1:2762775/1‑150 (MQ=255)

                                                            ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCa                                                                                           >  2:1847025/1‑150 (MQ=255)

                                                                   ccAGGCATCAAATTAGAAAAACCCCGGTCCATAATGCCGGGGTTTTTAGAATATCAATTATTTGCATGATGAAGGGAATCTCTTGTCAGTTCGGTATATCCAAATTCGTAGTAATCAAATTACTGTTCTCGATATTGAAAGCAAACGTGa                                                                                    >  1:6176445/1‑150 (MQ=255)

                                                                   ccAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGa                                                                                    <  2:3444715/150‑1 (MQ=255)

                                                                    cAGTCATCAAATTAGAAAAACCACTGTCCATAAGGCCTGGGTTTTTTGCATATCATTTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATATTGAAAGCAAACGTGaa                                                                                   <  2:4851873/150‑1 (MQ=255)

                                                                     aGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAg                                                                                  >  1:6444463/1‑150 (MQ=255)

                                                                     aGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAg                                                                                  >  1:4976107/1‑150 (MQ=255)

                                                                           cAACTTAGAAAAACCCCGGTCCATCAGGACGGGGTTTTTTGCATATCAATTATGTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                                                                            <  2:4436687/150‑1 (MQ=255)

                                                                           cAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                                                                            <  2:547366/150‑1 (MQ=255)

                                                                            aaaTTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAgg                                                                           <  2:7707267/150‑1 (MQ=255)

                                                                             aaTTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGga                                                                          <  1:5388403/150‑1 (MQ=255)

                                                                             aaTCAGAAAAACCCCCCTCCAACAGGCCGGCGATTTTTGCATAACAATTATTTGAATTATGAAGGCAACCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTTTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGga                                                                          <  2:5901502/150‑1 (MQ=255)

                                                                               ttAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGgaga                                                                        >  1:4603818/1‑150 (MQ=255)

                                                                               ttAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGgaga                                                                        >  1:2941980/1‑150 (MQ=255)

                                                                                  gAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc                                                                     >  1:4256829/1‑150 (MQ=255)

                                                                                  aaaaaaccccggccttatggaccggggtttttTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTAAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc                                                                     >  1:4119527‑M2/33‑150 (MQ=255)

                                                                                    AAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTg                                                                   >  1:7352327‑M1/1‑150 (MQ=255)

                                                                                     aaaCCCCGGACCAAAAGGCCGGGGTTATTTGAATAACAATTATTTGAATGATGAAGGGAATCTCATGTCAGTTCTGAATATCCAAAAACGTCGGAATCAAATTACAGTTAGCTATCTTGAAAGCAAACGAGAAGTGTCAGTAGACGCAGc                                                                  >  2:1978353/1‑150 (MQ=255)

                                                                                     AAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGGAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGc                                                                  >  1:188258‑M1/1‑150 (MQ=255)

                                                                                      AACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGcc                                                                 >  2:143555‑M1/1‑150 (MQ=255)

                                                                                       ACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCt                                                                >  1:5939284‑M1/1‑150 (MQ=255)

                                                                                        CCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCtt                                                               >  2:2050308‑M1/1‑150 (MQ=255)

                                                                                        CCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCtt                                                               <  2:2797882‑M1/150‑1 (MQ=255)

                                                                                         CCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGGCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTc                                                              <  2:602557‑M1/150‑1 (MQ=255)

                                                                                          CCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCa                                                             >  2:5212009‑M1/1‑150 (MQ=255)

                                                                                          CCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCa                                                             >  2:1618889‑M1/1‑150 (MQ=255)

                                                                                                 aTAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCa                                                      >  2:6914095/1‑150 (MQ=255)

                                                                                                   aaGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAgc                                                    >  1:8097473/1‑150 (MQ=255)

                                                                                                      gCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGcc                                                 <  1:2135405/150‑1 (MQ=255)

                                                                                                        cggggtttttTGCATATCAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTg                                               <  1:725724‑M2/140‑1 (MQ=255)

                                                                                                        cggggtttttTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTg                                               >  1:2070435‑M2/11‑150 (MQ=255)

                                                                                                        cggggtttttTGCATATCAATTATTTGCATGATGAAGGGAAACTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAATACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCCGCGCCTg                                               >  1:6922664‑M2/11‑150 (MQ=255)

                                                                                                         ggggtttttTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGt                                              <  1:4540289‑M2/141‑1 (MQ=255)

                                                                                                              ttttTGCATATCAATTATTTGCATGATGAAGGGAATCTCCTGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATc                                         <  1:3158185‑M2/146‑1 (MQ=255)

                                                                                                                  tGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCa                                     <  2:6284272/150‑1 (MQ=255)

                                                                                                                        tCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCt                               <  2:4155121/150‑1 (MQ=255)

                                                                                                                        tCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCt                               <  2:7323171/150‑1 (MQ=255)

                                                                                                                        tCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCt                               >  2:7128057/1‑150 (MQ=255)

                                                                                                                         cAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCtt                              <  1:5807829/150‑1 (MQ=255)

                                                                                                                          aaTTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCttt                             >  2:6679451/1‑150 (MQ=255)

                                                                                                                          aaTTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCttt                             <  1:1720969/150‑1 (MQ=255)

                                                                                                                          aaTTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCttt                             >  2:792067/1‑150 (MQ=255)

                                                                                                                              attTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTttgttg                         <  1:8011301/150‑1 (MQ=255)

                                                                                                                                ttGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCa                       <  1:3566910/150‑1 (MQ=255)

                                                                                                                                ttGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATGACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCa                       <  2:4763871/150‑1 (MQ=255)

                                                                                                                                   cATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGaa                    >  1:4231373/1‑150 (MQ=255)

                                                                                                                                   cATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGaa                    >  2:293118/1‑150 (MQ=255)

                                                                                                                                      gatgaAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGaaaaa                 >  1:5159437/1‑150 (MQ=255)

                                                                                                                                        tgaAGGGAATCTCAAGTCAGATCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGAGATCTTGAAAGCAAAAGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGAAGAAAATGt               >  2:7571527/1‑150 (MQ=255)

                                                                                                                                         gaAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGtt              >  1:7506384/1‑150 (MQ=255)

                                                                                                                                           aGGGAATCTCCTGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCt            <  1:4830575/150‑1 (MQ=255)

                                                                                                                                           aGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCt            >  1:7772318/1‑150 (MQ=255)

                                                                                                                                            gggAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTg           >  1:2607036/1‑150 (MQ=255)

                                                                                                                                            gggAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTg           <  2:4238762/150‑1 (MQ=255)

                                                                                                                                               aaTCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCaa        >  1:6826681/1‑150 (MQ=255)

                                                                                                                                                 tctcATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCTAGa      >  1:2444322/1‑150 (MQ=255)

                                                                                                                                                 tctcATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGa      >  2:2286731/1‑150 (MQ=255)

                                                                                                                                                 tctcATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGa      >  1:888696/1‑150 (MQ=255)

                                                                                                                                                  ctcATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAACAAGATCttgcagaa     <  2:678259/150‑8 (MQ=255)

                                                                                                                                                  ctcATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGAt     >  1:1015478/1‑150 (MQ=255)

                                                                                                                                                   tcATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATc    >  2:3613743/1‑150 (MQ=255)

ACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACCAGCGCCTGTTAATCGCCAATTTCTTTGTTGCAGAAAAAGTTCTGCAAGATCTT  >  NC_000913/4213085‑4213383

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